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Abstract 


Multivariable,  output  feedback  digital  control 
laws  are  designed  for  the  short  take-off  and  landing  F-15 
aircraft  in  the  landing  configuration.  The  design  is  based 
on  the  methods  developed  by  Professor  Brian  Porter  of  the 
University  of  Salford,  England,  and  was  accomplished  using 
a  computer-aided  design  and  simulation  program  called  MULTI. 

The  STOL  F-15  landing  configuration  includes  canards 
and  reversable  thrust  in  addition  to  the  conventional  F-15 
control  surfaces.  The  additional  controls  allow  decoupling 
of  the  output  variables  in  the  longitudinal  plane.  Longi¬ 
tudinal  aircraft  dynamics  are’ derived  from  data  provided 
by  McDonnell-Douglas ,  the  prime  contractor  for  the  STOL 
F-15,  and  are  presented  in  linearized  state  space  form  for 
the  design  procedure. 

Control  laws  are  developed  to  stabilize  the  air¬ 
craft  to  perform  longitudinal  landing  maneuvers  (flight 
path  control  and  flare)  at  six  flight  conditions.  The 
design  encompasses  actuator  dynamics,  computational  delay, 
sensor  dynamics,  sensor  noise,  and  plant  nonlinearities. 
Proportional  plus  integral  controller  designs  for  each  of 
the  flight  conditions  demonstrate  good  time  response  char¬ 
acteristics.  The  designs  of  two  of  the  flight  conditions 
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are  sufficiently  insensitive  to  plant  variations  to  be 
used  at  all  but  one  of  the  remaining  flight  conditions. 

The  technique  of  multivariable  output  feedback, 
through  the  use  of  the , program  MULTI,  is  shown  to  provide 
good  robust  designs  for  the  STOL  F-15.  Additional  areas 
of  research  on  this  aircraft  are  discussed  as  well  as  sug 
gested  enhancements  to  the  MULTI  program. 
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MULTIVARIABLE  OUTPUT  CONTROL  LAW  DESIGN  FOR  THE 
STOL  F-15  IN  LANDING  CONFIGURATION 

I .  Introduction 

1. 1  Background 

Short  takeoff  and  landing  (STOL)  and  thrust  vector¬ 
ing  technology  are  not  new  concepts  in  the  aerospace  indus¬ 
try.  STOL  techniques  such  as  high  lift  devices  or  wing 
flaps  have  long  been  employed  on  aircraft  of  all  types  and 
sizes.  In  the  past  few  decades  even  more  sophisticated 
STOL  devices  have  been  developed  for  special  purpose  air¬ 
craft  both  in  the  civilian  and  military  aviation  field. 
Vectored  thrust  and,  more  commonly,  thrust  reversing  are 
routinely  used  on  commercial  airliners  to  facilitate  shorter 
landing  rolls.  Thrust  reversing  propellers  on  the  Lockheed 
C-130  Hercules  provide  excellent  short  field  characteris¬ 
tics  for  that  nineteen-fifties  vintage  aircraft.  It  was 
not  until  the  late  sixties  and  early  seventies,  however, 
that  these  techniques  proved  successful  on  a  fighter  air¬ 
craft.  The  British  Aerospace  Harrier  aircraft  design 
allowed  not  only  short  field  operations,  but  also  vertical 
takeoff  capability.  That  design  features  moveable  engine 
nozzles  that  can  be  pointed  backward  for  normal  takeoff  and 
flight,  downward  for  hover  or  vertical  flight,  and  even 


forward  for  rapid  deceleration.  Deployed  by  the  United 
Kingdom  to  the  Falkland  Island  conflict  in  1982,  the 
Harrier  proved  itself  in  aerial  combat  against  the  Argentine 
air  forces.  The  maneuverability  and  operational  flexibil¬ 
ity  of  operating  from  short  Icinding  facilities  proved 
invaluable  to  the  success  of  the  air  war  in  the  Falklands 
(4:38).  Despite  the  successes  of  the  Harrier,  it  is  fairly 
slow  for  a  fighter  and  lacks  the  range,  payload,  and  sophis¬ 
ticated  avionics  of  modern  fighter  aircraft.  The  United 
States  Air  Force  recognizes  the  advantages  of  STOL/ thrust 
vectoring  technology  and  in  September  1983  issued  a  "request 
for  proposal"  for  advanced  development  of  a  STOL  demon¬ 
strator  fighter  aircraft.  Plans  call  for  an  aircraft  cap¬ 
able  of  takeoff  and  landing  in  1500  feet  or  less  in  all 
weather  conditions.  Furthermore,  the  Air  Force  seeks  an 
aircraft  with  supersonic  capability  and  advanced  maneuver¬ 
ing  technology,  including  thrust  vectoring  and  integrated 
flight  and  propulsion  controls  (11:30). 

Current  technology  fighters  require  extremely  long 
runway  surfaces,  particularly  for  landing.  The  McDonnell 
Douglas  F-15,  the  premier  air  superiority  aircraft  in  the 
Air  Force  inventory,  requires  8000  feet  to  safely  land  in 
an  all  weather  environment.  Runways  are  easily  targetable 
facilities  and  it  is  merely  prudent  to  assume  that  substan¬ 
tial  battle  damage  will  be  incurred  in  a  conflict.  A  run¬ 
way  requires  little  damage  to  reduce  its  usable  length  to 
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just  a  few  thousand  feet,  rendering  the  fighter  force 
useless  if  on  the  ground,  and  unrecoverable  if  in  the  air. 
Moveable  barrier  cables  similar  to  the  fixed  barrier  cur¬ 
rently  used  by  the  Air  Force  and  Navy  have  been  proposed  as 
a  simple  solution  to  the  landing  problem.  However,  such 
systems  detract  from  an  aircraft's  autonomy  by  increasing 
its  dependence  on  ground  support  equipment.  Reducing  the 
aircraft's  runway  requirements  opens  up  a  number  of  reason¬ 
able  alternatives.  Fighter  aircraft  could  operate  from 
usable  portions  of  battle  damaged  military  runways,  smaller 
civilian  fields,  or  even  stretches  of  unobstructed  highway. 
In  general,  three  factors  drive  the  amount  of  runway 
required  by  an  aircraft:  approach  velocity,  touchdown  dis¬ 
persion  and  braking  capability.  STOL  and  thrust  vectoring 
technology  can  reduce  the  approach  speed  and  provide  more 
precise  flight  path  control.  Most  importantly,  thrust 
reversing  substantially  decreases  the  stopping  distance  of 
a  high  speed  fighter  regardless  of  runway  braking  coeffi¬ 
cients.  The  benefits  extend  to  shorter  takeoff  rolls  as 
well,  although,  as  is  the  case  with  most  high  thrust  air¬ 
craft  like  the  F-15,  takeoff  distance  is  generally  much 
shorter  than  landing. 

The  design  of  an  aircraft  that  has  these  capabili- 


ties  often  requires  destabilizing  aerodynamic  surfaces,  and 
almost  without  exception  a  greater  number  of  controllable 
surfaces.  It  has  been  shown  in  numerous  previous  papers 


and  theses,  as  well  as  in  experimental  aircraft  like  the 
General  Dynamics  AFTI/F-16  and  the  Grumman  X-29,  that 
on-board  flight  control  computers  are  capable  of  compen¬ 
sating  for  instability  and  efficiently  controlling  addi¬ 
tional  control  surfaces  to  achieve  precise  decoupled  con¬ 
trol.  This  thesis  investigates  the  design  of  flight  con¬ 
trol  laws  for  the  STOL/F-15  as  a  preliminary  determination 
of  the  ability  to  achieve  precise,  decoupled,  low  speed 
flight  path  control. 

1 . 2  Problem 

The  objective  of  this  thesis  is  to  design  longi¬ 
tudinal  control  laws  for  the  STOL/F-15,  based  on  a  well 
developed  mathematical  model,  that  provide  stabilization  as 
well  as  precise  decoupled  flight  path  control  at  landing 
airspeed.  Using  the  multivariable  output  feedback  control 
law  design  techniques  of  Professor  Brian  Porter  of  the 
University  of  Salford,  England,  control  laws  are  developed 
to  perform  longitudinal  landing  maneuvers  at  six  flight 
conditions.  A  build-up  design  approach  is  used,  starting 
with  the  basic  aircraft  and  then  including  in  succession 
the  actuator  dynamics,  computational  time  delay,  sensor 
dynamics,  plant  nonlinearities,  and  sensor  noise. 

1 . 3  Scope 

This  thesis  accomplishes  the  following  objectives 
toward  the  ultimate  goal  of  practical  control  of  the  STOL/ 


1.  Successful  control  of  the  linear  aircraft  model 
including  actuator  dynamics,  sensor  dynamics,  and  computa¬ 
tional  time  delay. 

2.  Identification,  simulation,  and  successful 
compensation  of  a  specific  nonlinearity  arising  from  large 
control  surface  deflections. 

3.  Simulation  of  output  measurement  noise  and  its 
effects  on  the  performance  of  the  closed  loop  performance 
of  the  aircraft. 

4 .  Numerous  enhancements  to  the  computer  aided 
design  and  simulation  program  MULTI  (9) . 

1. 4  Overview 

A  general  description  of  the  aircraft  used  for  this 
study  is  presented  in  Chapter  II,  followed  in  Chapter  III 
by  a  more  detailed  discussion  of  the  mathematical  represen¬ 
tation  and  the  simplifying  assumptions  made  to  obtain  the 
model  used  for  design  and  simulation.  Having  defined 
the  mathematical  model,  the  details  of  the  design  procedure 
are  covered  in  Chapter  IV,  including  some  preliminary 
results  that  demonstrate  the  effects  of  each  of  the  design 
variables.  Chapter  V  contains  the  results  of  applying  the 
design  procedure  of  Chapter  IV  to  the  various  flight  con¬ 
ditions  and  levels  of  model  complexity  presented  in 
Chapter  III,  culminating  in  a  simulation  of  the  aircraft 
with  actuator  dynamics,  sensor  dynamics,  computational 


time  delay,  control  surface  nonlinearity,  and  sensor  noise 
The  chapter  also  presents  a  demonstration  of  controller 
robustness  to  plant  parameter  variations.  Finally, 

Chapter  VI  summarizes  these  results  and  recommends  poten¬ 
tial  improvements  and  recommended  topics  for  future  study. 
Four  appendices  are  included  as  supplementary  material  to 
augment  the  material  presented  in  the  body  of  the  thesis. 
Appendix  A  details  the  revisions  and  additions  made  to 
MULTI  in  the  course  of  the  thesis.  Appendix  B  presents  a 
brief  overview  of  the  theory  behind  the  Porter  method,  and 
Appendices  C  and  D  contain  modeling  data  for  the  actuators 
sensors  and  aircraft  in  general. 


II.  The  STOL  F-15  Aircraft 


2 . 1  Introduction 

The  STOL  F-15,  illustrated  in  Figure  2.1,  is  a 
technology  demonstrator  aircraft  and  is  the  object  of  a 
program  to  investigate,  develop,  and  validate  several  tech 
nological  developments  related  to  STOL  capability  for 
fighter  aircraft.  There  are  two  principal  objectives  of 
this  program:  to  demonstrate  the  use  of  two-dimensional 
thrust  vectoring/reversn.g  nozzles,  integrated  flight/ 
propulsion  control,  rough  field  STOL  Landing  Gear,  and 
an  advanced  pilot/ vehicle  interface;  and  to  provide  design 
options  for  future  fighter  aircraft,  specifically  the 
Advanced  Tactical  Fighter  (ATF) .  In  pursuit  of  these  objec 
tives  the  McDonnell-Douglas  Aircraft  Corporation  (MCAIR) 
developed  the  demonstrator  aircraft  with  the  following 
design  features: 

1.  Additional  control  surfaces  (canards) 

2.  Two-dimensional  engine  nozzles 

3.  Thrust  reversing  vanes 

4.  Improved  landing  gear 

5.  Advanced  avionics  and  cockpit  instrumentation 
This  chapter  is  devoted  to  describing  the  basic  airframe 
and  the  first  three  items  on  this  list  in  detail  and 


discussing  their  effect  on  the  control  of  the  aircraft. 
Since  this  thesis  is  an  investigation  of  the  physical  con¬ 
trol  of  this  vehicle  in  the  approach  phase  of  flight,  the 
improved  avionics  and  landing  gear  are  not  relevant  topics 
for  further  discussion. 

2 . 2  General  Description 

The  original  F-15  aircraft,  a  modern  and  sophisti¬ 
cated  airframe,  has  been  proven  by  more  than  a  decade  of 
safe  and  reliable  peacetime  cind  combat  service.  MCAIR  took 
advantage  of  this  by  using  the  same  basic  airframe  as  a 
baseline  for  the  demonstrator  aircraft  design.  The  out¬ 
ward  appearance  and  dimensions  of  the  STOL  aircraft 
(Figure  2.2)  are  very  similar  to  that  of  the  F-15B  (tandem 
seat  version)  except  for  the  addition  of  the  canards  on 
the  engine  inlets.  Structurally,  the  airframe  is  essen¬ 
tially  unchanged  other  than  the  canard  torque  shaft  and 
fairings  for  both  the  canard  and  nozzles.  Internally,  the 
flight  control  system  has  been  replaced  by  an  integrated 
digital  fly  by  wire  control  system.  The  STOL  F-15  features 
all  of  the  same  control  surfaces  as  the  original  aircraft. 
These  include  ailerons,  trailing-edge  flaps,  horizontal 
stabilators,  two  rudders,  and  a  speed  brake.  In  general, 
the  basic  F-15  is  a  two-engine  fighter  aircraft  charac¬ 
terized  by  an  unusually  high  thrust  to  weight  ratio 
(greater  than  1.0  in  certain  conditions)  and  exceptional 


maneuverability.  Although  untested  in  flight,  MCAIR 
anticipates  substantial  improvement  in  the  STOL  F-15  in 
both  performance  and  maneuverability  based  on  their  calcu¬ 
lations  and  wind  tunnel  tests. 

2 . 3  Canards 

Rather  than  design  and  manufacture  a  canard,  MCAIR 
chose  to  adapt  the  stabilator  of  the  F-18  aircraft,  an 
airplane  smaller  than  the  F-15  but  of  comparable  performanc 
and  sophistication.  The  canards  are  located  on  the  engine 
intakes  (Figure  2.1)  just  aft  of  the  variable  inlets.  As 
is  often  the  case  with  the  addition  of  canards,  the 
destabilizing  effect  of  the  surfaces  forward  of  the  center 
of  gravity  results  in  static  instability  in  some  flight 
conditions.  Although  this  situation  requires  active  con¬ 
trol  of  the  aircraft,  zero  or  negative  static  stability  is 
often  desirable  to  improve  the  maneuvering  capability  of 
the  aircraft.  Instability  is  of  no  advantage  in  the  land¬ 
ing  phase  of  flight  and  of  course  requires  effective, 
reliable  stability  augmentation.  Even  though  the  canard 
is  a  relatively  fast  surface,  MCAIR  plans  to  schedule  the 
canard  with  angle  of  attack,  accomplishing  the  bulk  of  the 
stability  augmentation  through  the  use  of  the  stabilators. 
For  the  purposes  of  this  thesis,  however,  the  canard  will 
be  treated  as  a  fully  controllable  surface.  Like  the 
stabilators,  the  left  and  right  canards  can  be  actuated 


either  together  or  independently,  allowing  their  use  in 
both  longitudinal  and  lateral  control  of  the  aircraft. 

The  canards  are  installed  with  fifteen  degrees  of  dihedral, 
enhancing  lateral  stability  as  well  as  affording  addi¬ 
tional  control  in  the  lateral  mode.  Deflection  of  the 
canard  is  limited  to  +15  degrees  (leading  edge  up) ,  and 
-35  degrees,  at  a  rate  of  23  degrees  per  second. 

2 . 4  Two-Dimensional  Nozzle  and 

Thrust  Reversing  Vanes 

The  conventional  nozzles  of  the  F-15  have  been 
replaced  with  the  two-dimensional  nozzles  depicted  in 
Figures  2.3,  2.4,  and  2.5.  Each  nozzle  features  four 
flaps  and  ten  vanes  to  control  the  thrust  of  the  jet  engine 
exhaust.  As  shown  in  these  figures,  the  nozzles  have  two 
modes  of  operation,  primary  jet  vectoring  and  rotating 
vane  vectoring.  In  the  primary  jet  vectoring  mode  the  ten 
rotating  vanes  remain  closed  and  the  direction  of  the  pri¬ 
mary  jet  is  turned  as  much  as  20  degrees  with  four  hydrauli¬ 
cally  actuated  flaps.  The  four  flaps  (two  convergent  and 
two  divergent)  produce  the  desired  force  and  moment  while 
simultaneously  maintaining  the  required  pressure  gradient 
and  area  ratio  for  the  nozzle.  In  this  mode,  the  principal 
means  for  controlling  the  magnitude  of  the  thrust  is 
through  fuel  flow  to  the  engine  as  in  a  conventional  air¬ 
craft.  The  two  nozzles  (left  and  right)  can  be  controlled 
symmetrically  and/or  differentially,  influencing  both  the 
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longitudinal  and  lateral  motion  of  the  aircraft.  In  the 
rotating  vane  vectoring  mode,  the  ducting  for  the  primary 
jet  exhaust  closes  down  entirely,  diverting  the  flow 
through  the  rotating  vanes.  The  vanes  are  also  hydrauli¬ 
cally  controlled  and  their  variable  deflection  angle  pro¬ 
duces  controllable  forces  and  moments.  Of  the  ten  vanes 
(five  each  top  and  bottom)  eight  are  dedicated  to  vector¬ 
ing  the  thrust  in  the  desired  direction  with  deflection 
limits  of  45  degrees  to  135  degrees  at  a  rate  of  180  degrees 
per  second.  The  remaining  two  (one  each,  top  and  bottom) 
are  deflected  independently  to  maintain  the  equivalent 
nozzle  throat  area  at  the  optimal  value.  The  top  and  bottom 
vanes  are  controlled  independently  except  for  the  two  vanes 
reserved  for  throat  area  control.  In  this  mode,  the 
resultant  force  produced  can  be  controlled  without  changing 
the  RPM  of  the  engines,  reducing  wear  and  tear  on  the  engine 
and  allowing  rapid  transition  to  a  full  thrust  landing 
abort.  Also,  since  the  vanes  are  significantly  faster 
than  the  engine  response  time,  much  more  precise  thrust  con¬ 
trol  is  possible.  Like  the  nozzles,  the  left  and  right 
engine  vanes  are  independent  and  contribute  to  both  the 
longitudinal  and  lateral  motion.  As  a  result,  the  avail¬ 
able  thrust  ranges  from  70  percent  of  its  maximum  military 
thrust  aft  (the  afterburner  is  disabled  in  this  mode)  to 
the  same  amount  of  thrust  reversal,  with  an  infinite  range 
of  control  in  between.  The  precise  thrust  control  and 


thrust  reversal  of  this  mode  are  well  suited  to  the  require 
ments  of  landing,  so  it  is  the  principal  mode  of  operation 
for  approach  and  landing  flight  conditions.  Since  this 
thesis  is  limited  to  the  study  of  the  landing  performance 
of  this  aircraft,  it  is  always  assumed  that  the  aircraft  is 
in  the  rotating  vane  vectoring  mode. 

2 .5  Summary 

The  STOL  F-15  is  a  demonstrator  aircraft  design 
derived  from  an  operational,  high  performance  fighter  air¬ 
craft.  The  addition  of  canards  and  thrust  vectoring/ 
reversing  affords  the  STOL  F-15  enhanced  control  authority 
and  projected  improvements  in  performance.  This  thesis 
is  a  study  of  one  control  law  design  technique  to  make  use 
of  the  STOL  F-15's  unique  capabilities  in  the  landing  con¬ 
figuration. 


III.  Aircraft  Model 


3 . 1  Introduction 

The  Porter  method  of  multivariable  output  feedback 
control  is  founded  in  the  principals  of  linear  control 
theory  and  requires  that  the  system  to  be  controlled  be 
expressed  in  linear  state  space  form.  Like  all  physical 
systems,  airplanes  are  not  linear,  but  in  most  cases  their 
equations  of  motion  are  adequately  approximated  by  a  linear 
system  of  perturbation  equations  around  a  trim  operating 
point.  This  chapter  describes  the  linear  aircraft  model, 
beginning  with  the  fundamental  assumptions  and  developing 
M  the  longitudinal  equations  for  six  flight  conditions.  In 

addition,  the  equations  for  modeling  the  dynamics  of  the 
surface  actuators,  output  sensors  and  noise,  and  computa¬ 
tional  time  delay  are  defined.  Finally,  a  limitation  to 
the  linearized  equations  is  identified,  and  a  nonlinear 
approximation  of  the  solution  is  presented. 

3. 2  Fundamental  Assumptions 

In  the  development  of  linear,  time-invariant  equa¬ 
tions  of  motion  for  aircraft,  a  number  of  assumptions  have 
been  adopted  as  standard  and  are  found  in  nearly  every  text 
that  develops  these  equations.  Following  is  a  summary  of 
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the  results  of  the  rigorous  descriptions  of  these  assump¬ 
tions  found  in  Etkin  (7:121-189)- 


I 

I 


I  IS 

I 


1.  The  surface  of  the  earth  is  flat  and  is  a  sta¬ 
tionary  inertial  reference  frame. 

2 .  The  air  is  stationary  with  respect  to  the 
earth's  surface. 

3.  The  physical  dimensions  of  the  aircraft  do  not 
change  in  time,  neglecting  changes  in  mass,  and  any  bending 
of  the  airframe. 

4.  The  aerodynamic  characteristics  of  the  aircraft 
are  fixed  for  a  given  flight  condition. 

5.  Airflow  variations  that  result  from  the  air¬ 
craft's  maneuvers  occur  instantaneously. 

6.  Perturbation  of  the  aircraft  from  equilibrium 
is  sufficiently  small  to  use  the  first  order  Taylor  series 
approximation  to  the  sine  and  cosine  of  perturbation  angles 

7.  The  motion  is  constrained  to  the  longitudinal 
plane  and  is  assumed  to  be  uncoupled  from  all  lateral 
motion.  This  requires  the  existence  of  a  plane  of  symmetry 
and  no  gyroscopic  effects. 

3 . 3  Linearized  Longitudinal  Equations 

The  application  of  the  preceding  assumptions  to  the 
generalized  equations  of  motion  result  in  a  system  of  four 
linear,  time-invariant  differential  equations  (7:163). 

These  equations  are  expressed  in  terms  of  the  dimensional 
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aerodynamic  stability  derivatives  in  the  stability  axis 
system  (X  axis  aligned  with  relative  wind  at  equilibrium) . 
Although  these  equations  are  useful  for  determining  the 
effects  of  specific  derivatives  on  the  motion  of  the 
aircraft,  it  is  more  convenient  for  controller  design  pur¬ 
poses  to  express  the  relationships  in  the  principal  (body) 
axis  system  (Figure  3.1)  and  in  dimensional  state  space 
for  (2:37).  The  STOL  F-15  data  provided  by  MCAIR  are 
expressed  in  nondimensional  stability  axis  coefficient 
form,  which  necessitates  computation  of  the  matrix  elements 
and  rotation  through  the  angle  of  attack  to  obtain  a  body 
axis  system  of  linear  equations.  These  computations  and 
results  are  contained  in  Appendix  D.  The  sign  conventions 
for  control  surface  deflections  are  also  shown  in  Figure  3.1 

3 . 4  Flight  Conditions 

One  of  the  drawbacks  of  linearized  aircraft  equa¬ 
tions  is  that  their  validity  is  entirely  dependent  on  the 
assumption  that  a  number  of  quantities  remain  constant, 
even  though  they  are  in  fact  variables.  The  process  of 
linearization  is  therefore  valid  at  only  one  particular 
design  condition  and  must  be  reaccomplished  whenever  the 
conditions  are  changed.  This  suggests  that  the  equations 
may  be  time  variant,  which  is  true,  but  in  the  case  of  the 
conditions  of  flight  it  is  assumed  that  the  quantities 
change  slowly  enough  in  time  that  they  are  constant  over 


the  time  period  of  interest.  To  perform  a  thorough  con¬ 
troller  design  it  is  necessary  to  determine  the  extent  of 
the  variability  of  these  pseudo-constants  and  linearize  at 
flight  conditions  throughout  the  parameter  space.  For  this 
thesis,  six  landing  flight  conditions,  listed  in  Table  3.1, 
are  chosen  from  the  available  data.  All  of  these  flight 
conditions  are  in  landing  configuration  (gear  dovm,  flaps 
20  degrees) .  With  1500  lbs  of  fuel  remaining,  the  aircraft 
weighs  33576  lbs,  a  condition  approximating  a  landing  at 
the  end  of  a  mission.  At  this  weight  and  at  sea  level  on 
a  standard  day,  the  recommended  landing  speed  is  200  ft/sec 
(condition  1) ,  20  percent  higher  than  the  stall  speed  of  168 
ft/sec  (condition  2) .  With  11435  lbs  internal  fuel,  the 
gross  weight  is  43511  lbs,  a  simulated  heavy-weight  landing. 


TABLE  3.1 

FLIGHT  CONDITIONS 

Flight 

Altitude 

Aircraft  Weight 

True  Velocity 

Condition 

(Ft) 

(lbs) 

(Ft/Sec) 

1 

0 

33576 

200 

2 

0 

33576 

168 

3 

0 

43511 

200 

4 

0 

43511 

304 

5 

10000 

43511 

304 

10000 


33576 


200 


The  velocities  chosen  in  flight  conditions  3,  4,  5,  and  6 
are  not  chosen  for  any  particular  significance  other  than 
that  the  data  are  available  and  represent  reasonable  land¬ 
ing  conditions. 

3 . 5  Representative  Plant 

Introduction.  In  the  course  of  this  thesis,  sample 
calculations  and  examples  are  included  to  illustrate  the 
concepts  and  techniques  used  throughout  the  design  process. 
For  consistency  and  brevity,  one  plant  model  is  chosen  as 
a  representative  plant  for  this  purpose.  The  plant  model 
of  flight  condition  1  is  chosen  because  it  represents  the 
conditions  most  commonly  encountered  in  an  approach  and 
landing  situation.  The  remainder  of  this  section  details 
the  development  of  each  of  the  plant  models  using  the 
representative  plant  as  an  excimple. 

Plant  Matrix.  The  plant  matrix,  also  referred  to 
as  the  "A"  matrix,  is  calculated  as  described  in  Section  3.3 
of  Appendix  D.  The  result  for  flight  condition  1  is  shown 
in  Figure  3.2.  The  states  in  this  model,  as  with  all  the 
models  of  this  thesis  are  pitch  angle  (6),  X-axis  component 
of  velocity  (u) ,  pitch  rate  (q) ,  and  the  angle  of  attack 
(a) .  Note  that  the  first  state  (6)  is  chosen  such  that  the 
kinematic  equation  is  at  the  top  of  the  A  matrix.  This  is 
a  requirement  of  the  computer  program  MULTI. 


X  =  Ax  +  Bu 


y  =  Cx 


A  = 


0 

-31.54 

0 

-0.03232 


0 

1 

0 

-0.06909 

-40.16 

0.3350 

-3.603x10"'^ 

-0.9912 

1.366 

0.9925x10"^ 

0.9796 

-0.6392 

B  = 


0 

-2.546 

0.8407 

-0.02060 


0 

-3.237 

-1.578 

-0.07660 


0 

-21.80 

-0.0250 

0 


C  = 


0 

0 

1 


1 

0 

0 


0 

0 

0 


0 

1 

-1 


Fig.  3.2.  Open  Loop  State  Space  Model,  Flight  Condition  1 


Input  Matrix.  The  input  ("B")  matrix  is  also  calcu¬ 
lated  according  to  the  procedure  of  Section  3.3  except 
that  some  additional  manipulation  is  required  to  put  it  in 
a  usable  form.  The  STOL  F-15  has  nine  independently  actu¬ 
ated  longitudinal  control  surfaces.  They  are  the  canard, 
the  stabilator,  the  trailing  edge  flaps,  the  ailerons 
(actuated  together  like  the  flaps) ,  each  of  the  four 
rotating  vanes,  and  the  speed  brake.  Since  it  is  a  funda¬ 
mental  requirement  of  multivariable  control  system  theory 
that  there  be  a  unique  solution  to  the  equations,  the  inputs 
must  be  mathematically  independent  and  there  can  be  no  more 
inputs  than  there  are  independent  states.  Pitch  rate  is 
simply  the  time  derivative  of  pitch  angle,  leaving  only 
three  independent  states.  It  should  also  be  obvious  that 
with  only  two  force  equations  and  a  mc'-'ent  equation  in  the 
longitudinal  plane,  a  fourth  input  would  clearly  be  a  linear 
combination  of  the  other  three.  Without  additional  equa¬ 
tions  only  three  independent  inputs  are  possible.  Optimal 
control  theory  allows  the  definition  of  cost  functions  that 
weight  the  inputs  according  to  user  definable  optimality 
criteria.  These  cost  functions  can  provide  the  necessary 
additional  equations  to  use  extra  inputs;  however,  the 
current  capabilities  of  the  Porter  method  do  not  include 
their  use.  The  first  simplification  is  made  by  eliminating 
the  speed  brake  for  consideration.  Although  the  rotating 
vanes,  when  operating  independently,  provide  significant 


independent  contributions  to  all  three  equations,  their 
utility  is  best  realized  as  an  input  into  the  horizontal 
force  equation.  By  summing  the  contributions  of  all  four 
vanes,  their  combined  effect  can  be  treated  as  a  single 
input.  This  same  mathematical  effect  could  be  achieved 
by  physically  prohibiting  independent  actuation  of  the 
rotating  vanes.  The  two  remaining  excess  inputs  are 
removed  by  summing  the  effects  of  the  flaps,  ailerons,  and 
canards.  To  make  maximum  use  of  each  of  these  surfaces 
they  are  weighted  prior  to  the  summation  such  that  they 
reach  maximum  surface  deflection  simultaneously.  Weighting 
the  surfaces  in  this  manner  is  in  effect  a  type  of  cost 
function  that  is  invariant  and  can  be  evaluated  in  advance. 
The  weighting  process  is  described  in  Appendix  C.  There 
are  now  three  inputs  and  they  are  labeled  in  Figure  3.2  as 
the  combined  canard,  flaps,  and  ailerons  acting  as  one 
equivalent  surface  (6^),  the  stabilator  (6„),  and  the 
rotating  vanes  (6,^,)  .  The  "B"  matrix  is  composed  of  partial 
derivatives  of  the  state  equations  with  respect  to  control 
surface  deflections,  and  as  such  are  accurate  only  for 
small  deflections  of  the  surface  away  from  the  equilibrium 
position.  In  certain  circumstances  it  is  necessary  to 
modify  the  "B"  matrix  as  control  surface  deflections  become 
large  to  account  for  the  inaccuracies  of  linearization. 
Paragraph  3.9  of  this  chapter  describes  these  circumstances 
and  the  modifications  made  for  this  thesis. 


Output  Matrix .  The  output  matrix  ("C")  combines 
and  weights  the  states  to  obtain  the  output  variables  which 
are  to  be  controlled.  With  three  independent  inputs,  three 
outputs  can  be  controlled.  Since  any  three  independent 
outputs  can  be  used,  the  designer  must  choose  the  desired 
outputs  based  on  the  type  of  maneuver  to  be  accomplished 
and  the  accessibility  of  the  states  that  must  be  measured. 
The  key  objective  of  a  landing  maneuver  is  to  control  the 
flight  path  of  the  aircraft.  Since  the  flight  path  angle 
is  simply  the  difference  of  the  pitch  angle  and  the  angle 
of  attack,  this  is  a  readily  measured  quantity  as  well. 
Control  of  velocity  and  angle  of  attack  is  also  important 
since  landing  is  usually  at  critically  slow  airspeed  and 
maximum  angle  of  attack.  Pitch  angle,  although  important 
as  visual  feedback  to  the  pilot,  is  subordinate  to  angle 
of  attack  in  this  flight  condition  and  must  be  controlled 
indirectly  using  flight  path  angle  and  angle  of  attack. 

The  result  of  these  considerations  is  the  "C"  matrix  of 
Figure  3.2  and  is  the  same  for  all  of  the  flight  conditions 
of  this  thesis.  The  variables  are  velocity  (u) ,  angle  of 
attack  (a) ,  and  flight  path  angle  (y) • 

3 . 6  Ac tuator  Dynamics 

The  equations  of  Figure  3.2  relate  the  motion  of 
the  aircraft  to  the  deflections  of  the  control  surfaces. 
Since  the  control  input  quantities  are  signals  generated 


from  the  measurement  of  the  outputs,  the  dynamic  response 
of  the  actuators  (converting  electrical  signals  to  surface 
movement)  must  be  considered.  The  data  provided  by  MCAIR 
(Appendix  C)  specifies  the  actuator  dynamics  for  each 
surface  as  LaPlace  domain  transfer  functions,  most  of  which 
have  third  order  characteristic  equations.  The  program 
MULTI  allows  entry  of  actuator  transfer  functions  without 
augmenting  the  plant  matrix  but  is  limited  to  second  order 
characteristic  equations  and  zero  order  numerators.  Since 
the  surfaces  have  been  combined  as  described  in  paragraph 
3.5,  the  actuator  data  provided  by  MCAIR  must  be  altered 
to  satisfy  these  constraints.  Appendix  C  demonstrates  that 
in  all  cases  the  dynamics  are  well  modeled  by  second  order 
transfer  functions  with  no  zeros,  even  after  combining 
surfaces.  Figure  3.3  contains  the  results  of  these  calcu¬ 
lations  and  illustrates  the  placement  of  the  actuator 
dynamics  in  block  diagram  form.  The  actuator  dynamics  are 
constant  and  independent  of  flight  condition.  Note  that 
for  each  actuator  a  variable  gain  is  included  for  potential 
design  enhancement.  The  use  of  these  gains  is  demonstrated 
in  the  following  chapter. 

3.7  Computational  Time  Delay 

The  program  MULTI  features  an  option  of  including 


a  delay  of  one  or  more  sample  periods  to  simulate  the  delay 
in  computing  the  control  inputs  from  the  measurements  of 


Actuator  Dynamics 


Fig.  3.3.  Open  Loop  Block  Diagram 


the  outputs.  The  computational  time  delay  for  this  thesis 
is  chosen  to  be  one  sample  period  (0.025  sec)  and  its 
incorporation  into  the  controller  design  is  discussed  in 
Chapter  IV. 

3.8  Sensor  Dynamics  and  Noise 

In  the  design  of  aircraft  controllers,  output  quan¬ 
tities  being  fed  back  must  be  measured  and  as  a  result  the 
feedback  signals  are  subject  to  the  limitations  of  imperfect 
measurement.  Two  of  these  limitations  are  considered  in 
this  thesis;  sensor  dynamics  and  measurement  noise.  Like 
any  physical  system,  sensors  are  incapable  of  responding 
instantaneously  to  an  input  (the  quantity  being  measured) 
but  their  response  is  usually  well  modeled  with  linear 
transfer  functions  similar  to  actuator  dynamics.  As  with 
the  actuator  dynamics,  transfer  functions  for  the  sensors 
are  entered  into  the  MULTI  program  independent  of  the 
plant,  input,  or  output  matrices  and  are  also  limited  to 
second  order  dynamics  with  no  zeros.  The  data  provided  by 
MCAIR  again  must  be  altered  (Appendix  C)  to  conform  to  these 
constraints  and  the  results  are  shown  in  Figure  3.3. 
Measurement  noise  arises  from  a  variety  of  sources,  each 
with  its  own  stochastic  characteristics.  The  simplest 
model,  and  the  one  used  in  this  thesis,  is  a  zero-mean, 
white,  Gaussian  noise  model.  Each  measured  output  is  cor¬ 
rupted  with  independent  noise,  added  immediately  after  the 


sensor  dynamics,  resulting  in  a  vector  of  measured  outputs 
as  shown  in  Figure  3.3.  As  with  the  actuators,  the  sensor 
dynamics  and  noise  are  constant  over  all  flight  conditions 


3.9  Control  Surface  Nonlinearity 

In  the  derivation  of  the  linear  perturbation  equa¬ 
tions  for  the  aircraft,  it  is  assumed  that  the  deviations 
from  equilibrixam  are  small  enough  for  the  forces  and 
moments  generated  to  be  linear  functions  of  the  deviation. 
For  most  of  the  parameters  the  deviations  can  be  quite 
large  (as  much  as  20  degrees)  before  significant  error 
results.  One  particular  exception  is  the  drag  (force  in 
the  longitudinal  direction)  resulting  from  changes  in  the 
angle  of  attack  on  the  various  aerodynamic  surfaces.  Drag 
variation  with  angle  of  attack  is  usually  dominated  by  the 
induced  drag,  particularly  at  high  angles  of  attack 
(5:149).  Induced  drag  varies  with  the  square  of  the 
lift  coefficient,  generally  considered  a  linear  function 
of  the  angle  of  attack.  The  change  in  drag  of  a  lifting 
surface  at  high  angles  of  attack  (like  the  wing  of  an  air¬ 
craft  in  a  landing  situation)  is  therefore  a  parabolic 
function  of  the  angle  of  attack.  Since  the  angle  of  attack 
of  the  wing  and  fuselage  of  the  aircraft  varies  at  most 
only  a  few  degrees  from  the  equilibrium  value  (see  Chapter 
V)  it  is  reasonable  to  treat  the  wing/body  drag  as  a  linear 
function  of  angle  of  attack.  At  angles  of  attack  near  zero 


(like  control  surfaces  nearly  aligned  with  the  relative 
wind)  the  slope  of  the  function  is  rapidly  changing  and  the 
function  is  no  longer  accurately  modeled  as  linear  unless 
only  a  very  small  range  of  values  is  considered.  The  con¬ 
trol  surfaces,  however,  can  vary  as  much  as  50  degrees  in 
some  cases.  Of  particular  importance  is  that  the  minimum 
value  of  induced  drag  occurs  at  zero  angle  of  attack 
(symmetric  airfoil)  and  in  no  case  is  the  drag  less  than 
zero  (it  is  important  to  note  that  this  discussion  is  deal¬ 
ing  with  the  absolute  drag  and  angle  of  attack,  not  per¬ 
turbation  values) .  A  linear  model  results  in  negative  drag 
(thrust)  at  negative  angles  of  attack,  a  situation  that 
grossly  misrepresents  the  physical  system.  This  phenomenon 
requires  an  accurate  model  to  have  "B"  matrix  elements  that 
are  a  function  of  each  surface's  angle  of  attack.  From  a 
design  perspective  a  time  dependent  B  matrix  would  be 
impractical  since  the  Porter  method  is  dependent  upon 
linear  constant  coefficient  equations.  However,  it  is 
possible  to  recreate  this  phenomenon  during  the  simulation 
of  the  aircraft's  performance  with  a  controller  designed 
using  linear  methods.  If  the  data  were  available,  one 
method  of  simulating  this  effect  would  be  to  have  a  table 
of  values  for  the  drag  coefficients  at  various  angles  of 
attack  for  each  of  the  aerodynamic  control  surfaces.  The 
simulation  program  could  then  use  the  appropriate  entry  in 


such  a  table  as  the  current  value  of  that  element  of  the 


"B"  matrix.  Since  the  data  are  not  available  from  MCAIR 
at  this  time,  an  alternative  solution  is  required.  The 
principal  objection  to  the  linear  model  is  the  apparent 
creation  of  thrust  by  aerodynamic  surfaces.  In  reality, 
an  aerodynamic  surface  has  the  capability  of  reducing  drag 
by  decreasing  its  angle  of  attack  to  zero.  However,  the 
linear  approximation  allows  the  mathematical  equivalent  of 
thrust  to  result  from  large  surface  deflections  in  the 
opposite  direction  of  the  equilibrium  deflection  for  which 
the  linear  derivative  is  accurate.  Observing  that  the 
derivative  of  a  parabolic  drag  function  is  positive  for 
positive  angles  of  attack  and  negative  for  negative  angles 
of  attack  suggests  a  simple  solution  to  this  error.  By 
testing  the  position  of  the  surface  at  each  step  in  the 
simulation,  the  angle  of  attack  is  determined  and  if  neces¬ 
sary  the  sign  of  the  appropriate  control  surface  derivative 
is  reversed  to  insure  that  the  drag  is  always  positive. 

This  results  in  an  absolute  value  function  with  slopes  of 
plus  or  minus  the  equilibrium  control  surface  derivative. 

This  is  a  crude  approximation  to  the  parabolic  function, 
but  precludes  the  erroneous  results  possible  with  the  linear 
model.  Since  drag  is  defined  to  be  parallel  to  the  rela¬ 
tive  wind,  this  force  should  technically  be  rotated  through 
the  angle  of  attack  for  the  aircraft  prior  to  being  included 
in  a  body  axis  representation  of  the  aircraft  equations  of 
motion.  This  would  result  in  a  modification  of  every  element 


of  the  "B"  matrix  and  would  require  knowledge  of  equilibrium 
parameters  not  currently  needed  in  the  simulation.  For 
simplicity,  the  effects  of  this  phenomenon  are  assumed  to 
be  limited  to  the  elements  of  the  "B"  matrix  that  influence 
the  longitudinal  velocity  of  the  aircraft  (u) .  Since  the 
rotation  angle  {a)  is  small,  this  is  a  reasonable  approxi¬ 
mation  of  reality.  This  approach  to  the  problem  is  used 
for  this  thesis  and  its  implementation  is  described  in 
Appendix  A. 


IV.  Design  Procedure 


4 . 1  Introduction 

In  the  preceding  chapters,  the  background  and 
motivation  for  this  thesis  effort  and  the  physical  and 
mathematical  characteristics  of  the  STOL  F-15  are  presented 
as  introductory  material  for  the  design  procedure,  results, 
and  conclusions.  This  chapter  is  the  first  of  three  chap¬ 
ters  that  deal  with  the  key  elements  of  the  thesis  and  is 
devoted  to  the  detailed  description  of  the  procedure  used 
to  design  the  longitudinal  controllers  for  the  STOL  F-15 
during  landing  operations. 

The  techniques  developed  by  Professor  Porter,  out¬ 
lined  in  Appendix  B,  although  founded  in  the  deterministic 
mathematics  of  linear  algebra,  require  considerable  trial 
and  error  and  qualitative  assessment  of  the  system  response 
to  achieve  a  satisfactory  controller  design.  The  determinis¬ 
tic  portion  of  the  design  technique  is  accomplished  through 
use  of  the  computer  program  MULTI.  MULTI  is  an  indispens¬ 
able  tool  in  the  design  process,  but  like  any  computer 
program,  the  quality  of  the  output  is  merely  a  reflection 
of  the  quality  of  the  input.  There  are  a  number  of  design 
parameters  in  the  Porter  method  that  are  left  to  the  user 
to  define,  based  on  theoretical  insight  and  experience 
gained  during  the  trial  and  error  process.  As  the  designer 
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becomes  more  familiar  with  this  method,  the  program  MULTI, 
and  the  characteristics  of  the  system  to  be  controlled, 
trends  can  be  identified  during  design  parameter  selection 
that  lead  to  a  satisfactory  design.  This  chapter  describes 
the  design  procedure  used  for  this  thesis,  and  an  attempt 
is  made  to  justify  design  decisions  based  on  quantitative 
results.  Past  thesis  efforts  have  suggested  that  their 
subjective  design  procedures  can  be  applied  to  a  general 
class  of  problems,  presenting  a  generalized  approach  (2;  5) . 
Many  of  the  observations  contained  in  these  theses  prove 
valuable  in  this  design  effort,  but  in  general  the  design 
procedures  are  found  to  be  unique  to  the  system  being 
controlled  and  the  logical  thought  processes  of  the 
designer.  As  a  result,  the  approach  of  this  thesis  is 
intended  to  be  a  description  of  the  route  followed  in 
obtaining  the  specific  results  of  Chapter  V. 

This  chapter  begins  with  mathematical  discussions 
intended  to  supplement  the  information  in  Appendix  B,  as 
well  as  to  further  justify  the  choice  of  output  vectors. 

The  design  process  is  then  covered  in  detail,  beginning  with 
the  basic  plant  and  expanding  the  design  to  account  for  the 
actuator  dynamics,  computational  delay,  sensor  dynamics, 
and  control  surface  nonlinearities.  Finally,  a  demonstra¬ 
tion  of  the  controller's  sensitivity  to  parameter  variation 
and  output  measurement  noise  is  presented. 


4.2 


Mathematical  Considerations 


Controllability  and  Observability.  It  is  a  funda¬ 
mental  requirement  of  any  method  that  the  system  be  com¬ 
pletely  controllable  and  observable.  Complete  controlla¬ 
bility  implies  that  every  state  in  the  state  vector  can  be 
driven  to  any  finite  value,  in  a  finite  time,  with  a  control 
input  of  finite  magnitude.  Controllability  is  a  function  of 
the  system  plant  and  input  matrix  and  can  be  determined  by 
evaluating  the  system  zeros  (6:443).  If  there  are  no 
input-decoupling  zeros  the  system  is  completely  con¬ 
trollable.  The  program  ZERO  is  used  to  determine  the  sys¬ 
tem  zeros  for  this  thesis  and  reveals  that  there  are  no 
input-decoupling  zeros  for  any  flight  condition.  Thus, 
all  of  the  aircraft  models  are  completely  controllable 
systems.  Complete  observability  requires  that  every  mode 
of  each  state  appear  in  at  least  one  element  of  the  output 
vector.  Like  controllability,  observability  is  also  deter¬ 
mined  from  the  system  zeros  and  the  absence  of  output¬ 
decoupling  zeros  indicates  complete  observability.  Output¬ 
decoupling  zeros  are  a  function  of  the  plant  and  output 
matrices.  Observability  is  influenced  by  the  choice  of 
the  output  vector  elements  and  it  is  imperative  to  estab¬ 
lish  the  existence  of  elements  that  result  in  a  completely 
observable  system.  Fortunately,  ZERO  shows  that  there  are 
no  output-decoupling  zeros  regardless  of  the  choice  of 
physically  meaningful  outputs.  Therefore,  all  of  the  states 


include  each  of  the  system  modes  and  the  output  vector  is 
not  restricted  due  to  observability  considerations.  Porter 
and  Bradshaw  have  shown  that  the  addition  of  proportional 
plus  integral  feedback  control  does  not  change  the  system 
controllability  and  observability  provided  that  some  well 
defined  conditions  are  satisfied  (13). 

Porter  Design  Requirements.  Although  a  number  of 
output  vectors  may  be  chosen  without  causing  an  unobservable 
mode,  the  principals  of  multivariable  output  feedback  con¬ 
trol  in  general  and  specifically  the  Porter  method  place 
additional  restrictions  on  the  desired  output  variables. 
Transmission  zeros,  defined  as  zeros  of  the  equivalent 
transfer  function  representation  that  block  transmission 
of  a  particular  exponential  input,  are  regions  to  vdiich 
some  of  the  slow  roots  of  a  system  migrate  as  the  gain 
approaches  infinity.  Although  transmission  zeros  in  the 
left  hand  "s"  plane  do  not  guarantee  stability  at  all  gains, 
transmission  zeros  in  the  right  half  plane  always  result 
in  a  region  of  a  gain  for  which  the  system  is  unstable. 

Since  determination  of  gain  boundaries  in  multivariable 
control  is  difficult,  it  is  desirable  to  have  all  trans¬ 
mission  zeros  in  the  left  half  plane.  The  existence  of  a 
transmission  zero  at  the  origin,  a  condition  that  results 
from  including  pitch  rate  in  the  output  vector,  is  a  situa¬ 
tion  that  can  complicate  the  search  for  a  satisfactory 
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design  since  this  may  cause  a  region  of  gain  for  which  the 
system  is  unstable.  The  location  of  the  transmission 
zeros  affects  only  the  migration  pattern  of  certain  roots 
as  a  function  of  gain.  At  any  particular  gain  the  roots 
of  the  system  are  fixed  and  unaffected  by  as  asymptotic 
properties  of  the  system.  Therefore,  the  location  of  trans¬ 
mission  zeros  (slow  closed-loop  system  roots)  and  asymptotes 
(fast  roots)  are  significant  only  during  the  design  process 
while  choosing  the  various  gains  (Appendix  B) .  When  the 
number  of  outputs  is  equal  to  the  number  of  inputs  (a 
requirement  of  the  program  MULTI)  the  number  of  transmission 
zeros  is  given  by  (8) : 

#Z^  =  (n  -  m)  -  r  (4-1) 

where 

#Z^  =  number  of  transmission  zeros 
n  =  number  of  states 
m  =  number  of  outputs 
r  =  rank  defect  of  £2^2 

The  use  of  pitch  rate  as  an  output  results  in  one  trans¬ 
mission  zero  (at  the  origin)  because  £282  is  of  full  rank. 
Failure  to  feed  back  pitch  rate  results  in  no  transmission 
zeros.  However,  without  pitch  rate  in  the  output  vector, 
£2^2  has  a  rank  deficiency  of  one,  and  a  measurement  matrix 
is  required  to  use  Porter's  method  of  determining  the 
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controller  gain  matrices.  The  use  of  a  measurement  matrix 
introduces  transmission  zeros  whose  locations  are  dependent 
on  the  matrix  elements  chosen  by  the  designer.  It  is 
interesting  to  note  that  when  the  minimum  measurement 
matrix  is  chosen,  the  result  is  equivalent  to  a  minor  loop 
pitch  rate  feedback  with  the  feedback  gain  determined  by 
the  one  measurement  matrix  entry  required.  The  measurement 
matrix  in  this  case  introduces  one  transmission  zero, 
located  at  the  negative  reciprocal  of  the  pitch  rate  feed¬ 
back  gain.  These  findings  demonstrate  two  important  facts. 
First,  pitch  rate  feedback  is  required  for  the  design  method 
but  pitch  rate  need  not  be  actively  controlled.  Second, 
there  must  be  an  additional  transmission  zero,  but  its  loca¬ 
tion  can  be  chosen  by  the  designer.  Chapter  III  outlines 
the  practical  reasoning  for  choosing  velocity,  angle  of 
attack,  and  flight  path  as  controlled  quantities.  This 
section  also  reveals  that  feeding  pitch  rate  back  in  a  minor 
loop,  rather  than  as  an  output,  affords  additional  control 
over  the  asymptotic  properties  of  the  system.  It  is  impor¬ 
tant  to  realize  that  changing  the  gain  of  the  pitch  rate 
feedback  affects  not  only  system  response  to  high  gain,  but 
also  affects  the  transient  response  of  the  system  at  finite 
gain.  For  this  reason  the  designer  does  not  have  unlimited 
authority  to  arbitrarily  choose  the  location  of  the  trans¬ 


mission  zero. 


4 . 3  Design  Process 

General .  The  generalized  procedure  for  designing 
a  control  law  for  the  most  realistic  STOL  F-15  model 
(actuators,  time  delay,  sensors,  noise,  and  nonlinear  sur¬ 
faces  all  included)  is  to  begin  with  the  basic  model  as  a 
baseline  design  and  to  redesign  for  each  of  these  addi¬ 
tional  factors  one  at  a  time.  The  effects  of  actuators  and 
computational  time  delay  with  respect  to  the  Porter  method 
have  been  investigated  in  the  past  and  are  demonstrated  in 
several  documents  (2;  5).  Therefore,  the  simplest  model 
used  for  design  in  this  thesis  is  that  of  the  plant  with 
actuators  and  computational  time  delay.  However,  the  basic 
plant  is  used  for  computation  of  open  loop  transfer  func¬ 
tions  and  demonstration  of  the  asymptotic  characteristics. 
Unfortunately,  once  the  actuators  and  computational  time 
delay  are  included  in  the  simulation,  the  roots  of  the 
complete  system  are  not  easily  calculated.  From  that  point 
onward,  the  design  process  relies  on  evaluation  of  the  simu 
lated  time  response  of  the  system  to  determine  the  require¬ 
ments  for  changes  in  the  design  parameters. 

Design  Variables.  Following  is  a  list  of  the  vari¬ 
ables  that  are  available  for  assigning  the  asymptotic 
properties  of  the  system.  These  variables  are  defined  in 
Appendix  B. 


1.  sigma  (j;)  matrix  diagonal  elements. 

2.  Epsilon  (e) ,  the  sigma  matrix  multiplier. 

3.  Alpha  (ct)  ,  ratio  of  integral  to  proportional 
control . 

4.  K  ,  measurement  matrix  element  (pitch  rate 
fledback  gain) . 

An  additional  parameter,  stabilator  actuator  gain  (K^) , 
is  used  as  a  design  variable  even  though  actuator  dynamics 
are  actually  plant  parameters. 

Design  for  the  Basic  Plant.  The  state  space  model 
of  Figure  3.2  (flight  condition  1)  results  in  the  following 
LaPlace  transfer  functions: 

CE  =  (s  -  .3468) (s  +  1.998) (s  +  .02428  -  j. 07015) (s  +  .02428  +  j. 07015) 


(4-2) 

u(s)/6^(s)  =  (s  +  31.76) (s  +  .6421  -  j.l493) (s  +  .6421  +  j.l493)  /  CE 

(4-3) 

u(s)/6g(s)  =  (s  +  15.35) (s  +  .03081) (s  -  .06574)  /  CE  (4-4) 

u(s)/5^(s)  =  (s  +  .06900) (s+  .6661  -  j3.068) (s  +  .6661  +  j3.068)  /  CE 

(4-5) 

a(s)/£^(s)  =  (s  -  19.49) (s  +  .9976) (s  +  .5559)  /  CE  (4-6) 

a(s)/6„(s)  =  (s+21.21)(s  -  .03936)  (s  +  .07535)  /  CE  (4-7) 

a(s)/5^(s)  =  (s  +  4.198) (s  -  3.627) (s  +  .06949)  /  CE  (4-8) 

> (s)/6^(s)  =  (s  +  1.979) (s  -  .3659) (s  -  .02897)  /  CE  (4-9) 

Y(s)/6g(s)  =  (s  -  .2343)  (s  -  .03284)  /  CE  (4-10) 

y(s)/6„,(s)  =  (s  +  53.59)  (s  -  .2307)  /  CE  (4-11) 


These  transfer  functions  are  typical  of  a  statically 
unstable  aircraft,  characterized  by  a  lightly  damped  complex 
pair  (phugoid)  and  two  real  roots  (short  period) ,  one  of 
which  is  positive.  Naturally,  the  first  objective  of  the 
controller  is  to  stabilize  the  aircraft,  and  then  the  con¬ 
troller  is  tailored  to  obtain  an  acceptable  response. 

Porter  and  Bradshaw  (13)  have  shown  that  the  roots  of  the 
closed  loop  system  approach  very  predictable  values  as  the 
gain  approaches  infinity.  Specifically,  in  the  four  state, 
three  output  system  of  this  thesis,  there  are  seven  closed 
loop  roots  divided  into  two  categories,  fast  roots  and  slow 
roots.  The  fast  roots  (Z^)  depend  on  the  value  of  gain  and 
are  a  function  of  the  elements  of  diagonal  sigma  matrix 
which  are  selected  by  the  designer. 

Z3  =  {X: 1  XI  -  gE 1  =0}  (4-12) 

where  E  =  FBKO . 

The  slow  roots  of  the  system  consist  of  the  union  of  two 
sets  of  roots  (Z^  Z^) ,  defined  by  Equations  (4-13)  and 
(4, -14)  : 

Z,  =  {X;  (X+K  ■^)  =  0}  (4-13) 

z  q 

Z^  =  {X;]XK0  +  Kl|  =  0}  (4-14) 

Selecting  the  controller  matrices  KO  and  K1  so  that  they 
differ  only  by  a  constant  (a) ,  Equation  (4-14)  reduces  to: 


=  {X:  (X-a)  =  0} 


(4-15 


By  trial  and  error,  and  the  observations  of  Courtheyn 
(5:51)  relating  the  steady  state  transfer  function,  G(0), 
to  the  sigma  matrix,  the  following  design  parameters  are 
found  to  provide  a  satisfactory  time  response  to  a  step 
input.  As  described  in  Appendix  B,  these  parameter  define 
controller  matrices,  and  they  result  in  the  system  illus¬ 
trated  in  Figure  4.1. 

a  =  .01 

e  =  .05 


K  =  .25 

q 

M  =  [  0  K  0 

q 

These  parameters  result  in  the  following  asymptotic  roots: 
Fast  Roots 

=  {-.75g,  -.02g,  -.02g} 

Slow  Roots 

Z.  =  {-4.0} 


Fig.  4,1.  Discrete  Closed  Loop  Block  Diagram 


Although  the  gain  (g)  is  fixed  for  the  STOL  F-15  model  by 
the  reciprocal  of  the  sampling  time  (40.0),  using  gain  as 
a  variable  in  this  discussion  allows  demonstration  of  the 
asymptotic  properties  of  the  system.  Table  4.1  compares 
the  asymptotic  roots  to  the  actual  roots  at  two  different 
gains  to  show  the  progression  with  higher  gain. 

Porter  and  Bradshaw  have  also  shown  that  as  the 
gain  approaches  infinity,  the  input/output  transfer  func¬ 
tion  matrix  asymptotically  approaches  a  diagonal  or  near 
diagonal  form,  depending  on  the  measurement  matrix  and 
choice  of  output  variables  (13) .  if  no  measurement  matrix 
is  used  (pitch  rate  in  the  output  vector)  this  matrix  is 
always  diagonal  and  contains  only  the  fast  roots  of  the 
system,  thereby  exhibiting  increasingly  tight  and  decoupled 
control  with  increasing  gain.  However,  with  a  measurement 
matrix  the  asymptotic  transfer  function  may  not  be  diagonal 
and  always  includes  the  transmission  zero  as  a  mode  that 
dominates  one  of  the  output  responses.  This  results  in  at 
least  one  transfer  function  that  does  not  exhibit  increas¬ 
ingly  tight  control,  and  if  it  is  not  diagonal,  the  off- 
diagonal  terms  define  coupling  terms  that  also  contain  the 
transmission  zero  as  a  mode,  a  situation  that  does  not 
change  with  increasing  gain.  It  is  important  to  note  that 
if  the  transmission  zero  is  not  due  to  a  measurement  matrix 
the  asymptotic  transfer  function  matrix  is  diagonal  and  the 
mode  associated  with  the  transmission  zero  does  not  appear 
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TABLE  4.1 

COMPARISON  OF  ASYMPTOTIC  AND  FINITE  SYSTEM  ROOTS 


Asymptotic  Roots 

-.01  (a) 

-.01  (a) 

-.01  (a) 

-4.0  (Kg~^) 
-.80  (.02g) 
-.80  (.02g) 
-30.0  (.75g) 
-.01  (a) 

-.01  (a) 

-.01  (a) 

-4.0  (K^“^) 
-1.6  (.02g) 
-1.6  (.02g) 
-60.0  (.75g) 


System  Roots 

-.00776  +  .001839j 
-.00776  -  .001839j 
-.01154 
-1.873 

-.6700  +  1.200j 
-.6700  -  1.200j 
-30.06 

-.0120  +  .00436j 
-.0120  -  .00436j 
-.004388 
-2.417 

-1.197  +  1.893j 
-1.197  -  1.893j 


-60.06 


in  the  matrix.  Since  the  controller  integrates  the  feed¬ 
back  error,  the  off-diagonal  elements  have  the  character¬ 
istics  of  a  disturbance  rejection  transfer  function  and  a 
final  value  of  zero.  Using  the  procedure  of  Reference  13, 
the  asymptotic  transfer  function  matrix  (F)  is  computed  for 
the  design  parameters  listed  previously. 


r(X)  = 


.75g/ (A  +  .75g) 


4.0/(A  +  4.0)  -A/(A+4.0)  (4-16 

0  .02g/ (A  +  .02g) 


It  is  evident  from  (4-16)  that  the  cross-coupling  term  is 
not  affected  by  gain  and  is  present  regardless  of  the 
choice  of  design  variables,  although  its  location  can  be 
changed  with  the  pitch  rate  feedback  gain.  The  signifi¬ 
cance  of  this  element  is  that  regardless  of  design  param¬ 
eters  chosen,  there  is  coupling  of  flight  path  angle  and 
angle  of  attack.  At  less  than  infinite  gain  this  element 
has  an  additional  term  that  is  a  function  of  gain  and  the 
third  diagonal  element  of  the  sigma  matrix  (o^)* 

12,3  (A)  =  ■  £^3  g)  ^  /  (^  +  4)  (A  +  EC3  g)  (4-17 


Judicious  selection  of  may  help  to  reduce  the  coupling 
but  might  also  degrade  the  desired  tracking  response,  of 
course,  the  actual  system  cannot  operate  with  infinite 
gain,  and  Figures  4.2  through  4.10  are  the  step  responses 


for  the  actual  transfer  functions  and  show  that  at  finite 


gains  the  system  exhibits  reasonably  tight  tracking  in  the 
diagonal  elements.  Also,  that  there  is  significant  coupling 
that  is  not  predicted  by  the  asymptotic  properties,  par¬ 
ticularly  in  the  output  responses  corresponding  to  the 
velocity  input  command. 

Actuator  Dynamics  and  Computational  Time  Delay. 

Once  actuators  and  time  delay  are  included  in  the  system 
simulation,  the  program  does  not  compute  all  of  the  closed- 
loop  system  roots.  To  design  a  controller  under  these  and 
all  subsequent  considerations,  the  linear  characteristics 
of  the  system  are  used  to  estimate  an  acceptable  control 
law,  and  then  a  complete  simulation  is  performed.  An 
iterative  technique  is  used  to  find  a  controller  that 
results  in  an  acceptable  time  response.  Rather  than  dis¬ 
cuss  the  results  of  each  step  in  the  iterative  process, 
this  thesis  presents  examples  of  the  effects  of  the 
principal  parameters  used  to  achieve  the  desired  results. 
Since  the  inclusion  of  the  various  "real  world"  effects 
does  modify  the  performance,  it  is  important  to  consider 
the  practical  aspects  of  aircraft  flight  control  systems. 

The  first  of  these  practical  aspects  to  limit  the  flexi¬ 
bility  of  controller  design  is  control  surface  saturation. 
The  previous  section  shows  the  possible  methods  of  influ¬ 
encing  the  transient  response  by  manipulating  the  system's 


asymptotic  properties.  However,  the  control  surface  author¬ 
ity  limits  the  choice  of  design  variables  dramatically. 

In  fact,  after  stabilizing  the  aircraft,  the  most  diffi¬ 
cult  problem  encountered  is  obtaining  a  reasonably  fast, 
critically  damped  response  without  surface  deflections  that 
exceed  the  physical  limits.  In  all  design  cases,  the  input 
to  the  system  is  that  of  a  six  degree  decrease  in  flight 
path  angle  from  level  flight,  and  the  desired  output  is  a 
smooth  and  non-oscillatory  flight  path  response  with  little 
or  no  change  in  angle  of  attack  and  velocity.  This  input, 
although  not  the  ma.rimum  possible  based  on  the  steady-state 
transfer  function,  is  chosen  as  a  reasonably  challenging 
input  command  for  a  jet  aircraft.  Table  4.2  lists  the  maxi¬ 
mum  steady-state  flight  path  change  and  the  limiting  surface 
for  each  flight  condition.  The  maximum  possible  command 
input  for  this  aircraft  is  not  limited  by  the  steady  state 
control  surface  requirements  but  by  the  transient  overshoot 
in  the  control  sur.‘"aces.  Although  not  shown  in  Table  4.2, 
the  transient  control  surface  deflections  for  the  maximum 
maneuver  at  any  flight  condition  must  not  exceed  the  deflec¬ 
tion  limits.  To  assist  in  minimizing  initial  control  sur¬ 
face  overshoot  the  commanded  input  is  slightly  ramped 
(0.8  seconds)  and  smoothed.  Two  additional  inputs  are  con¬ 
sidered  to  simulate  an  aborted  landing  situation  and  the 
final  roundout/ flare  for  touchdown.  The  designs  are  not 
altered  for  these  maneuvers,  rather  the  adaptability  of 


Flight 

Condition* 


Maximum  Flight 
Path  Change 


Limiting  Surface 


2 

3 

5 

6 


48.7  degrees 

25.7  degrees 
14.3  degrees 

15.5  degrees 
17.2  degrees 

11.5  degrees 


rotating  vanes 
stabilator 
stabilator 
canard 
canard 
canard 


*See  page  22  for  flight  conditions. 

landing  approach  designs  to  other  maneuvers  appropriate  to 
the  flight  condition  is  demonstrated. 

1.  Sigma  Weighting  Matrix.  It  is  difficult  to 
identify  one  parameter  to  use  to  achieve  a  particular 
effect,  but  several  trends  are  evident  from  analysis  of  the 
system  response.  The  following  paragraphs  describe  these 
trends  and  present  graphical  justification.  Analysis  of 
the  asymptotic  transfer  function  indicates  that  at  high 
gain  the  first  output  (velocity)  should  be  uncoupled  from 
the  other  two  and  the  principal  control  variable  is  the 
corresponding  sigma  matrix  element.  The  velocity  can  be 
controlled  relatively  independently  of  the  other  two  out¬ 
puts  as  is  demonstrated  in  Figures  4.11  through  4.14. 


As  the  sigma  element  corresponding  to  velocity  is  increased 
the  velocity  transient  diminishes  without  significantly 
changing  the  other  two  outputs.  As  the  asymptotic  trans¬ 
fer  function  matrix  suggests,  the  angle  of  attack  and 
flight  path  are  not  decoupled  and  changing  either  of  the 
two  corresponding  sigma  elements  influences  both  outputs 
(Figures  4.15,  4.16,  and  4.17).  Changing  the  element 
to  minimize  coupling  (Equation  (4-17))  proves  to  be  imprac¬ 
tical  due  to  the  control  surface  deflection  and  rate 
requirements.  Thus,  choosing  =  2  theoretically  elimi¬ 
nates  coupling  but  is  prohibitive  from  a  practical  stand¬ 
point  as  it  results  in  system  instability.  Obviously, 
roots  that  do  not  appear  in  an  asymptotic  analysis  do  influ 
ence  the  system  response. 

2.  Integrator  Gain.  As  discussed  earlier,  the 
ratio  of  integral  to  proportional  control  defines  three 
slow  roots  that  do  not  appear  in  any  of  the  asymptotic 
transfer  functions.  However,  at  finite  gains,  the  slow 
roots  are  observable  and  significantly  influence  the  out¬ 
puts.  Figure  4.18  shows  how  increasing  integral  gain 
causes  overshoot  in  the  commanded  output.  This  situation 
can  cause  considerable  variation  in  the  apparent  steady 
state  response  of  the  system  as  flight  conditions  change, 
since  the  step  response  may  not  settle  within  the  time 
period  of  interest.  This  phenomenon  is  not  shown  n.jre  but 
can  be  seen  in  the  following  chapter.  It  is  important  to 
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realize  that  pilots  expect  changes  in  the  response  of  their 
aircraft  with  changing  airspeed,  altitude,  and  gross  weight, 
and  it  is  a  fundamental  aspect  of  learning  to  fly  a  par¬ 
ticular  aircraft. 

3.  Actuator  Gain.  The  gain  of  each  of  the  surface 
actuators  may  also  be  adjusted  and,  in  particular,  the 
stabilator  actuator  gain  (K^)  proves  useful  in  reducing 
surface  and  output  oscillations.  Figures  4.19  through  4.24 
show  that  increasing  the  stabilator  actuator  gain  by  25 
percent  (K^  =  1.25)  damps  the  oscillation  and  further 
increase  results  in  an  overdamped  response. 

4.  Pitch  Rate  Feedback.  As  discussed  previously 
in  this  chapter,  the  measurement  matrix  results  in  a  minor 
loop  feedback  of  the  pitch  rate.  It  is  important  to  point 
out  that  this  is  not  pitch  rate  or  pitch  damping  control 
because  the  pitch  rate  is  not  compared  to  a  commanded  pitch 
rate  or  pitch  angle.  It  is  included  as  a  signal  to  augment 
the  angle  of  attack  feedback  signal.  Therefore,  the  pitch 
rate  feedback  gain  should  not  be  expected  to  affect  pri¬ 
marily  pitch  damping  as  predicted  by  conventional  methods 
(10:59).  However,  it  is  obvious  from  Figures  4.25  through 
4.30  that  the  pitch  rate  feedback  gain  most  definitely 
influences  the  system's  oscillation,  and  increasing  gain 
results  in  increased  damping. 
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Fig.  4.20.  Surface  Deflection,  K_  =  1 
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Sensor  Dymamics,  The  addition  of  sensor  dynamics 
tends  to  destabilize  the  system  response  and  is  generally 
compensated  for  by  reduction  of  the  gain  (sigma  matrix 
multiplier)  and  minor  modification  of  the  sigma  weighting 
elements.  The  same  qualitative  trends  observed  in  the 
design  of  the  basic  plant  with  actuators  and  time  delay  hold 
true  with  the  addition  of  sensors  but  usually  changing 
other  parameters  is  not  required.  A  complete  presentation 
of  the  results  of  applying  these  observations  are  contained 
in  Chapter  V. 

Control  Surface  Nonlinearity.  Simulating  the  con¬ 
trol  surface  nonlinearity  described  in  the  preceding 
chapter  also  produces  some  mild  instability  in  some  cases. 
The  situation  is  easily  remedied  by  a  reduction  in  either 
the  first  sigma  matrix  element  or  epsilon.  The  source  of 
this  instability  is  probably  the  point  discontinuity  in 
the  control  surface  derivatives  as  they  pass  through  zero 
angle  of  attack.  In  reality,  the  function  is  a  smooth 
parabolic  arc  and  should  not  cause  the  erratic  behavior 
of  Figure  4.32.  Even  with  the  discontinuity  of  the  sur¬ 
faces,  the  outputs  (Figure  4.31)  are  not  significantly 
affected.  Figures  4.33  and  4.34  show  the  surface  deflec¬ 
tions  and  outputs  after  the  sigma  elements  are  adjusted. 

Parameter  Variation.  To  demonstrate  the  ability 
of  the  design  method  to  provide  satisfactory  control  in  the 
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presence  of  plant  uncertainty,  each  controller,  tailored 
for  one  particular  flight  condition,  is  tested  at  each  of 
the  remaining  five  flight  conditions.  No  changes  are  made 
to  any  of  the  controllers  during  this  analysis  since  the 
objectives  are  strictly  demonstrative. 

Sensor  Noise .  To  demonstrate  the  influence  of 
sensor  noise,  the  noise  values  of  Appendix  C  are  added  to 
the  simulation.  In  addition,  each  noise  value  is  increased 
individually  to  its  threshold  level,  defined  as  the  level 
above  which  the  system  diverges  beyond  control  surface 
limits  within  the  twenty  second  simulation  period.  A 
Monte  Carlo  analysis,  consisting  of  five  independent  simula¬ 
tions,  is  conducted  for  each  of  the  noise  configurations 
to  obtain  a  mean  response.  MULTI  does  not  currently 
include  the  capability  to  calculate  the  variance  for  a 
Monte  Carlo  simulation.  As  a  result,  the  variance  must  be 
inferred  from  the  qualitative  analysis  of  each  of  the  five 
simulations.  All  noise  analysis  is  conducted  at  flight 
condition  1  and  no  changes  are  made  to  the  controller  design 
The  results  of  this  analysis  are  contained  in  the  following 
chapter . 

Additional  Maneuvers.  Two  additional  maneuvers 
are  performed  at  flight  condition  1  to  demonstrate  that  the 
controllers  are  not  tuned  to  just  this  one  maneuver,  and 
they  are  capable  of  performing  other  maneuvers  that  must  be 


considered  in  the  landing  phase  of  flight.  Assuming  that 
the  aircraft  can  establish  a  steady  state  descent  on  a 
desired  flight  path  in  preparation  for  landing  (the  result 
of  the  first  maneuver  studied) ,  there  are  basically  two 
possible  follow-on  maneuvers:  a  roundout  and  flare  to 
touchdown,  and  a  go-around  or  aborted  landing.  The  first 
maneuver  entails  a  smooth  return  to  a  level  flight  attitude 
just  as  the  landing  gear  touch  down.  It  is  necessary  to 
insure  that  there  is  a  significant  nose  high  pitch  attitude 
at  touchdown  to  prevent  nose  wheel  damage,  and  the  velocity 
should  be  minimized  to  reduce  the  landing  roll  and  touch¬ 
down  dispersion.  At  the  velocity  of  flight  condition  1 
(120  knots)  the  level  flight  pitch  attitude  is  already  more 
than  eleven  degrees  nose  high.  This  is  more  than  sufficient 
for  the  flare  and  no  increase  in  angle  of  attack  is  com¬ 
manded.  It  is  also  assumed  that  the  equilibrium  velocity 
will  be  maintained,  rather  than  commanding  it  to  a  level 
dangerously  close  to  a  stall  condition.  The  maneuver  is 
thereby  simplified  to  a  commanded  return  to  zero  flight 
path  angle  with  no  input  to  either  the  angle  of  attack  or 
flight  path.  The  second  maneuver  assumes  that  for  some 
reason  the  landing  is  refused  and  the  descent  must  be 
immediately  terminated.  It  is  further  assumed  that  after 
stabilizing  briefly  at  a  level  or  near  level  condition  the 
landing  is  resumed.  The  purpose  of  this  maneuver  is  two¬ 
fold;  to  demonstrate  that  the  aircraft  will  smoothly  and 


quickly  return  to  near  equilibrium  conditions,  and  that 
despite  the  relatively  long  theoretical  settling  time 
resulting  from  low  integral  gain,  the  aircraft  can  be 
rapidly  maneuvered  without  inducing  unacceptable  flight 
characteristics.  This  analysis  is  conducted  with  the  model 
of  flight  condition  1  with  actuators,  computational  time 
delay,  sensor  dynamics  and  surface  nonlinearities  included 
in  the  simulation.  No  changes  are  made  to  the  controller 
design  during  this  demonstration  and  the  results  are  pre¬ 
sented  in  Chapter  V. 

4. 4  Summary 

In  general,  it  is  very  difficult  to  describe  one's 
thought  processes  as  he  proceeds  through  a  design  of  a 
complex  system  by  a  method  that  requires  some  trial  and 
error.  There  are  a  number  of  clues  contained  in  the  theory 
used  to  develop  the  design  method,  but  often  these  relation 
ships  are  not  apparent  until  after  successful  results  are 
achieved.  This  chapter  is  an  attempt  to  present  sys¬ 
tematically  the  abstract  procedures  that  are  personal  to 
an  individual  and  peculiar  to  a  specific  mathematical  model 
in  a  quantitative  and  objective  manner  that  may  be  useful 
to  future  users  of  the  Porter  method  of  multivariable  out¬ 
put  feedback  control. 


V.  Results 


5 . 1  Introduction 

This  chapter  presents  the  results  of  applying  the 
design  procedure  of  Chapter  IV  to  the  STOL  F-15  aircraft 
model.  The  chapter  begins  with  a  discussion  of  the  general 
format  of  the  data  presented  and  then  proceeds  to  the  spe¬ 
cific  results,  starting  with  the  plant  plus  actuators  and 
computational  delay,  and  progressing  to  a  simulation  of  the 
plant  with  actuators,  computational  time  delay,  sensor 
dynamics,  nonlinear  control  surface  effects  and  sensor 
noise.  Next,  simulations  of  the  response  of  each  con¬ 
troller  to  plant  parameter  variation  are  presented  to 
establish  the  robustness  of  the  controllers  designed  at 
each  flight  condition.  Finally,  two  additional  maneuvers 
are  simulated,  demonstrating  the  controlled  aircraft's 
ability  to  perform  other  essential  maneuvers  characteristic 
of  a  landing  situation. 

5.2  Format  of  the  Results 

Commanded  Input.  Except  for  the  results  of  Section 
5.8  of  this  chapter,  all  simulations  represent  the  response 
of  the  complete  closed-loop  system  to  a  smoothed,  ramped, 
step  command  input  to  the  flight  path  angle  channel  of  the 
controller  (Figure  5.1).  The  commanded  flight  path  angle 


input  is  ramped  over  0.8  seconds  and  has  a  final  constant 
value  of  -6.0  degrees.  This  input,  along  with  a  commanded 
input  of  zero  in  the  other  two  channels,  represents  a  con¬ 
stant  velocity,  constant  angle  of  attack  descent  command  at 
an  angle  of  six  degrees  to  the  horizontal.  Section  5.8 
presents  two  additional  maneuvers  described  in  detail  in 
that  section. 

Plotted  Data .  The  graphical  data  presented  con¬ 
sists  of  two  plots  for  each  simulation.  The  first  plot  is 
a  time  history  of  the  three  output  variables  (u,  a,  and  y) 
expressed  in  units  of  degrees.  The  second  plot  is  a  time 
history  of  the  deflection  of  each  of  the  control  surfaces 
(6^,  6gf  and  6^)  also  expressed  in  degrees.  All  variables 
plotted  are  perturbation  values  and  the  magnitude  plotted 
represents  the  deviation  from  an  equilibrium  value. 

Tabular  Data.  In  addition  to  graphical  data,  the 
result  of  the  design  procedure  at  each  flight  condition 
and  the  corresponding  time  responses  are  presented  in 
tabular  form.  The  design  data  tables  (Tables  5.1,  5.3, 
and  5.5)  contain  the  pertinent  design  parameters  and  result¬ 
ing  controller  matrix  at  each  flight  condition.  There  are 
two  intentional  omissions  from  these  tables.  First,  since 
the  integral  control  matrix  (Kl)  is  simply  equal  to  the 
product  of  the  proportional  control  matrix  (Kl)  and  a,  Kl 
is  not  presented.  Also,  the  sigma  weighting  matrix 


multiplier  (e)  is  merely  a  design  tool  and  the  true  sigma 
weighting  matrix  is  the  product  of  "e"  and  the  diagonal 
sigma  matrix  (I).  The  resulting  diagonal  sigma  matrix 
elements  are  included  in  the  tables.  Presenting  the  sigma 
elements  in  this  manner  allows  direct  comparison  between 
flight  conditions  without  regard  to  "e."  The  simulation 
results  tables  (Tables  5.2,  5.4,  and  5.6)  contain  the 
figures  of  merit  calculated  by  MULTI  for  each  of  the  out¬ 
puts  at  the  various  flight  conditions.  The  tables  include 
figures  for  the  final  value  of  the  flight  path  angle. 

With  integral  control  the  system  is  guaranteed  to  have  zero 
steady  state  error  and  therefore  the  final  value  is  the 
same  as  the  commanded  input.  However,  since  the  gain  of 
the  integral  control  is  so  low  (one  hundredth  of  the  propor 
tional  gain) ,  the  simulation  may  not  reach  steady  state 
within  the  simulation  time.  The  final  value  is  therefore 
taken  to  be  the  value  at  the  end  of  the  simulation.  The 
settling  time  is  also  computed  with  respect  to  this  final 
value.  These  values  are  presented  to  indicate  the  quasi¬ 
steady  state  response  of  the  system  within  the  time  period 
of  interest.  Provided  the  system's  transient  response  is 
essentially  completed  in  the  simulation  time  these  values 
provide  reasonable  quantification  of  the  settling  time  of 
the  response. 


5.3 


Plant  Plus  Actuators  and 
Computational  Delay 

The  design  parameters  obtained  for  the  six  flight 
conditions  are  shown  in  Table  5-1.  In  general,  the  sigma 
element  corresponding  to  velocity  is  at  least  one  order  of 
magnitude  higher  than  the  other  two  elements.  The  ratio  of 
integral  to  proportional  control  is  0.01  for  all  flight 
conditions.  The  simulation  results  using  these  controllers 
are  contained  in  Table  5.2,  and  Figures  5.2  through  5.13. 

The  rise  times  and  settling  times  of  Table  5.2 
clearly  indicate  that  the  slowest  responses  are  for  the 
flight  conditions  in  which  the  aircraft  weight  is  the  high¬ 
est.  This  is  an  expected  result  since  the  higher  mass  and 
pitch  moment  of  inertia  certainly  affects  the  open  loop 
bandwidth  of  the  plant.  These  results  are  substantiated 
by  qualitative  inspection  of  the  corresponding  figures 
as  well.  The  speed  of  the  response  is  also  apparently 
affected  by  dynamic  pressure  as  the  response  at  a  particu¬ 
lar  weight  is  the  slowest  at  the  lowest  dynamic  pressure 
(flight  conditions  3  and  6).  This  is  no  surprise  either, 
since  the  controlling  forces  generated  by  the  aerodynamic 
surfaces  are  directly  proportional  to  the  dynamic  pressure. 
It  is  important  to  keep  in  mind  that  these  comparisons  are 
made  between  different  controllers  at  different  flight 
conditions.  Each  controller  is  chosen  subjectively  as  the 
best  one  found  for  that  given  flight  condition. 


TABLE  5.1 


Flight 

Condition 


Sigma 


-0  Matrix 


DESIGN  PARAMETERS  FOR  THE  PLANT  PLUS  ACTUATORS 
AND  COMPUTATIONAL  TIME  DELAY 


1.25 

0.025 

0.025 

-7.401E-4 

5.033E-4 

-5.740E-2 

5.146E-2 

-3.500E-2 

2.569E-3 

3. 160E-1 
1.115E-1 
-3.273E-2 

1.25 

0.025 

0.025 

6.656E-4 

-1.262E-1 

-5.210E-2 

6.049E-2 
-4 .040E-2 
-1.067E-3 

2. 923E-1 
8. 379E-1 
-2.701E-2 

1.11 

0.027 

0.018 

-1.647E-3 

-1.777E-2 

-6.869E-2 

5. 737E-2 
-4.333E-2 
-2.599E-3 

3.044E-1 

1.138E-1 

-2.082E-2 

1.20 

0.056 

0.034 

-1.472E-4 

-1.075E-2 

-6.635E-2 

5.141E-2 

-4.743E-2 

-3.995E-3 

3.312E-1 

1.543E-1 

-2.887E-2 

1.313 

0.0338 

0.0188 

-2.194E-2 

-1.487E-3 

-9.024E-2 

4.866E-2 

-2.677E-2 

-3.434E-3 

2.768E-1 

1.579E-1 

-4.028E-2 

1.000 

0.015 

0.015 

-1.042E-2 

-5.522E-3 

-6.397E-2 

5.418E-2 
-2.020E-2 
-4 .620E-3 

2.975E-1 
1. 078E-1 
-3.201E-2 
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TABLE  5.2 

FIGURES  OF  MERIT  FOR  PLANT  PLUS  ACTUATORS 
AND  COMPUTATIONAL  TIME  DELAY 


Flight 

Condition 

Output 

Peak 

Value 

Time 

To  Peak 

Final 

Value* 

Settling 

Time 

1 

u 

0.158 

2.60 

★  * 

*  * 

a 

-0.558 

3.20 

** 

*  * 

Y 

-6.305 

6.00 

-6.196 

4.60 

2 

u 

0.196 

2.00 

*  * 

*  * 

a 

-1.287 

2.80 

** 

*  * 

't' 

-6.746 

5.20 

-6.491 

6.60 

3 

u 

0.177 

4.60 

*  ★ 

*  it 

a 

-0.426 

3.80 

★  * 

*  * 

y 

-7.631 

11.2 

-7.437 

14.4 

4 

u 

0.210 

7.60 

*  * 

*  ★ 

a 

1.376 

1.20 

★  ★ 

★  ★ 

Y 

-7.162 

12.2 

-7.055 

7.40 

5 

u 

0.123 

4.00 

★  * 

** 

a 

0.479 

1.20 

*  * 

*  * 

Y 

-6.209 

17.4 

-6.209 

7.60 

6 

u 

0.194 

3.40 

** 

*  ★ 

a 

-8.876 

4.20 

** 

★  * 

y 

-6.812 

8.80 

-6.711 

6.60 

Fig.  5.6.  Outputs  with  Actuators  and  Computational 
Delay,  Flight  Condition  3 
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Fig.  5.8.  Outputs  with  Actuators  and  Computational 
Delay,  Flight  Condition  4 


Fig.  5.9.  Surface  Deflections  with  Actuators  and 
Computational  Delay,  Flight  Condition  4 
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The  plots  at  each  flight  condition  demonstrate 
relatively  good  decoupling  of  the  outputs  at  all  flight 
conditions,  with  the  greatest  coupling  occurring  at  low 
weight  and  low  dynamic  pressure  (Figures  5.4  and  5.12). 

In  general,  the  responses  are  smooth,  non-oscillatory ,  and 
reasonably  fast.  In  most  of  the  surface  deflection  plots 
it  is  evident  that  the  canard  deflects  almost  to  its  deflec 
tion  limit  (-35  degrees) .  This  is  by  design  to  achieve 
the  maximum  performance  of  the  system.  However,  in  the 
performance  of  the  maneuver  the  canard  also  is  at  the  limit 
of  and  sometimes  exceeds  its  maximum  deflection  rate  (23 
degrees  per  second) .  This  is  a  characteristic  of  high  gain 
controllers.  Future  design  efforts  with  this  aircraft 
should  consider  implementing  rate  limits  in  the  MULTI  simu¬ 
lation  to  account  for  this  realistic  limitation. 

5 . 4  Plant  Plus  Actuators  Delay 

and  Sensors 

The  addition  of  sensor  dynamics  tends  to  destabilize 
the  system  response  using  the  controllers  of  Section  5.3. 
Therefore,  the  controllers  are  tailored  to  optimize  the 
response  including  the  sensors  and  the  results  of  this 
adjustment  are  tabulated  in  Table  5.3.  The  new  controllers 
generally  feature  a  reduction  in  the  sigma  element  corres¬ 
ponding  to  velocity  and  minor  adjustments  to  the  other  two 
channels.  The  simulation  results  are  contained  in  Table  5.4 
and  Figures  5.14  through  5.25. 


TABLE  5.3 


DESIGN  PARAMETERS  FOR  THE  PLANT  PLUS  ACTUATORS 
COMPUTATIONAL  TIME  DELAY  AND  SENSOR  DYNAMICS 


Flight 
Condition 


Sigma 


Matrix 


0.75 

0.02 

0.02 

-4.441E-4 
3.020E-4 
-3. 444E-2 

4. 117E-2 
-2.800E-2 
2.055E-3 

2.528E-1 
8. 921E-2 
-2.618E-1 

0.750 

0.030 

0.025 

3.994E-4 

-7.575E-3 

-3.126E-2 

7.259E-2 
-4 .848E-2 
-1.280E-3 

2.923E-1 

8.379E-2 

-2.701E-3 

0.809 

0.023 

0.015 

-1.138E-3 

-1.228E-2 

-4.746E-2 

4.909E-2 

-3.707E-2 

-2.224E-3 

2.511E-1 

9.386E-2 

-1.718E-2 

0.900 

0.056 

0.033 

-1.104E-4 

-8.066E-3 

-4.976E-2 

5.141E-2 

-4.743E-2 

-3.995E-3 

3.214E-1 

1.497E-1 

-2.802E-2 

0.825 

0.023 

0.018 

-1.379E-2 

-9.344E-4 

-5.672E-2 

3.224E-2 

-1.784E-2 

-2.289E-3 

2. 584E-1 
1.474E-1 
-3.759E-1 

1.200 

0,056 

0.034 

1.476E-4 

1.078E-2 

6.652E-2 

5.143E-2 

-4.614E-2 

4.006E-3 

3.313E-1 

1.637E-1 

2.894E-2 

u 


The  same  observations  that  are  discussed  in  the 
preceding  section  regarding  aircraft  weight  and  dynamic 
pressure  apply  when  sensor  dynamics  are  included.  Com¬ 
paring  the  peak  values  of  the  angle  of  attack  and  velocity 
from  Tables  5.2  and  5.4,  it  is  evident  that  the  addition  of 
sensor  dynamics  tends  to  increase  coupling  between  the 
input/output  channels.  The  increase  is  relatively  small 
and  appears  to  be  within  acceptable  limits.  Overall,  it 
is  possible  to  obtain  responses  nearly  as  fast  as  without 
sensors  at  the  expense  of  a  slight  increase  of  coupling  and 
increased  overshoot  of  the  commanded  input.  The  responses 
are  still  smooth  and  non-oscillatory . 

5.5  Plant  Actuators,  Delay,  Sensors, 

and  Surface  Nonlinearity 

The  nonlinearity  of  the  drag  on  the  control  surface 
as  modeled  in  this  thesis  induces  a  slightly  erratic 
behavior  of  the  control  surfaces  as  the  surface  passes 
through  zero  angle  of  attack.  This  is  most  likely  caused 
by  the  discontinuity  in  the  modeled  drag  derivative  that 
does  not  exist  in  nature.  However,  this  problem  is  largely 
resolved  by  a  reduction  in  one  or  more  of  the  sigma  ele¬ 
ments  as  necessary.  The  resulting  controller  designs  are 
shown  in  Table  5.5.  The  simulation  results  with  the  modeled 
surface  nonlinearities  included  are  contained  in  Table  5.6 


and  Figures  5.26  through  5.37. 


TABLE  5.5 


DESIGN  PARAMETERS  FOR  THE  PLANT  PLUS  ACTUATORS, 
COMPUTATIONAL  TIME  DELAY,  SENSOR  DYNAMICS 
AND  SURFACE  NONLINEARITIES 


Flight 


Condition 

a 

Sigma 

Kq  Matrix 

1 

.01 

0.450 

0.025 

0.025 

-2.665E-4 

1.812E-4 

-2.066E-2 

5.095E-2 

-3.465E-2 

2.543E-3 

3.128E-1 

1.104E-1 

-3.240E-2 

2 

.01 

0.750 

0.030 

0.025 

3.994E-4 

-7.575E-3 

-3.126E-2 

7.259E-2 

-4.848E-2 

-1.280E-3 

2.923E-1 
8. 379E-2 
-2.701E-3 

3 

.01 

0.578 

0.023 

0.015 

-8.129E-4 

-8.877E-3 

-3.390E-2 

4. 909E-2 
-3.707E-2 
-2.224E-3 

2.511E-1 

9.386E-2 

-1.718E-2 

4 

.01 

0.750 

0.056 

0.033 

~9.200E-6 

-•6.721E-4 

-4.147E-3 

5.141E-2 

-4.743E-2 

-3.995E~3 

3.214E-1 
1.497E-1 
-2 .802E-2 

5 

.01 

0.600 

0.025 

0.017 

-1.003E-2 

-6.796E-4 

-4.125E-2 

3.244E-2 
-1. 784E-2 
-2.289E-3 

2.584E-1 

1.474E-1 

-3.759E-2 

6 

.01 

1.200 

0.056 

0.034 

1.476E-4 

1.078E-2 

6.652E-2 

5.143E-4 

-4.614E-2 

4.006E-3 

3.313E-3 
1.637E-1 
2. 894E-2 

K  =  .25 


TABLE  5.6 


«• 


FIGURES  OF  MERIT  FOR  PLANT  PLUS  ACTUATORS, 
COMPUTATIONAL  TIME  DELAY  SENSORS  AND 
SURFACE  NONLINEARITY 


Flight 

Condition 

Output 

Peak 

Value 

Time 

To  Peak 

Final 

Value* 

Settling 

Time 

1 

u 

0.467 

0.80 

it  it 

** 

ct 

-0.651 

3.20 

it* 

** 

Y 

-6.388 

5.40 

-6.197 

6.60 

2 

u 

0.338 

1.80 

** 

** 

a 

-1.190 

2.80 

** 

*  * 

Y 

-6.619 

4.80 

-6.494 

3.80 

3 

u 

0.432 

9.80 

** 

** 

a 

-0.484 

4.20 

** 

** 

Y 

-8.521 

13.6 

-8.323 

15.8 

4 

u 

2.732 

9.80 

** 

** 

a 

1.589 

1.40 

** 

Y 

-7.103 

12.6 

-7.015 

7.60 

5 

u 

0.284 

6.40 

*  * 

a 

0.729 

1.40 

** 

** 

Y 

-6.366 

10.2 

-6.152 

13.4 

6 

u 

0.387 

3.40 

** 

** 

a 

-0.952 

4.20 

★  * 

** 

Y 

-6.820 

8.20 

-6.686 

8.60 

% 


100 


v>y  •rli 


if 


The  output  responses  are  basically  unaffected  by  the 
surface  nonlinearities  other  than  the  expected  result  of 
altered  sigma  matrix  elements.  Increased  coupling,  par¬ 
ticularly  with  respect  to  velocity,  is  evident  from  the 
reduction  in  the  corresponding  sigma  elements  for  flight 
conditions  1,  3,  4,  and  5.  The  flight  path  angle  response 
is  not  significantly  degraded  from  the  adjustments  for  the 
surface  nonlinearity.  Control  surface  deflections,  however, 
still  exhibit  a  slight  jump  (particularly  the  rotating 
vanes)  that  cannot  be  avoided  without  further  reduction  of 
the  velocity  sigma  element  or  possibly  a  smoothing  of  the 
surface  drag  derivative  at  angles  of  attack  near  zero. 

This  erratic  behavior  does  not  appear  in  the  outputs  with 
sufficient  magnitude  to  warrant  redesign  to  attempt  to 
remove  it . 

5.6  Addition  of  Sensor  Noise 

The  controller  design  of  Section  5.5  for  flight 
condition  1  is  next  subjected  to  simulations  in  which  inde¬ 
pendent  white,  gaussian  noise  is  injected  into  each  of  the 
four  measured  quantities  (u,  a,  y,  and  q) .  The  first  simu¬ 
lation  includes  realistic  noise  levels  for  a  comparable  air¬ 
craft  (Grumman  F-14  Tomcat,  Reference  12) .  This  data  is 
presented  in  Table  5.7  as  derived  in  Appendix  C  (units  are 
radians) .  Note  that  no  realistic  values  of  velocity  measure 
ment  noise  are  ^-ovided  and  it  is  therefore  omitted  from  the 
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TABLE  5.7 


SENSOR  NOISE  DATA 


Measured  Quantity 

Noise  Mean 

Variance 

a 

0.0 

1.220E-5 

Y 

0.0 

1.309E-5 

q 

0.0 

3.220E-5 

simulation.  To  establish  statistical  validity,  the  simula¬ 
tion  is  performed  five  times  with  a  different  random  noise 
seed.  The  results  of  the  first  four  simulations  and  a  mean 
response  are  shown  in  Figures  5.38  through  5.47.  Com¬ 
paring  these  results  to  Figures  5.26  and  5.27,  the  response 
is  indistinguishable  from  the  simulation  without  noise. 
Although  this  is  a  satisfying  result,  it  does  not  yield  any 
information  on  the  system's  response  to  higher  noise  levels. 
Therefore,  the  noise  levels  are  increased  individually  to  a 
value  that  appears  to  be  the  threshold  that  the  system  can 
withstand.  The  threshold  is  defined  in  this  case  as  the 
highest  noise  level  that  does  not  result  in  divergence 
beyond  control  surface  limits  within  the  twenty-second  simu¬ 
lation  time.  Table  5.8  contains  the  threshold  values  of 
noise,  expressed  in  units  of  radians  (velocity  noise  units 
are  feet  per  second) .  Comparing  the  threshold  variance  on 
the  basis  of  signal  to  noise  ratio,  the  system  is  most 
sensitive  to  noise  in  the  velocity  measurement  and  least 


Fig.  5.42.  Outputs,  Realistic  Noise,  Simulation  3 


6* 

Fig.  5.43.  Surface  Deflections,  Realistic  Noise,  Simulation  3 


Fig.  5.45.  Surface  Deflections,  Realistic  Noise,  Simulation  4 
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sensitive  to  measurement  of  the  pitch  rate.  Since  the 
velocity  feedback  channel  has  the  highest  gain  (diagonal 
sigma  element)  it  is  to  be  expected  that  the  system  would 
be  most  sensitive  to  noise  in  that  channel.  The  signal  to 
noise  ratios  presented  here  are  calculated  based  on  the 
maximum  absolute  value  of  the  signal  of  interest  encoun¬ 
tered  during  the  simulation.  As  with  the  realistic  noise, 
five  independent  simulations  are  accomplished  at  each 
threshold  value  to  obtain  a  statistically  valid  sample  of 
the  responses.  Figures  5.48  through  5.87  contain  the  first 
four  responses  and  a  mean  response  for  each  of  the  four 
types  of  noise.  It  is  apparent  in  many  of  the  individual 
and  mean  simulations  that  the  noise  causes  the  control 
surfaces  to  gradually  diverge  in  time.  This  result  is  pre¬ 
dictable  since  the  white  gaussian  noise  input  is  being 
integrated  by  the  controller,  resulting  in  Brownia  motion 
or  random  walk  (10:154).  This  divergence  tends  to  be  hidden 


Fig.  5,54.  Outputs,  Velocity  Noise,  Simulation  4 


.  5.55,  Surface  Deflections,  Velocity  Noise,  Simulation  4 
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Fig,  5.68.  Outputs,  Flight  Path  Noise,  Simulation  1 
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Fig.  5.74.  Outputs,  Flight  Path  Noise,  Simulation  4 


Fig.  5.75.  Surface  Deflections,  Flight  Path 
Noise,  Simulation  4 
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Outputs,  Pitch  Rate  Noise,  Simulation  2 


Fig.  5.82.  Outputs,  Pitch  Rate  Noise,  Simulation  3 


Surface  Deflections,  Pitch  Rate  Noise, 
Simulation  4 


in  the  mean  response  because  of  the  effect  of  averaging. 

It  would  be  desirable  to  compute  the  variance  of  the 
responses  as  a  function  of  time  in  future  noise  investiga¬ 
tions  to  aid  in  the  analysis  of  noise  effects.  In  any  case, 
the  plots  of  the  effect  of  high  levels  of  noise  are  presented 
to  demonstrate  their  effects  on  the  linear  model  and  it 
should  be  noted  that  in  many  cases  the  surface  deflection 
rate  limits  are  exceeded,  thereby  invalidating  the  linear 
model.  The  significance  of  this  result  is  that  in  general 
it  is  not  advisable  to  drive  an  integral  controller  with 
unfiltered  noise.  A  practical  solution  is  to  place  a 
cascade  low  pass  filter  between  the  noise  source  and  the 
controller;  however,  this  thesis  does  not  pursue  solutions 
to  the  problems  associated  with  noisy  measurements. 

5 . 7  Parameter  Variation 

To  establish  the  robustness  of  the  designed  con¬ 
trollers  with  respect  to  plant  parameter  variation,  the 
control  matrices,  KO  and  Kl,  from  each  flight  condition 
are  used  to  control  the  plants  of  the  other  five  flight  con 
ditions.  Naturally,  not  all  of  these  simulations  are  suc¬ 
cessful  and  in  fact  some  are  unstable.  Table  5.9  indicates 
which  controllers  are  satisfactory  at  the  various  flight 
conditions.  The  simulations  that  are  stable  are  plotted 
in  Figures  5.88  through  5.111.  The  unstable  responses  are 
not  presented  because  they  quickly  exceed  the  boundaries  of 
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Fig.  5,102.  Outputs,  Controller  3,  Flight  Condition  4 


Fig.  5.103.  Surface  Deflections,  Controller  3 
Flight  Condition  4 
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TABLE  5.9 


PARAMETER  VARIATION  RESULTS 


Flii^t 

Condition 

Controller 

1 

2 

3 

4 

5 

6 

1 

Good 

Stable 

Stable 

Unstable 

Unstable 

Stable 

2 

Stable 

Good 

Unstable 

Unstable 

Unstable 

Stable 

3 

Stable 

Stable 

Good 

Unstable 

Unstable 

Stable 

4 

Stable 

Stable 

Stable 

Good 

Unstable 

Stable 

5 

Unstable 

l&TS  table 

Unstable 

Unstable 

Good 

Unstable 

6 

Stable 

Unstable 

Unstable 

Unstable 

Unstable 

Good 

a  valid  linear  model  and  are  meaningless  other  than  to 
establish  that  the  response  is  divergent.  All  of  the  stable 
responses  contain  a  lingering  transient  error  in  the  angle 
of  attack,  the  principal  reason  only  the  diagonal  entries 
in  Table  5.8  are  considered  "Good."  In  all  of  the  stable 
simulations  the  surface  deflection  limits  are  within  limits 


and  only  two  responses  (Figures  5.94  and  5.98)  exhibit 
objectionable  oscillatory  characteristics.  In  a  practical 
system,  it  is  often  convenient  to  use  one  controller  at  a 
variety  of  flight  conditions  to  avoid  the  complications  of 
gain  scheduling.  These  results  indicate  that  either  of 
two  controllers  (those  designed  for  flight  conditions  1 
and  6)  provide  acceptable  control  over  all  but  one  of  the 
remaining  flight  conditions  (flight  condition  5) .  The 


only  satisfactory  controller  for  flight  condition  5  is  the 
one  designed  at  that  condition. 

5 . 8  Special  Maneuvers 

For  the  reasons  outlined  in  Chapter  IV,  two  addi¬ 
tional  maneuvers  are  simulated  to  demonstrate  the  system's 
ability  to  perform  other  maneuvers  that  are  representative 
of  the  landing  flight  condition.  Both  of  these  maneuvers 
are  demonstrated  at  flight  condition  1  with  actuators, 
computational  time  delay,  sensor  dynamics,  and  surface  non- 
linearities  included.  The  first  maneuver  simulated  is  that 
of  an  approach  and  flare  as  the  aircraft  nears  the  ground. 
The  commanded  input  (Figure  5.112)  begins  with  a  command 
identical  to  the  input  of  the  previous  section.  After  the 
aircraft  has  sufficient  time  to  establish  a  steady-state 
descent  (ten  seconds)  the  flight  path  is  ramped  back,  to 
zero  over  four  seconds  to  simulate  the  roundout  and  flare. 

A  profile  of  the  resultant  flight  path  response  (altitude 
vs  time)  is  shown  in  Figure  5.113.  Figures  5.114  and  5.115 
contain  the  outputs  and  surface  deflections  during  the 
maneuver  and  demonstrate  smooth  performance  of  the  intended 
maneuver. 

The  second  maneuver  performed  simulates  an  aborted 
landing  in  which  the  initial  descent  is  commanded,  followed 
by  a  return  to  level  flight,  and  finally  the  initial  descent 
is  resumed  (Figure  5.116).  Figures  5.117  and  5.118  contain 


Sir^ULflTED  RP  =  SOflCh  RND  FlPRE  COMMAND  INPUT 

Fig.  5.112.  Flare  Command  Input 


Flight  Path  Command  Input 


the  system  response  to  this  maneuver  and  demonstrate  the 
controller's  ability  to  satisfactorily  respond  to  sequential 
inputs  with  only  a  minimum  time  to  respond  between  inputs. 


5.9  Summary 

The  results  contained  in  this  chapter  demonstrate 
the  application  of  the  design  procedure  of  Chapter  IV  to 
the  aircraft  model  of  Chapter  III.  The  unstable  open  loop 
plants  of  six  different  flight  conditions  are  stabilized 
and  satisfactory  results  are  obtained  despite  the  simulation 
of  the  destabilizing  effects  of  actuators,  computational 
time  delay,  sensor  dynamics  and  control  surface  nonlinear¬ 
ity.  The  effect  of  measurement  noise  is  simulated  and  the 
responses  are  found  to  be  unaffected  by  realistic  values  of 
noise.  Also,  threshold  values  of  noise  are  determined  and 
indicate  that  the  controlled  system  at  flight  condition  1 
is  most  sensitive  to  noise  in  the  velocity  measurement. 

The  controllers  at  two  flight  conditions  are  found  to  be 
sufficiently  insensitive  to  plant  parameter  variation  to 
provide  satisfactory  control  at  all  but  one  of  the  remaining 
flight  conditions.  Finally,  two  additional  maneuvers 
representative  of  the  landing  scenario  are  satisfactorily 
performed. 


readily  chosen.  Early  in  the  design  process  it  may  be 
difficult  to  obtain  user  definable  parameters  that  yield 
stable  results,  particularly  in  the  case  of  unstable  open 
loop  plants.  However,  after  some  experimentation  it  soon 
becomes  apparent  what  measures  must  be  taken  to  stabilize 
the  plant.  For  this  aircraft,  reducing  all  of  the  elements 
of  the  diagonal  weighting  matrix  and  the  ratio  of  integral 
to  proportional  gain  ( a  )  is  necessary  to  achieve  closed 
loop  stability.  However,  different  plants  may  well  have 
different  requirements  for  obtaining  stability.  The  design 
method  lacks  an  algorithmic  approach  to  initial  stabiliza¬ 
tion,  but  since  controller  design  is  a  simple  calculation, 
an  iterative  approach  using  simulation  results  is  accept¬ 
able.  Therefore,  recommendations  for  improvement  are  con¬ 
cerned  exclusively  with  the  simulation  program  MULTI. 

6 . 3  Improvements  to  MULTI 

There  are  several  recommendations  for  improving  the 
design  process  through  the  expansion  of  the  capabilities 
of  the  computer  aided  design  and  simulation  program  MULTI . 
The  previous  section  points  out  that  the  design  procedure 
relies  on  the  repetitive  use  of  MULTI  design  and  simulation 
features  to  achieve  the  desired  control  design.  The 
options  in  MULTI  that  perform  design  calculations  are  rela¬ 
tively  fast  and  use  computer  time  sparingly.  These  calcu¬ 
lations  include  the  computation  of  continuous  closed  loop 


roots,  provided  actuator  and  sensor  dynamics  are  not 
included  in  the  simulation.  Inspection  of  these  roots  pro¬ 
vides  an  initial  evaluation  of  stability  prior  to  a  complete 
simulation,  an  option  that  requires  considerable  central 
processor  time.  Since  actuator  and  sensor  dynamics  are 
assumed  linear,  their  effects  should  be  included  in  the 
calculation  of  closed  loop  transfer  functions,  allowing  the 
rapid  identification  of  unstable  designs. 

Currently,  the  control  surface  deflection  data  is 
inaccessible  for  plotting  when  actuator  dynamics  are 
included.  Rather,  the  data  plotted  data  is  the  input  to 
the  surface  actuators.  The  algorithm  that  creates  the  plot 
file  data  should  be  modified  to  obtain  the  actual  surface 
deflections.  Also,  control  surface  rate  limits  should  be 
implemented  in  the  program.  However,  position  and  rate 
saturations  often  result  in  instability  with  integral  con¬ 
trollers  and  it  is  recommended  that  the  implementation  of  a 
compensating  algorithm  be  considered. 

MULTI  features  algorithms  that  smooth  the  commanded 
inputs  prior  to  comparing  them  with  the  feedback  signal. 
Often  this  type  of  smoothing  is  accomplished  by  passing  the 
input  signal  through  a  linear  low  pass  filter.  Implementa¬ 
tion  of  this  type  of  filter  in  addition  to  the  smoothing 
may  help  to  reduce  the  high  initial  surface  deflections 
and  deflection  rates. 


It  is  further  recommended  that  the  simulation  of 
noise  be  expanded  to  include  a  Monte  Carlo  evaluation  of 
the  variance  as  a  function  of  time,  as  well  as  an  option 
to  include  filtering  of  the  measurement  noise. 

Finally,  the  recently  developed  methods  of  Professor 
Porter  et  al.  (17)  for  designing  proportional  plus  integral 
plus  derivative  control  should  be  added  as  an  option  to 
MULTI . 

6 . 4  Proposed  Future  Work 

This  thesis  presents  an  initial  evaluation  of  the 
effects  of  sensor  noise  on  the  longitudinal  control  of  the 
STOL/F-15  in  landing  configuration.  To  completely  evaluate 
the  noise  effects  and  practical  implementations  of  the  con¬ 
trollers,  it  is  necessary  to  implement  filter  algorithms 
and  quantitatively  determine  the  variance  of  the  output 
signals.  Furthermore,  the  ability  to  simulate  noise  in 
the  form  of  disturbance  inputs  exists  in  MULTI,  but  dis¬ 
turbance  effects  on  the  STOL/F-15  with  output  control  are 
currently  untested.  Further  study  of  this  aircraft  and 
design  method  should  include  both  of  these  evaluations. 

Longitudinal  control  of  the  F-15/ST0L  using  the 
Porter  techniques  of  multivariable  output  control  is  inves¬ 
tigated  in  this  thesis  as  well  as  in  a  parallel  effort 
using  the  combat  configuration  of  the  F-15/ST0L  (16) .  It 


is  recommended  that  the  Porter  methods  be  applied  to  the 
lateral  control  of  the  STOL/F-15. 

Finally,  recent  improvements  to  the  Porter  method 
allow  control  of  systems  with  rank  deficient  first  Markov 
parameters  without  the  use  of  the  measurement  matrix  (17) . 
As  a  result,  there  is  no  fundamental  finite  limit  to  either 
the  speed  or  decoupling  of  the  responses.  It  is  recom¬ 
mended  that  these  methods  of  proportional  plus  integral 
plus  derivative  control  be  applied  to  the  STOL/F-15. 
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Appendix  A;  Additions  to  MULTI 


Introduction 

During  the  course  of  this  thesis,  a  number  of 
changes  and  additions  were  made  to  the  computer  program 
MULTI  to  facilitate  current  and  future  research  efforts. 
This  appendix  describes  each  of  these  changes  and  additions 
as  well  as  providing  an  outline  of  the  entire  program  for 
the  benefit  of  future  programmers. 

List  of  Changes  and  Additions 

A.  Gaussian  noise  option 

B.  Custom  input  option 

C.  Suppression  of  actuators  and  sensors 

D.  Saving  memory  files  without  exit 

E.  Convert  input  vector  "u"  from  radians  to 
degrees 

F.  Plot  combination  of  states  and  inputs 

G.  Simulation  of  nonlinearities  peculiar  to 
aircraft 

H.  Calculate  initial  integrator  state  Z(0)  vector 

I.  Program  outline 

A.  Gaussian  Noise  Option 

1.  Description .  This  addition  gives  the  user  the 


option  of  simulating  zero  and  non- zero  mean,  white. 


gaussian,  noise  inputs  to  the  system  during  execution  of 
option  26 .  There  are  three  types  of  noise  inputs  avail¬ 
able,  distinguished  by  the  place  in  which  the  noise  is 
injected  into  the  linearized  model.  The  first  type,  out¬ 
put  measurement  noise,  is  that  noise  which  is  introduced 
by  the  sensors  used  to  measure  the  output  variables  being 
fed  back  to  the  controller.  The  second,  measurement  matrix 
noise,  is  identical  to  output  measurement  noise,  except 
that  it  is  the  noise  associated  with  measuring  the  quanti¬ 
ties  required  to  augment  a  rank  defective  CB  matrix.  The 
third  type  of  noise,  disturbance  noise,  allows  the  user  to 
add  disturbance  inputs  directly  into  the  state  equations 
in  the  form 

i  =  Ax  +  Bu  +  Gw  (A-1) 

where  w  is  a  vector  of  random  variables  representing  the 
disturbance  input  and  G  is  a  matrix  that  governs  the  dis¬ 
tribution  of  the  noise  into  the  state  equations. 

This  addition  also  provides  the  user  the  option  of 
making  multiple  simulations  to  statistically  determine  the 
influence  of  noise  through  the  use  of  a  Monte  Carlo  analysis 

2.  User ' s  Guide.  Option  25  is  selected  to  enter 
the  data  for  the  simulation  of  noise.  Prior  to  entry  into 
option  25  the  user  must  have  provided  the  number  of  states, 
outputs  and  inputs,  by  option  2,  9,  or  199.  In  this  case 
there  are  3  states,  2  inputs,  and  2  outputs.  The  following 


prompt  appears  upon  selection  of  option  25 


OPTION,  PLEASE  >  # 

?  25 

THIS  OPTION  ALLOWS  SIMULATION  OF  INDEPENDENT  GAUSSIAN 
DISTUPSANCES  AND  SENSOR  NOISE. 

ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS-. 


ENTER. SUPPRESS  OR  RESET  DISTURBANCE  INPUT . "0“ 

ENTER, SUPPRESS  OR  RESET  OUTPUT  MEASUREMENT  NOISE..."!" 
ENTER, SUPPRESS  OR  RESET  MEASUREMENT  MATRIX  NOISE... "2" 

DEFINE  MONTE  CARLO  SIMULATION . . . "3" 

TO  QUIT  OPTION  25 . "4" 

■’  0 


At  the  prompt  the  user  selects  "0"  to  operate  on  the  dis¬ 
turbance  noise.  In  this  case  the  user  desires  to  enter  new 
noise  data,  makes  the  proper  selection,  and  enters  the  data 
at  the  prompts. 


THIS  option  allows  SIMULATION  OF  A  DISTURBANCE  OF  THE 
POPM  XDOT  =  AX  +  3U  *  GW,  WHERE  W  IS  A  VECTOR  OF  N 
INDEPENDENT  GAUSSIAN  RANDOM  VARIABLES. 

3  IS  A  MATRIX  THAT  IS  N  BY  N  WHERE  N  IS  THE  NUMBER  OF 
states,  forming  a  LINEAR  COMBINATION  OF  THE  RANDOM 
VARIABLES. 

ENTER  VOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS: 


ENTER  NEW  DISTURBANCE  PARAMETERS . "0" 

SUPPRESS  DISTURBANCE  INPUT . "1" 

RESET  DISTURBANCE  INPUT . "2" 

'  0 

ENTER  THE  MEAN  AND  STANDARD  DEVIATION  OF  WC) 

n  T 

W(2) 

"  2.3 
WC) 

ENTER  THE  3  MATRIX  3Y  ROW,  3  ELEMENTS  PER  ROW 
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ENTER, SUPPRESS  OR  RESET  DISTURBANCE  INPUT . "O" 

ENTER, SUPPRESS  OR  RESET  OUTPUT  MEASUREMENT  NOISE..."!" 
ENTER, SUPPRESS  OR  RESET  MEASUREMENT  MATRIX  NOISE... "2" 

DEFINE  MONTE  CARLO  SIMULATION . "3" 

TO  QUIT  OPTION  25. . “4“ 

?  2 

THIS  OPTION  ALLOWS  SIMULATION  OF  A  NOISY  MEASUREMENT  OF 
THE  STATE  DERIVATIVES  IN  THE  CASE  OF  AN  IRREGULAR  PLANT 
THE  NOISE  IS  MODELLED  AS  INDEPENDENT  GAUSSIAN  RANDOM 
VARIABLES  WITH  MEAN  AND  VARIANCE  OF  YOUR  CHOICE  ADDED  TO 
ANY  OR  ALL  OF  THE  DERIVATIVES  OF  THE  X2  VECTOR 

ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS: 


TO  ENTER  NEW  MEASUREMENT  NOISE  PARAMETERS . "0" 

TO  SUPPRESS  MEASUREMENT  MATRIX  NOISE . "I" 

TO  RESET  MEASUREMENT  MATRIX  NOISE . 

^  0 


ENTER  THE  MEAN  AND  STANDARD  DEVIATION  OF  THE  NOISE 
ASSOCIATED  WITH  MEASURING  STATE  DERIVATIVE  1 
^  3,4 

STATE  DERIVATIVE  2? 

^  6,5 

ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OP-^IQNS: 


ENTER. SUPPRESS  OR  RESET  DISTURBANCE  INPUT . "O" 

ENTER.  SUPPRESS  OR  RESET  OU'.PUT  MEASUREMENT  NOISE..."!" 
ENTER. SUPPRESS  OR  RESET  MEASUREMENT  MATRIX  NOISE...':" 

DEFINE  MONTE  CARLO  SIMULATION . "C" 

TO  SUIT  QP"-ION  25 . '4" 

ENTER  NUMBER  OF  SIMULATION  RUNS  DESIRED  FOR  MONTE  CARLO 
ANAL  'SIS . > 


The  user  has  selected  a  Monte  Carlo  simulation  that  is  to 
be  comprised  of  five  independent  noise,  simulations.  Now, 
every  time  the  simulation  option  (option  26)  is  selected 
the  user  will  be  asked  if  that  simulation  is  to  be  included 
in  the  Monte  Carlo  analysis.  The  user  may  exit  MULTI, 
log  off,  or  run  as  many  "non-Monte  Carlo"  sim.ulatiors  as 


IM>-R164  916  NULTIVRRIRBLE  OUTPUT  CONTROL  LRN  DESIGN  FOR  THE  STOL 
(SHORT  TRKEOFF  RND  L. .  (U>  RIR  FORCE  INST  OF  TECH 
URIGHT-PRTTERSON  RFG  OH  SCHOOL  OF  ENOI. .  G  H  RCKER 
UNCLRSSIFIED  DEC  89  RFIT/’GE/ENG/89D-1  F/G  1/2 


he  chooses.  Once  five  Monte  Carlo  simulations  have  been 


executed  the  program  returns  to  normal.  To  avoid  unneces¬ 
sary  expenditure  of  computer  resources,  the  total  number  of 
Monte  Carlo  simulations  is  limited  to  twenty-five.  Finally 
the  user  decides  that  the  output  measurement  noise  should 
be  suppressed  temporarily.  This  allows  the  noise  to  be 
eliminated  without  having  to  re-enter  the  noise  parameters 
when  the  noise  input  is  required.  Entering  a  "1"  at  the 
prompt  suppresses  the  noise,  entering  a  "2"  will  reset 
suppressed  noise.  After  suppressing  the  noise  the  user 
exits  to  the  main  program. 


ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS! 


ENTER. 3IJPPRE3S  OR  PEEE’  OIST'JRSANCE  INPUT . "0" 

ENTER. SUPPRESS  OR  RESET  OUTPUT  NEASUREMENT  NOISE... 
ENTER. SUPPRESS  OR  RESET  NEASUREflE.NT  MATRIX  NOISE...'':" 

DEFINE  MONTE  CARLO  SIMULATION . '0" 

YO  our  OP’-ION  05 . "4" 

• 

"HIS  OPTION  ALLOWS  SIMULATION  OF  NOISY  OUTPUT  SENSORS. 
CORRUPTING  *HE  SIGNAL  3EING  FED  SACK.  INDEPENDENT 
SAUSSIAN  NOISE  IS  ADDED  TO  EACH  ELEMENT  OF  'HE  OUrU' 
VECTOR  WITH  MEAN  AND  STANDARD  DEVIATION  OF  'CUR  CHCIO: 

ENTER  YOUR  CHOICE  OF  THE  POL-OWING  OPTIONS: 


TO  ENTER  NEW  OUTPUT  NOISE  PARAMETERS . ’0" 

TO  SUPPRESS  OUTPUT  SENSOR  NOISE . ■' 1" 

TO  RESET  OUTPUT  SENSOR  NOISE . "C" 

1 

ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS: 


ENTER, SUPPRESS  OR  RESET  DISTURBANCE  INPUT . "0" 

ENTER, SUPPRESS  OR  RESET  OUTPUT  MEASUREMENT  NOISE... "1" 
ENTER, SUPPRESS  OR  RESET  MEASUREMENT  MATRIX  NOISE. .. ’‘2" 

DEFINE  MONTE  CARLO  SIMULATION . "Z” 

TO  QUIT  OPTION  23 . *4“ 

?  4 


At  this  point  the  user  desires  to  verify  the  inputs  he  made 
in  option  25.  This  is  accomplished  with  option  125,  which 
displays  the  current  noise  parameters.  Notice  under  the 
data  for  output  measurement  noise  the  word  "(SUPPRESSED)", 
indicating  that  this  noise  is  not  currently  being  used. 


OPTION,  PLEASE  >  4 


DISTURBANCE  NOISE  PARAMETERS 
3  MATRIX 

.lOooE+oi  .:oooE-oi  .;oooE-oi 
.:oooE*oi  .ooooE^oi  .loooE-n 
.:oooE+oi  .ioooE+01  .:oooE-oi 


NOISE  MEANS  AND  STANDARD  DEVIATIONS 
.2000E+01  .3000E+01 

.2000E*01  .3000E+01 

.DOOOE-fOl  .3000E  +  01 

OUTPUT  MEASUREMENT  NOISE 
MEAN  AND  STANDARD  DEVIATION 
.2000E+01  .3000E-01 

.lOOOE^Ol  .2000E*0l 

(SUPPRESSED) 


riEASUREflENT  MATRIX  NOISE 
MEAN  AND  STANDARD  DEVIATION 
.3000E+01  .4000E+01 

.6000E+01  .5000E+01 


3.  Proqrarmner  '  s  Guide .  The  following  FORTRAN  code 
is  located  in  PROGRAM  OPT20,  a  subprogram  of  the  executive 
program  MULTI.  This  portion  of  the  noise  option  addition 
is  an  interactive  routine  in  which  the  user  enters  the 
desired  noise,  associated  parameters,  and  the  number  of 
runs  desired  in  the  Monte  Carlo  analysis.  The  following 
variables  have  been  introduced  in  this  section  of  code: 


Variable 
WRMEAN ( I ) 

WSIGMAiI) 


G(I,J) 


DISTURB 


Description 

Vector  containing  the  means  of  each  disturbance 
to  be  added  to  the  state  equations. 

Vector  containing  standard  deviations  of  dis¬ 
turbances  . 

Matri.x  distributing  disturbances  into  state 
equation . 

Integer  flag  indicating  e.xistence  c;  iisturbance 
noise . 


PG(I, J) 


Matrix  where  G(I,J)  is  permanently  stored  when 
the  disturbance  noise  is  suppressed. 


RMEAN(I)  Vector  of  output  measurement  noise  means. 

RSIGMA(I)  Vector  of  output  measurement  noise  standard 

deviations. 

PRMN(I)  Vector  where  RMEAN(I)  is  permanently  stored 

when  the  output  measurement  noise  is  suppressed 

PSIG(I)  Vector  where  RSIGMA(I)  is  permanently  stored. 

NOISE  Integer  flag  indicating  existence  of  output 

measurement  noise. 

MRMEAN(I)  Vector  of  measurement  matrix  noise  means. 

MSIGMA(I)  Vector  of  measurement  matrix  noise  standard 

deviations. 

PMRMN(I)  Vector  where  MRMEAN(I)  permanently  stored  when 
measurement  matrix  noise  is  suppressed. 

PMSIG(I)  Vector  where  MSIGMA(I)  is  permanently  stored. 

MNOISE  Integer  flag  indicating  existence  of  measure¬ 

ment  matrix  noise. 

MONTC  Integer  indicating  number  of  Monte  Carlo  Simula 

tions  desired. 

DAT4  =80,  output  device  assignment  for  local  file 

MEM30,  the  file  that  contains  the  running  sum 
of  the  simulation  data. 

MCOUNT  Integer  counter  indicating  the  number  of 

Monte  Carlo  simulations  already  run.  This 
variable  is  only  initialized  and  stored  in  this 


«««# ftf «f f 

Ci»«««t«44*»  OPTION  25  IS  THE  NOISE  INPUT  OPTION*********************** 
C********************************************************************** 

2025  PRINT*,  'THIS  OPTION  ALLOWS  SIHULATION  OF  INDEPENDENT  3AUSSIAN’ 
PRINT*,  'DISTURBANCES  AND  SENSOR  NOISE.' 

PRINT*,  '  ' 

IF  (IFLAG(2) .EQ.O)  THEN 

PRINT*,  '#  OF  STATES,  INPUTS  Si  OUTPUTS  MISSING. .. SEE  OPTION  #2' 
80  TO  8007 
ENDIF 

5000  PRINT*,  'ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS:' 

PRINT*,  '  ' 

PRINT*,  '  ' 

PRINT*,  'ENTER, SUPPRESS  OR  RESET  DISTURBANCE  INPUT . "O'" 

PRINT*,  'ENTER, SUPPRESS  OR  RESET  OUTPUT  MEASUREMENT  NOISE..."!'" 

PRINT*,  'ENTER, SUPPRESS  OR  RESET  MEASUREMENT  MATRIX  NOISE. .. "2'" 

PRINT*,  'DEFINE  MONTE  CARLO  SIMULATION . "3'" 

PRINT*.  'TO  QUIT  OPTION  25 . ''4" ' 

READ*,  ISKIP 

IF  (ISKIP. EQ. 4)  30  TO  9007 

IFLAG(2:)al 

IF  (ISKIP. EQ.O)  THEN 

PRINT*,  'THIS  OPTION  ALLOWS  SIMULATION  OF  A  DISTURBANCE  OF  THE' 
PRINT*.  'FORM  XDOT  »  AX  ♦  BU  +  GW,  WHERE  W  IS  A  VECTOR  OF  N' 
“RINT*.  'INDEPENDENT  GAUSSIAN  RANDOM  VARIABLES.' 

PRINT*.  Q  IS  A  MATRIX  THAT  IS  N  BY  N  WHERE  N  IS  THE  NUMBER  OF' 
9RINT*,  'STATES,  FORMING  A  LINEAR  COMBINATION  OF  THE  RANDOM' 
SRINT*.  'VARIABLES.' 

SRINT*.  ■  ■ 

=RINT*,  'ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS:' 

5RINT*,  '  ' 

PRINT*.  ■  ' 

PRINT*,  ENTER  NEW  DISTURBANCE  PARAMETERS . "O"  ' 

PRINT*.  'SUPPRESS  DISTURBANCE  INPUT . "I"' 

PRI'JT*.  'RESET  DISTURBANCE  INPUT . "2'" 

READ*.  ISKIP 
IF  (ISKIP. EQ.O)  THEN 

: . --ENTER  DI3TUR9ANCE  PARAMETERS . 

PRINT*.  'ENTER  ^HE  MEAN  AND  STANDARD  DEVIATION  OF  W(l! 

READ*.  'WRMEAN(l)  .WSISMAd) 

IF  (N.EQ. 1)  60  TO  5002 
DO  5001  1=2. N 

PRINT*,  ■W(',I.')' 

500;  READ*,  WRMEAN(I)  .WSISMAd) 

5002  PRINT*.  'ENTER  THE  3  MATRIX  BY  ROW.  ',N,'  ELEMENTS  PER  ROW' 

DO  5004  1  =  1  . N 
PRINT*.  'ROW  ',! 

READ*.  (Gd  .J)  ,J  =  1  ,N) 

DO  50'03  ;  =  !,:) 
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5003  P8(I,J)»G(I,J) 

5004  CONTINUE 

OISTURBal 

BO  TO  5000 

ENDIF 

C444444444444444444444444SUPPRESS  DISTURBANCE  PARAMETERS4444444444444444 

IF  (ISKIP.EQ.l)  THEN 
DO  5007  I«1,N 
DO  5006  Jal,N 

5006  G(I,J)»0 

5007  CONTINUE 
DISTURB»0 
GO  TO  5000 

ENDIF 

C444444444444444444***444R£SET  DISTURBANCE  PARAMETERS4444444444444444444 

IF  (ISK:P.cQ.2)  then 
DO  5010  1  =  1, N 
DC  5009  J=1.N 

5009  3(I.J)»P5(I.J) 

5010  CONTINUE 
0I3TUR9=1 

ENDIF 

SO  TO  5000 
ENDIF 


IF  >  I  SKI -.ZD.;.  ’HEN 

SRINT4,  'THIS  DP’’ION  ALLOWS  SIMULATION  OF  NOISY  OUTPUT  SENSORS.' 
PRINT*.  CORRUP'^ING  *HE  SIGNAL  BEING  FED  BACK.  INDEPENDENT  ' 
=>RIN:*,  'GAUSSIAN  NOISE  IS  ADDED  TO  EACH  ELEMENT  OF  THE  OUTPUT' 
PRINT*,  'VEC’OR  AI’H  MEAN  AND  STANDARD  DEVIATION  OF  YOUR  CHOICE 
=RINT*.  ■ 

PRINT*,  ENTER  YOUR  CHOICE  OF  THE  FOLLOWING  OPTIONS:’ 

PRINT*,  ' 

PRINT*, 

PRINT*,  ”0  ENTER  NEW  OUTPUT  NOISE  ^aRAMETERS . "0" 

PRINT*.  'TO  SUPPRESS  OUTPUT  SENSOR  NOISE . '1"' 

PRINT*,  -0  RESET  OUTPUT  SENSOR  NOISE . "C" 

READ*,  ISKIP 

C******************  enter  OUTPUT  N0ISE*****4***4*4**44**************->*** 

IF  (ISKIP. EQ.C)  then 

PRINT*,  'ENTER  THE  MEAN  AND  STANDARD  DEVIATION  OF  THE  NOISE' 
PRINT*,  'ASSOCIATED  WITH  MEASURING  OUTPUT  1: 

READ*,  RMEAN(l)  .SSIGMAd) 

?RMN(  1 )  =RI1EAN(  1 ) 

PS!3(1)=RSI3MA(1) 


00  5012  I»2,P 
PRINT  *,  'OUTPUT'  , I 
READ*,  RNEANd)  .RSIGIiAd) 

PRMNd)aRI1EANd) 

PSIGdlaRSIONAd) 

5012  CONTINUE 
NOISEal 
GO  TO  5000 
ENDIF 

C*****************  SUPPRESS  OUTPUT  NOISE  ******************************* 

IF  dSKIP.EQ.l)  THEN 
DO  5014  Ial,P 
RI1EANd)aO 
5014  RSIGI1Ad)aO 
NOISEaO 
GO  TO  5000 
ENDIF 

C»«**«44***44*4««4**  RESET  OUTPUT  NOISE  *4444444444444444444444444444444 

IF  (ISKIP.EQ.O)  THEN 
DO  5016  1=1, P 
RMEANd)  =PRMNd) 

5016  RS!3riAd)=PSI3(I) 

NQISE  =  1 
ENDIF 

30  TO  5000 
ENDIF 


.  ♦•*44 


"EASUREMEN:  HATPIX  noise 


IF  (ISKIP.EQ.:)  'HEN 


PRINT*, 

PRINT*. 

PRINT*. 

PRINT*. 

’RINT*. 

PRINT*, 

PRINT*. 

PRINT*. 

PRINT*. 

PRINT*, 

PRINT*, 

PRINT*, 


THIS  OP’ICN  AL-OyS  SIMULATION  OF  A  NCI3Y  MEASUREMENT  OF  ' 
THE  STATE  OERIVA’iyES  IN  THE  CASE  OF  AN  IRRE3ULAR  =LANT ' 
"HE  NOISE  IS  .’lODELLEO  AS  INDEPENDENT  SAUSSIAN  RANDOM' 
VARIABLES  NITH  MEAN  AND  VARIANCE  OF  'OUR  CHOICE  ADDED  TO' 
ANY  OR  AL-  of  the  derivatives  OF  'ME  d  VEC'OR 


ENTER  YOUR  CHOICE  OF  'HE  pOLLOHING  OP’IONS: 


TO  ENTER  NEH  MEASUREMENT  NOISE  PARAMETERS . "0" 

TO  SUPPRESS  MEASUREMENT  MATRIX  NOISE . ''I'” 

TO  RESET  MEASUREMENT  MATRIX  NOISE . 'O’" 

READ*.  I3KIP 

»4444*444*4**44444  ENTER  MEASUREMENT  MATRIX  NO  I SE********************** 

IF  (ISKIP.EQ. 0)  THEN 

PRINT*.  'ENTER  THE  MEAN  AND  STANDARD  DEVIATION  OF  THE  NOISE' 
PRINT*.  'ASSOCIATED  WITH  MEASURING  STATE  DERIVATIVE  I' 

READ  ♦.  MRMEANd!  .  MSISMAd) 


PI1RMN(1)»I1RMEAN(1) 

PI1S1G(1)*I1SIGI1A(1) 

DO  5018  I»2,P 

PRINT#,  'STATE  DERIVATIVE  ',1,'?' 

READ#,  MRMEANd)  ,nSIGMA(I> 

PHRMN(I)»I1RMEAN(I) 

5018  PMSlGdlanSIGNAd) 

I1N0ISE>1 
GO  TO  5000 
ENDIF 

C##*#*######*###*#*#*  SUPPRESS  MEASUREMENT  NOISE  «♦♦*########*#**###**♦ 

IF  dSKIP.EQ.l)  THEN 
DO  5020  1=1, P 
MRMEANd)=0 
MSIGMAd)=0 
5020  CONTINUE 
HN0I3E=0 
GO  TO  5000 
ENDIF 

C#**##*###*######**#*  RESET  MEASUREMENT  NOISE  *♦#**♦**#♦**#♦♦♦*♦♦**#♦♦♦ 

IF  <ISKIP.£Q.2)  THEN 
DO  5022  1=1.? 

MRMEANd)  sPMRMNd) 

MSIQMAd)=PMSI3d) 

5022  CONTINUE 
MNOISE»l 
ENDIF 

30  TO  ;000 
ENDIF 

0  Dc'INE  MONTS  CARLO  ANALYSIS  44444444444444444444*44*4444 

IF  (ISKIP.E2.:)  THEN 

5023  PRINT#,  'ENTER  NUMBER  OF  SIMULATION  RUNS  DESIRED  -OR  MON'^E  CARLO' 

PRINT#,  'ANALYSIS . >' 

READ*,  MONTC 
IF  (MONTC. 3T. ;5)  THEN 

PRINT*,  'YOU  HAVE  SOT  TO  3E  JOKING.  OBVIOUSLY  'OUR  NOT  =AYING' 
PRINT*,  'PQR  -HIS.  'HE  OUTPUT  NILL  BE  ROUTED  *0  'HE  IG  ^OR  ' 
PRINT*,  FRAUD.  HASTE  5.  ABUSE  INVESTIjATIDN.  ' 

PRINT*.  ' 

30  TO  5023 

ENDIF 

DAT4=30 

OPEN  (DATi.FILE=  MEM30') 

REWIND  DATA 
MCOUNT  =  0 

WRITE  (DATA.*)  MCOUNT 
CLOSE  (DATA.STATUS='KEE?' ; 

GO  TO  5000 


The  noise  described  by  the  entries  in  option  25  is 
entered  into  the  simulation  in  option  26  by  means  of 
several  subroutine  calls  to  SUBROUTINE  GPNML  (listed  below) 
This  subroutine  generates  a  seed  and  then  makes  a  call  to 
the  IMSL  Library  routine  GGNMC  which  returns  the  random 
vector  RDEV(I).  Variables  introduced  here  are: 


Variable 

DSEED 

RDEV(I) 


NR 


Description 

Real  seed  for  IMSL  routine  GGNML. 

Vector,  random  and  zero-mean,  returned  by  GGNML 
Integer  dimension  of  RDEV(I). 


SUBRC'JTINE  3PNML  ( RME AN , RS I SMA . N . SDEV ! 

C  DEVIATES  RE’URNEO  FROM  IMSL  IN  RDEV<I) 

INTE3ER  NR 

REAL  MRMEAN.MSISMA 

COMMON  ;3  27/  MONTC.MCOUNT 

DIMENSION  RMEANdS)  .RSIGMAdS)  .RDEVdS) 

D0U3L-  PRECISION  DSEED 
DATA  DSEED  /lOOl.DO/ 

NR  =  N 

C  REFETI'IVE  CALLS  TO  3GNMC  (IMSL:  4IL.  AU’OMAT I CAL.Y  CHANGE  DSEED 

:  33NML  RE'LRNS  A  NORMALIZED  ZERO  '’EAN  3AUS3IAN  N(0.d 
DSEED  DSEED  ♦  (1000  ♦  MCOONT) 

CAL.  3GNML  (DSEED, NR, RDEV) 

DO  5C:'  1=1, N 

3025  RDEVd:  =RDEV(I)*RSI3MA(i;  *  5ME;ANd; 

:  *  =  ANSFCRM  'HE  NORMALIZED  VEC'OR  *D  N  . RMEAN . RS 1 3MA : 

R£"JRN 

END 


A  number  of  minor  changes  were  made  throughout 
MULTI  to  accommodate  the  noise  input  option.  Previously, 


the  calculation  step  size  was  entered  in  option  25.  This 
function  is  now  accomplished  in  option  24.  Option  125  now 


prints  out  the  current  values  of  the  noise  parameters. 

All  noise  data  entered  in  option  25  is  stored  in  local 
file  MEM20,  and  as  a  result  the  options  which  affect  the 
reading  and  writing  of  MEM20  (options  29,  99,  and  199) 
are  changed  accordingly.  Finally,  option  26  includes  a 
section  of  code  that  reads,  operates  on,  and  writes  to 
local  file  MEM30  to  keep  a  running  total  of  the  simulation 
data  needed  to  perform  the  Monte  Carlo  analysis.  During 
each  simulation  the  current  data  is  added  to  the  values 
stored  in  MEM30  from  previous  simulations,  creating  a 
running  total  at  each  time  increment.  When  the  last  run 
is  complete,  the  running  totals  are  divided  by  the  total 
number  of  runs  to  obtain  an  "average"  run.  This  data  can 
then  be  plotted  in  the  same  manner  as  the  results  of  any 
other  simulation.  Currently  this  code,  listed  below,  only 
calculates  the  mean  value  of  multiple  simulation  runs.  It 
is  recommended  that  in  the  future  this  be  expanded  to 
include  a  calculation  of  the  standard  deviation  as  well. 

The  new  variables  in  this  section  are; 

Variable  Description 

MONTY  Logical  character  indicating  whether  user 

wishes  current  simulation  to  be  included  in  the 

Monte  Carlo  analysis. 

DATD  =90,  input  device  assignment  for  local  file 


MEM30 . 


DAT4  =80,  output  device  assignment  for  local  file 

MEM30. 

MYP(IJ,I)  Matrix  containing  a  running  sum  of  output 
data . 

MUP(IJ,I)  Matrix  containing  a  running  sum  of  input  (U) 
data . 

MVP(IJ,I)  Matrix  containing  a  running  sum  of  input  (V) 
data . 

MXP(IJ,I)  Matrix  containing  a  running  sum  of  state  data. 

B .  Custom  Input  Option 

1.  Description .  This  option  expands  the  input 
alternatives  to  include  a  wide  variety  of  possibilities  as 
defined  by  the  user.  By  selecting  the  custom  input  feature 
of  option  22,  the  user  can  select  ten  points  that  define 
the  input  magnitude  as  a  function  of  time.  The  points  are 
connected  with  straight  lines  by  the  program  and  if  desired 
the  corners  are  smoothed.  The  option  of  using  the  original 
input  routine  is  retained  and  its  use  is  recommended  when¬ 
ever  possible,  since  it  is  easier  to  use. 

2.  User ' s  Guide.  To  select  a  custom  input,  the 
user  enters  "22"  at  the  option  prompt.  Following  is  a 
sample  of  the  interactive  prompts  and  inputs. 
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PT.  5>> 

?  9,10 
PT.  6>> 

?  11,12 

PT.  7>> 

?  13,14 
PT.  8>> 

?  15,16, 

PT.  9>> 

?  17,18 
PT.  10>> 

?  19,20 

DO  YOU  WANT  TO  SMOOTH  THE  INPUTS?  Y  OR  N 
Y 


Notice  that  if  the  user  attempts  to  enter  the  data  in  other 
than  sequential  or  chronological  order,  the  program  inter¬ 
prets  this  as  going  backward  in  time  and  requests  corrected 
data.  Like  most  of  the  data  options  in  MULTI,  the  values 
may  be  verified  in  its  corresponding  100-series  option. 

JP-;3N,  PLEASE  4 

,  initial  S’A'EE.  .  . 

.C000E*00  .0000E*00  . OOOOE-00 
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In  order  to  make  effective  use  of  the  custom  input  feature 
it  is  imperative  that  the  user  understand  the  mathematical 
foundations  of  the  smoothing  routine  and  the  assumptions 
made  in  implementing  the  option.  The  specifics  of  the 
smoothing  algorithm  are  discussed  in  the  programmer's 
guide.  Following  is  a  summary  of  features  and  limitations 
t.hat  the  user  may  find  useful. 

a.  Step  inputs  cannot  be  smoothed.  It  is 
assumed  that  if  a  smoothed  input  is  desired  a  ramp  would 
be  selected  for  the  initial  step  up  or  down.  The  program 
defines  a  step  input  as  any  two  consecutive  points  having 
the  same  time  axis  coordinate.  If  any  part  of  any  of  the 
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inputs  is  a  step,  then  none  of  the  inputs  can  be  a 
If  the  user  desires  smoothed  step  inputs  within  the  custom 
input,  it  is  recommended  that  a  ramp  with  a  duration  of 
less  than  a  sample  period  be  entered.  It  is  very  likely 
that  a  ramp  of  such  short  duration  cannot  be  smoothed  with 
the  polynomial  techniques  used,  but  even  if  unable  to  smooth 
the  step  the  algorithm  will  continue  to  smooth  the  remainder 
of  the  input  normally. 

b.  It  is  important  that  the  input  be  defined 
for  at  least  the  longest  simulation  time  anticipated.  In 
most  cases,  failure  to  do  so  will  result  in  the  value  of 
the  magnitude  of  the  last  point  being  held  throughout  the 
undefined  region.  Obviously,  points  beyond  the  simulation 
time  will  never  be  encountered  in  option  26  but  they  may 
be  useful  for  shaping  the  input  prior  to  the  end  of  the 
simulation. 

c.  All  ten  points  must  be  defined.  Note  that 
there  are  no  trivial  inputs.  If  no  input  to  a  particular 
channel  is  desired,  then  an  input  that  is  specified  as  zero 
magnitude  for  the  entire  simulation  time  is  required,  that 
is,  at  each  of  the  ten  points.  (Simply  entering  zeros  at 
both  the  time  and  magnitude  prompts  will  result  in  an  input 
which  is  only  defined  at  the  origin.) 


d.  The  input  always  begins  at  the  origin. 
Unless  a  step  is  desired,  the  time  at  point  1  should  be 


greater  than  zero.  If,  however,  the  time  at  point  1  is 
chosen  to  be  zero,  the  magnitude  should  be  non- zero. 

e.  Clever  application  of  the  mathematical 
principles  used  for  the  smoothing  algorithm  can  produce 
nearly  any  input  desired.  The  duration  and  amount  of 
smoothing  can  be  varied  without  changing  the  basic  input 
by  inserting  extra  points  along  straight  line  segments. 

A  sample  input,  both  smoothed  and  unsmoothed,  is  shown  in 
Figures  A.l  and  A. 2. 

3.  Programmer*  s  Guide .  The  bulk  of  the  code  to 
accomplish  this  feature  is  located  in  one  of  two  places — 
in  PROGRAM  OPT20  under  option  22,  and  in  PROGRAM  OPT26 . 

The  code  in  option  22  is  where  the  data  is  entered  for  the 
custom  input  feature,  and  where  the  parameters  for  the 
smoothing  curve  are  calculated.  The  basic  structure  of  the 
algorithm,  as  shown  in  Figure  A. 3,  is  to  first  establish 
the  unsmoothed,  "dot-to-dot"  input  curve.  Then,  if  smooth¬ 
ing  is  desired,  a  third  order  polynomial  is  chosen  such 
that  the  slope  and  magnitude  of  the  polynomial  match  the 
basic  curve  at  the  beginning  and  end  of  smoothing.  Smooth¬ 
ing  occurs  in  the  last  20  percent  of  the  line  segment  before 
the  point  of  interest  and  the  first  20  percent  of  the  line 
segment  following  the  point.  Often,  the  curve  to  be 
smoothed  changes  slope  too  rapidly  to  be  adequately 
smoothed  by  a  third  order  polynomial. 
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This  condition  is  indicated  by  the  second  derivative  of  the 
polynomial  at  the  beginning  and  end  of  smoothing  being 
opposite  in  sign  from  the  desired  slope  change.  This  situ¬ 
ation  can  sometimes  be  corrected  by  using  a  fifth  order 
polynomial  and  improving  two  more  conditions;  specifically 
that  the  acceleration  at  the  beginning  and  end  of  smooth¬ 
ing  be  continuous  (i.e.,  zero).  Even  fifth  order  poly¬ 
nomials  can  have  unacceptable  smoothing  characteristics  if 
the  third  derivative  at  the  beginning  or  end  of  smoothing 
is  of  opposite  sign  of  the  desired  acceleration  change. 

In  this  event  the  algorithm  attempts  to  time  scale  the 
fifth  order  polynomial  to  satisfy  these  conditions.  If 
the  user  attempts  to  smooth  an  input  that  has  large  slope 
changes  with  short  line  segments  the  smoothing  will  be 
unsatisfactory.  The  only  way  to  identify  inadequate  smooth 
ing  is  to  plot  the  inputs  (V  vector)  using  a  calcomp  plo-^ 
option  (34  or  35) .  The  astute  programmer  will  find  that 
the  routine  used  to  find  an  acceptable  time  scaling  factor 
is  an  unsophisticated,  brute  force  sequential  search.  It 
was  found  that  determining  that  the  current  value  does  not 
satisfy  the  required  conditions  yields  no  information  on 
which  direction  to  search,  rendering  a  more  efficient 
approach,  like  a  binary  search,  impossible.  If  no  solution 
is  found,  eventually  a  matrix  that  must  be  inverted  becomes 
singular,  and  the  program  returns  a  message  to  that  effect 
and  does  not  smooth  that  particular  point.  The  math  used 


to  solve  these  problems  is  not  complex,  but  is  difficult 
to  follow  from  the  code  alone.  The  key  equations  are 
developed  below  to  facilitate  understanding  of  the  program. 
The  variables  used  in  this  development  are  as  follows: 
t  independent  variable,  time 
h  dependent  variable,  input  magnitude 
tg  time  at  previous  input  point 
t^  time  at  current  point 
time  at  next  point 
t  time  smoothing  begins 
t  duration  of  smoothing 

a^  (n=0,5)  coefficients  of  smoothing  polynomial 

hQ  magnitude  of  previous  input  point 

h^  magnitude  of  current  point 

h2  magnitude  of  next  point 

s^  slope  before  current  point 

S2  slope  after  current  point 

f  time  scaling  factor 

a.  Third  order  smoothing.  If  a  third  order 
polynomial  is  to  be  used,  the  general  form  of  the  input 

/N  /N 

between  t  and  (t+t)  is: 

h(t)  =  Uq  +  a^(t-t)  +  a2(t-t)^  +  a2(t-t)^  (A-2) 

To  solve  for  the  four  unknowns  (a^,  a^^,  a2»  and  a^)  ,  four 
constraints  must  be  satisfied.  In  this  case,  the  conditions 


are  chosen  in  order  to  match  the  magnitude  and  slope  of 
the  straight-line  input  at  the  beginning  and  end  of  smooth- 
ing  (t  =  t,  and  t  =  t+t) .  Applying  these  conditions  to  the 
polynomial  and  its  derivative  (slope)  and  then  solving  the 
four  simultaneous  equations  yields: 

aQ  =  ho  +  Sj^(t-to)  (A-3) 

=  Sj^  (A-4) 

a2  =  [3  (h^  +  .2t2S2  -  .2t^S2  ~  “  a^t)  -  S2t  +  a^^t]  /t^ 

(A- 5) 

a^  =  (S2  -  -  2a2t)/3t^  (A-6) 

b.  Fifth  order  smoothing.  The  general  form  of 
the  fifth  order  polynomial  used  to  smooth  more  difficult 
inputs  is: 

h(t)  =  aQ  +  a^(t-t)  +  a2(t-t)^  +  a2(t-t)^  +  a^(t-t)^ 

+  (A-7) 

Since  there  are  now  six  unknowns,  two  more  constraints  must 
be  applied  to  find  a  unique  solution  for  each  of  the 
coefficients.  These  conditions  are  chosen  so  that  the 
second  derivative  of  the  input  is  zero  at  the  beginning 
and  end  of  smoothing.  These  constraints  result  in  the 
following  equations: 


(A-8) 


^0  =  Si(t-to) 


ai  =  s^ 


(A- 9) 


^2  =  ° 


(A-10) 


a^  +  a^t  +  a^t  =  {h^^  +  S2  [ .  2  (t2-t^)  ]  -  a^  -  a^t}/t 

(A-11) 

a^  +  (4/3)a^t+  (5/3)a5t^  =  (s^  -  a^)/3t^  {A-12) 


a^  +  2a^t  +  (20/6)a^t  =  0 


(A-13) 


For  programming  convenience,  since  a2  is  always  zero,  a^, 
a^ ,  and  a^  are  changed  to  a2»  a^/  and  a^  respectively. 
Equations  (A-11) ,  (A-12) ,  and  (A-13)  are  solved  as  simul¬ 
taneous  equations  by  MULTI  in  PROGRAM  OPT20  under  option 
#22. 

c.  Time  scaled  fifth  order  polynonial.  The 
general  form  of  the  time  scaled  fifth  order  polynomial  is 
the  S£Lme  except  for  the  independent  variable: 

h(t)  =  a^  +  aj^t  +  a2t^  +  ^3^^  ^4^^  (A-14) 

t  =  (t  -  t)/f  (A-15) 

The  conditions  of  the  fifth  order  polynomial  are  again 
applied  with  the  additional  constraint  that  the  third 
derivative  at  the  beginning  of  the  smoothing  be  the  same 
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sign  as  the  change  in  slope  desired  and  of  opposite  sign 
at  the  end  of  smooth,  ng.  The  program  calculates  the  poly¬ 
nomial  coefficients  and  executes  an  iterative  search  for 
a  value  of  "f”  that  will  satisfy  these  constraints. 

d.  FORTRAN  code.  The  equations  presented  in 
the  last  three  paragraphs  are  only  the  basic  framework  for 
the  custom  input  routine.  The  programmer  will  also  notice 
a  number  of  conditional  statements  in  both  option  22  and 
option  26  to  avoid  overflow  conditions  that  result  from 
dividing  by  zero  and  other  discontinuities.  The  variables 
introduced  in  these  sections  of  code  are  contained  in  the 
comment  statement  preceding  option  22. 


El-  OP'I 

NP, EO. 

O'i  'HEN 

HRIN't, 

"HIE 

’OR'ION  OF  -HE 

PROGRAM  A,..0W5  *H£  USER 

’PIN'-* . 

-"1  •'3 

'w  w  - 

-INE  lO  ’01  NTS 

ALONG  A  CUSTOM  INPUT  'H 

HRIN'*  . 

AWE  ■ 

0  'SE  OONNEO'Ej 

3Y  S'RAIGH'  .ONES  AND 

PBIN'*  . 
00  1 320 

•  "nEN 

smgothe:  if  so 

DESIRED. 

PRINT*. 
00  .35: 

'EN'ER 
I  = .  I 

INPUT  ' .K. ■ ; 

time.  MAGNI'UDE 

=  ‘^INT*, 

■  J  T  ^ 

.  I  . 

:f  :'ip<=’  k .  - . :  ^  4rjD.  - .  j£. : '  '‘'“Ien 

=  ~  I  ^.E  -'IP  ,1.  jPEA'ZP  'ihN 

■>P:Nr*,  OP  EQUAL  'C  'RV  AGAI'J,  ■ 

30  '0  1321 

£nd:p 

•.  350  O'JNT'.NUt 

:3i:  JON'INUE 

PRINT*.  00  'OU  NANT  "0  SMOOTH  THE  INPUTS'’  ^  OR  N ' 

'READ  ■  A)  .  SMOP’ 

0  ”-!£  OOOE  'HAT  ■CL-O'.^S  OOMPUTES  THE  OOErPICIENTS  OF  'HE 

:  ’O.  NOMIAlS  'HA'  ARE  JSEO  AS  SMOOTHING  OURVES  3E’'H£EN  ..'NE  SEOMEN 
0  IN  EACH  INPUT.  THE  ’CLTNOMIALS  ARE  INITIAL-'^  OHCEEN  AS  'HIRE  IRC 


189 


POLYNOMIALS  SUCH  THAT  THEIR  MAGNITUDE  AND  SLOPE  MATCH  THE  LINE  ♦ 
SEGMENT  values  AT  THE  POINTS  WHERE  THE  SMOOTHING  STARTS  AND  STOPS.  * 
OFTEN,  HOWEVER,  THE  INPUT  CHANGES  SLOPE  TOO  RAPIDLY  TO  BE  SMOOTHED  * 
BY  A  THIRD  ORDER  POLYNOMIAL.  MATHEMATICALLY  IT  IS  POSSIBLE  TO  MEET  * 
ANY  BOUNDARY  CONDITIONS  (OTHER  THAN  INFINITE  SLOPES)  WITH  ONLY  A  ♦ 
THIRD  ORDER  POLYNOMIAL  BUT  THE  CURVE  SOMETIMES  INITIALLY  TURNS  IN  * 
THE  WRONG  DIRECTION.  THIS  PHENOMENON  IS  EVIDENT  IN  THE  SECOND  ♦ 
DERIVATIVE  OF  THE  FUNCTION,  WHICH  SHOULD  AT  LEAST  BE  THE  SAME  SIGN  ♦ 
AS  THE  change  in  SLOPE  AT  THE  POINT  OF  INTEREST.  THE  CODE  TES^^S  FOR* 
THIS  CQNDI'ION  AND  CALCULATES  COEFFICIENTS  FOR  FIFTH  ORDER  PQL'^S  ♦ 
IF  NECCE3ARY.  THIS  ALLOWS  TWO  MORE  CONDITIONS  TO  BE  IMPOSED  ON  THE  ♦ 
SMCCTHING  CURVE,  AND  ’’HEY  ARE  CHOSEN  SUCH  THAT  THE  ACCELERA’-IQN  IS  ♦ 

OONTINUCUS  and  zero  at  the  start  and  stop  of  the  smoothing.  I'  IS  ♦ 

=0S3I3LE  IN  EX-'REMELY  DIFFICULT  SMCCTHING  SITUATIONS  (LIKE  LARGE  ♦ 
SLOPE  changes  in  very  short  time  »ERI0DS)  that  the  same  problem  WIL-* 
ARISE  IN  -IE  third  derivative  (JERK)  OF  THE  FIRTH  ORDER  POL'^NCMIAL.  ♦ 
"HIS  CONDI'ION  IS  TESTED  AS  WEL-  AND  IF  NEOOESSARY  THE  FIF’H  ORDER  ♦ 
ROL'NOMIA;.  is  TIh£  scaled  3Y  a  -AC'GR  SINCE  THE  SCALE  paCTOR  ♦ 
SESLL'S  IN  non-line, AR  SIMUL^’ANEOUS  EDUATICNS  A  SOLUTION  (NC*  UNIOUE)* 
IS  -u'JND  ■(■'ROUGH  AN  I'ERATIVE  SEARCH.  A  SOLUTION  MAY  NOT  EXIST  IN  ♦ 
*^E  REGION  SEARCHED  (  =  .01  ;  * 
IN  ’HIS  E'.'ENT  THAT  ^ARTICULAR  -OINT  IS  NOT  SMCO’HED.  NOTE  THAT  AS  ♦ 
■HE  ORDER  OF  THE  POL'^NQMIAL  INCREASES  IT  APPROACHES  A  ’AYLDR  * 
SERIES  RE=PE5EN^ATICN  OF  *H£  INP'U*.  anD  -ESS  SMOOTHING  OCO-RS. 

AP'-R  SMCD'HING  THE  INPUT  SEGMENTS  ARE  STRAIGHT  IN  THE  MIDDLE  rl’'.  ’ 
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Hi,,. 

.  .  .  IMpP* 

Y 

•  •  '• 

.  MAG.  A-  : 

URRENT 

INP'C  P-. 

HZ... 

. . .  IN?  =  T 

■  ,  * 

; • .  mag,  a ' 

NEC 

INPUT  P'. 

i  1  .  .  . 

. . . SLOPE 

wP 

■0 

CURRENT  input  p’. 

5Z .  .  . 

. 3-ZPE 

AR' 

E  R 

C-RRE.N*  INP 

'  j  T  ^  , 

TH.  .  . 

...  'I. "IE 

SMCC 

H  * 

NG  S'AR'S. 

H  ,  .  , 

. . . DURA* 

ION 

J  * 

SMOOTHING 

<DC"H.  .  . 

. . . value 

IC 

5MCCTHE0  ACCEL 

EPATION  AT  ■= ' -w-  - 

CDT.  .  . 

.  .  . VAL-E 

smoothed  jerk 

AT  T=( 

TH*'HH 

fT 

. . . ’IME 

3CAL 

ING 

-AC'GR  P 

.  o: 

AC.AA.  ,  . 

. . . INPPT 

■■  K. ,  I 

,  J  ! 

,  ro 

L'NOMIAL  CCETPICIE 
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INPPT{K,I,B) 


C 

C 

C 

C 


0, 

3, 

5, 

f, 


WHEN  UNABLE  TO  SMOOTH  ♦ 
WHEN  THIRD  ORDER  SMOOTHING  USED  « 
WHEN  FIFTH  ORDER  SMOOTHING  USED  ♦ 
WHEN  FIFTH  ORDER  IS  TIME  SCALED  * 


IF  (SMOPT.EQ.  'Y')  THEN 
DO  185S  K=1,P 
DO  1054  1=1,9 
IF  (I.EQ.l)  THEN 
H0=0 


T0  =  0 


ELSE 

HO=INP»T(K, 1-1.2) 

TO=INPPT(K, 1-1,1) 

ENDIF 

H1=INPPT (K, I ,2) 

H2=INPPT(K. 1+1,2) 

Tl=rNPPT :k, I .1) 

’:=IMPPT ;k, :+i , 1) 

IF  (T1.EQ.T2,QR.T0.EQ.T1)  THEN 

PRINT*,  'YOU  HAVE  A  STEP  IN  INPUT  '  .  tC ,  '  THAT' 

PRINT*,  'CANNOT  3E  SMOOTHED.  DO  YOU  WISH  TO  <A>BCRT' 
PRINT*.  'THE  SMOOTHING  ROUTINE,  OR  <E;'NTER  NE'W 
PRINT*,  'INPUT  DATA'!’  ENTER  <A>  OR  <£>' 

READ ' (A) ' ,3MCPT 
IF(SMOP’.EQ. 'A';  THEN 
3M0?’’=  N' 

30  *0  :S55 

E-3EIP  ■  S^OP".  10.  E')  "'HEN 
30  *0  :c;: 

E-3E 

30  '0  ;o:: 

ENDIF 

EL3£ 

3'.  = 'HI -HO:  "'.-’01 
32= ;H2-H1 ; / (T;-ti  ; 

ENDIF 

'‘^=,3*T‘.  - 
^HH=,2*T2  -  ,0*T0 
A0=H0+51*f’H-T0) 

Al=31 

a:=  hi  ♦  .;♦';*=:  -  .;*’'i*s:  -  ao  -  thh*a' 
a;=  a2*3  -  32  +  ‘'+H  ♦  a:*''’hh 
A2=  A2/"+^**2 

A3=(S2  -  A1  -  2*A2*'^HH)  /  (3*THH**2) 


<DDTH=2*A2*3*A3*ThH 
INPO' (K. I , 3) =3 

^  -NO. a:. 0)  inp^t (k. i , 3' =5 


iF(Sl  .L'I'.O.  AND.  A2.GT.0)  INPPT  (K  ,  I  ,  Si  =5 
IF(S2.GT.O. AND.XDDTH.LT.O)  INPPT ( K , I , 8) =5 
IF(S2.LT.0.AND. XDDTH.GT.O)  I NPPT ( K , I , 8 ) =5 
IF(INPPT(K,I,8).EQ.5)  THEN 
PQLYMAT<1,1)=1. 

PQLYMATd  ,2>=THH 
PQLYMAT ( 1 ,3) =THH**2 
PQLYMAT(2,1)=1. 

PQLYMAT(2,2)=(4. /3. )*THH 
PQLYIAT (2,3) =(5. /3. )*THH#*2 
P0LYI1AT(3,l)  =  l. 

PQLYMAT(3,2)=2.«THH 
PQLYMAT(3,3)=(20./6.)*THH**2 
HH=H1+S2* (THH-Tl+TH) 

CA2='HH-A1*THH-A0) /THH**3 
CA:=(32-A1)/(3*THH**2) 

CA4=0, 

CAL.  INVEPKPQlYMAT,  IPaLY,3,3,»i:61) 

A2  =  CA2*  IPQLY ( 1 , 1) +CA3*IPOLY ( 1 . 2 ) +CA4* I  POL Y (1.3) 

A:=CA2.:?QLY(2,l)>CA:#IPQLY(2.2)*CA4*IPaL'd2.3> 
A4=CA2*IPaLY(3,l)-CA3*IPQLY(:.2)-CA4*I?0LY(:,3) 
<T0Tsi.*A2+24. ♦A3*THH*oO. ♦A4*THH**2 
:nP°’C<.  :.3)=5 

32-31 ) .  LT.  0.  AND.  A2.  3T.0)  :NP®’’;K.I.3)=1 
IF  '  'S:-3i; .3T.0. AND. AC.lT.O)  INP®'  : K . 1 , 3) =1 
IF  'A2. 3T. 0. AND. XTOT, 3T. 0)  !NP»’ ( K ,  1 . 3 ) = i 
IF  A2.L'.  0.  AND.  XTDT.LT.O)  I NP®’’ ( K .  1 .  3 )  s  1 
IF  ;np®’-.k.i.8).ed.S)  then 
;nP®^;k.;.:)=ao 
:np®’':k,:.4!sai 
INP®T  :  .5)  =A2 
IN?®*{K.:.j;=43 
I  NP®’  '  K  .  I  .  " ) =A4 
J '«  ■  J  .  -  2  . 

2  -  '  £ 

30  *0  1222 
END  IF 

INP®'  K  .  I  .  3)  =0 

PPINT*.  IJNABLE  TO  3N0C'H  INP'JT  '.n,'  AT  =■".  '  I 
INP®'^  (K  .  I  ,3)  =A0 
INP®’ (K. I . 4) =Al 
■F  ' INP®’  < , : . 3) . EC. 1 i  ’HEN 
F  =  .  01 
®  I  =0 

mh=h;-32«(THH-T1-'H) 

2At  =  '  '^H-41*THH  /  F-AO'  /  (THH/F  **Z 
CA3=i3:-Ai' ; ;3*(THH/F)**2) 

I  H  =  0  . 
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P0LYI1AT(1,2)=THH/F 
P0LYMAT(l,3)a(THH/F)**2 
P0LYI1AT(2, 1)  =  1. 

P0LYMAT(2,2)a(4./3.)*THH/F 

PQLYMAT(2,3)s(5./3.)*(THH/F)**2 

P0LYf1AT(3,nal. 

P0LYI1AT(3,2)»2.*THH/F 
PQLYI1AT<3,3)s(10./3.  )*(THH/F)**2 
CALL  INVERT(POLY«AT, IP0LY,3,3,*126l) 
A2=CA2*IPaLY(l,l)*CA3*IP0LY(l,2)-fCA4*IPQLY(l , 
A3=CA2*IP0LY(2,1)+CA3*IPQLY(2,2)+CA4*IP0LY(2, 
A4sCA2*IPQLY(3, l)+CA3*IP0LY(3,2)+CA4*IP0LY(3,3) 
XTDT=6. *A2+24.*A3*(THH/F)+60. *A4*(THH/F) **2 

TP  A  Akir  AO  OT  A\  PT-l 


IF  (Fl.ca.l)  THEN 
r=P*1.02 
30  TO  1263 
cL3E 

INPPT(K.I.:)=AO 
INPPTfK. I .4)=Al 
INPOT  (K. : .5) *A2 
INpoT  :k.  : .  j) *43 

;np°t  ;k. : ,7) =A4 
INPOTCK.I.ai^F 
■NDIF 
END  IF 
EL3E 

INP®T  'K. I.:)=A0 
iNP°';K,;.ii=Ai 
INPOT  ;k. 1 .5) =A2 
:np»-(k.:.3/=a3 
iNP‘’T(K,:,r)=o. 
INPPT{K,I,3)=3 
ENDIF 
CONTINUE 
CONTINUE 
ENDIF 
ENDIF 

iF'.AG.:::!  =  i 
30  -0  3007 


C-l 


f.3  fj* 


The  following  code  is  located  in  PROGRAM  OPT26 


DO  1861  K=1,P 

C »*44 *4 «**«*«««**« 4 *4 

C  THIS  IS  THE  CODE  WHERE  THE  INPUT  IS  GENERATED  WHEN  4 

C  A  CUSTOM  INPUT  HAS  BEEN  SELECTED.  AO , A1 , A2 , A3 , A4  ARE  4 

C  COEFFICIENTS  OF  THE  THIRD  OR  FIFTH  ORDER  POLYNOMIAL  USED  4 

C  TO  SMOOTH  THE  CURVE.  TO. Tl. HO. HI  ARE  THE  TIMES  AND  4 

C  MAGNITUDES  AT  THE  BEGINNING  AND  END  OF  THE  LINE  SEGMENT  4 

C  RESPECTIVELY  AND  31  ;=A1,'  13  "HE  SLOPE  OF  THE  LINE.  ♦ 

C  THE  VALUE  CONTAINED  IN  INPPT(K.I.0)  DETERMINES  WHETHER  * 

THE  CURVE  IS  SMOOTHED  BY  A  THIRD  ORDER  POLY,  FIFTH  ORDER  4 
POLY,  OR  A  TIME  SCALED  FIFTH  ORDER  POLY.  OPTION  *22  4 

CONTAINS  A  DETAILED  DESCRIPTION  OF  THE  SMOOTHING  ME'HODS  4 
AND  THE  definition  OF  THE  SMOOTHING  VARIABLES  IN  A  CCMMENT4 
STATEMENT  PRIOR  "0  THE  SMOOTHING  ALGORITHM.  * 

4444'*4*444444444'*4444*4«**«*444««*4*««*****4«**44444*'*«4444 

IF  (SMOPT.EQ. 'N'!  THEN 


1401  IF.I.^T.IO)  'HEN 

IF  ;;.gT.  IINPPT  :K.  ; .  1,  ■ ;  'HEN 
1*1  +  1 

30  TO  1401 
ENDI" 


END  IF 

IF  'I.EC.D  'HEN 
IF'INPP'^K.I.I;. 

^iKlsINPP’’:-;,! 

ELSE 

3i  =  i  ?ot;:<.;.:: 

ViKl*Sl*T 

ENDIF 

ELSE 

IF  f  (  INPO'  I  .i: 

V  'to  =  iNpp'  ;k.  : . 


iNpp  ■ 


'  H  E  N 


sisiNpp"  .  I . - 1 Npp’’ ' 

3i=si,  ' iNF='  -M , . ; -:nf=' 

J  :  t; s  :  N  p  s  T  ;  ^  •  .  ;  _  ;  .  .  5  I  »  '  -  I  N  P  °  T  :  ,  T  .  1  ,  ;  ;  , 

END  I  F 
END  IF 
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IF(I.LT.IO)  THEN 
IF  (T.GE. (.a*INPPT{K,I,l)+.2*INPPT(K,I+l,l)))  THEN 
1*1  +  1 

GO  TO  1402 
ENDIF 
ELSE 
1  =  10 
ENDIF 

IF  {I.EQ.D  THEN 
T0=0 
H0  =  0 
ELSE 

TO=INPPT(K, 1-1,1) 

H0=INPPT{K, I-l ,2) 

ENDIF 

T1=INPPT (K, I , 1 ) 

TH=. 3*T1+, 2*T0 
H1=INPPT(K,I,2) 

A0=:nPoT(K,I,3) 

A1=!NPBT (K, I ,4) 
fl2=INPPT{i(,I,5) 

A.:=!NPPT  ;k,  1 ,4) 

A4stNPPT(K,I,7' 
if;:.3t.th.and.i.le.9)  then 
IP  :npo''k,:.3) .S3.:)  then 

V  ('<:  =A0  +  A1«  ■,’-*H)-A:t  (T-TH)  ♦♦2*A.:t  (T-TH) 

SL3E:F:INPoT;K,:.3).S3.3)  then 

V  (k;  =ao+Ai*(T-TH)  *a:*(t-th)  ♦♦3+A:*a-TH)  **4 

+A4*(T-ph)**5 

SL3E;F(T.3E.T1.AND.  INPP'^fK.I.B)  .E3.0)  THEN 
32=INPP'^  (K.  I  +  l  .2)  -INPPT.K.  1 ,2) 

32  =  32/  ( INPPT  ;k.  1-1 .  l ) - 1  NPof ( K , 1 . 1 ) ) 

V!K)=H1  ♦  52+(T-T:) 

EL3E:F  (INPPT(K,:.3).NE.O)  THEN 
p=1./INPPT:K,:.3) 

y  (K) =A0  +  A1»  iT-TH) •P-A2* ' ( T-TH ) ♦?  '  +*3 
+  A3*  (  (T-tH) ♦?' t  +  A  +  AA* ( (T-TH) ♦F) *4  = 

ENDIF 

EL3E 

•J'K)  =H0*AU  (T-TO) 

ENDIF 

NDIF 

CNTINUE 


C.  Suppression  of  Actuators 
and  Sensors 

1.  Description.  It  is  sometimes  convenient  to 
eliminate  actuator  and  or  sensor  dynamics  from  the  simula¬ 
tion.  Previously,  this  would  involve  destroying  the 
actuator  and  sensor  data  and  then  re-entering  the  same  data 
when  the  dynaimics  are  desired.  The  actuator  and  sensor 
data  is  now  stored  in  permanent  variable  locations,  while 
temporary  variables  can  be  suppressed  and  reset  in  option 

4  or  option  5. 

2.  User ' s  Guide .  The  interactive  prompts  are  self 
explanatory  for  this  change  and  are  listed  below. 

■:p”:cn.  please  >  * 

* 

■'-ilS  ENTERS  *HE  ACTUATOR  STATE  EQUATION  DATA 

ENTER  '0"  *C  SUPPRESS  AC’UATORS 

en'er  *0  SET  ac'uator  values...'' 
ente=  ’O"  *0  USE  store:  actuator  values...; 

■  ) 

OPTION.  PLEASE  •'  # 

-HIS  OPTION  ENTERS  THE  SENSOR  STATE  EQUATION  DATA 

ENTER  "O"  TO  SUPPRESS  SENSORS 

ENTER  '1”  TO  SET  SENSOR  VALUES... > 

ENTER  '0"  "0  USE  STORED  SENSOR  VALUES...) 


It  is  important  to  note  that  if  the  actuator  or  sensor 
dynamics  are  suppressed,  they  will  not  be  saved  in  MEMO 
when  exiting  the  program.  A  warning  message  to  this  effect 
has  been  added  to  the  exit  routine  and  is  shown  in  the 
"Saving  Memory  Files  Without  Exiting"  section  of  this 
appendix. 

3.  Programmer*  s  Guide.  The  code  to  accomplish 
this  option  is  very  simple  but  is  spread  out  in  options  4, 

5,  9,  99,  104,  and  105.  For  these  two  reasons  it  is  not 
repeated  here.  The  following  variables  are  added  for  this 
feature: 

Variable  Description 

PNA(I)  Vector  of  "m"  integers  (m  =  number  of  inputs), 

each  being  the  number  of  states  in  the  actuator 
for  that  input.  This  variable  is  a  permanent 
storage  location  for  the  vector  variable  NA(I), 
the  quantity  used  by  the  simulation  for  actuator 
state  data.  NA(I)  is  set  to  zero  when  the 
actuators  are  suppressed  and  is  set  equal  to 
PNA(I)  when  the  actuators  are  reset. 

PNS(I)  Vector  analagous  to  PNA(I)  containing  the  num¬ 

ber  of  states  for  each  output  sensor.  NS (I) 
is  the  local  variable  used  by  the  simulation. 


D.  Saving  Memory  Files  Without  Exit 

1.  Description.  Upon  selection  of  option  99, 

MULTI  will  save  all  pertinent  data  in  local  files  MEMO, 
MEMIO ,  MEM20,  and  MEM30  and  then  the  program  will  either 
return  to  normal  execution  or  exit  according  to  the  user's 
desires. 

2.  User' s  Guide .  Option  99  allows  graceful 
termination  of  MULTI  and  automatically  saves  all  plant, 
actuator,  sensor,  design  and  simulation  data  in  local  files 
prior  to  exiting  the  program.  However,  as  all  MULTI  users 
will  inevitably  discover,  there  are  a  number  of  ways  to 
exit  MULTI  involuntarily,  leaving  the  user  with  the 
irritating  task  of  re-entering  all  data  that  had  not  been 
saved.  The  most  commonly  encountered  inadvertent  termina¬ 
tion  of  MULTI  occurs  when  the  user  enters  a  "RETURN"  at  the 
prompt  without  any  data  preceding  the  "RETURN".  The  com¬ 
puter  program  has  interrupted  execution  at  a  read  state¬ 
ment  and  is  expecting  input  from  the  terminal.  If  no  input 
±s  provided,  an  "END  OF  FILE"  is  encountered  and  the  pro¬ 
gram  aborts  execution.  Naturally,  this  phenomenon  is 
accompanied  by  the  loss  of  all  volatile  data,  which  may 
have  taken  hours  to  generate.  This  problem  has  not  been 
corrected,  but  if  the  user  is  cautious  to  save  data 
regularly,  the  frustration  of  re-entering  data  can  be 
avoided  and  one  is  likely  to  stay  motivated  toward  the 
ultimate  objective  for  a  considerably  longer  time.  The 


procedure  for  saving  data  is  quite  simple  as  is  demon¬ 
strated  below: 


OPTION,  PLEASE  >  # 

?  99 

ALL  PLANT  INPUT  DATA  HAS  BEEN  SAVED  IN  A  LOCAL  FILE 
CALLED  "MEMO" 

ALL  DESIGN  DATA  HAS  BEEN  SAVED  IN  A  LOCAL  FILE 
CALLED  "MEMIO" 

ALL  SIMULATION  DATA  HAS  BEEN  SAVED  IN  A  LOCAL  FILE 
CALLED  "MEM20'' 

ALL  MONTE  CARLO  SIMULATION  DATA  HAS  BEEN  SAVED 
IN  A  LOCAL  file  CALLED  "MEMCO" 

»*4***«**««'»4*4*4*«4**44*«*»«**«*'»«*4»4»*4 

4  ACTUATORS  AND  SENSCRS  WERF.  NOT  SAVED  * 

444444444444444444444444*44444444444444444 

DO  YCU  HANT  TO  EXIT  MULTI:  Y  OR  N  ? 

•'  N 

OPTION.  PLEASE  >  # 


3.  Programmer ' s  Guide .  The  code  changes  required 
to  accomplish  the  desired  changes  to  option  99  consist  of 
several  conditionals  to  determine  whether  sensors  and 
actuators  have  been  suppressed,  generation  of  a  warning 
based  on  that  determination  and  finally  a  question  asking 
the  user  if  termination  is  desired.  The  only  variable 
introduced  is  a  logical  character  "EXIT",  depending  on  the 
user's  desires.  The  exit  routine  code  reads  as  follows: 


c 


. ROUTINQ  FOR  OPTION  #99 . 

ELSEIF  (N0PT.EQ.99)  THEN 
IF  (IPLOT.GT.O)  THEN 
CALL  PLOTE  (BLK) 

PRINT' (A/)  ' ,  '  REMINDER:  ROUTE  "PLOT(S)"  BEFORE  LOGOUT 
ENDIF 

CALL  OVERLAY  (MULTI, 12,0) 

IF  (ACT.EQ. 'N' .OR.SEN.EQ. 'N' )  THEN 
PRINT*,  '*««*««««*4***««»*«*«*****»«««*»«**«««»««»* ' 

IF  (ACT.EQ.  'Y')  THEN 

PRINT*,  '*  NOTE:  SENSOR  DYNAMICS  WERE  NOT  SAVED  *' 
ELSEIF(SEN.EQ. 'Y')  THEN 

PRINT*,  '*  NOTE:  ACTUATOR  DYNAMICS  WERE  NOT  SAVED  *' 
ELSE 

PRINT*,  '*  ACTUATORS  AND  SENSORS  WERE  NOT  SAVED  *' 
ENDIF 

PRINT*,  '»*****************444***«*4****«*******4*«' 

ENDIF 

PRINT*,  'DO  YOU  WANT  TO  EXIT  MULTI:  Y  OR  N 
READ  • (A)  ■  ,  EXIT 
IF  (EXIT-EI]'.  'N' )  THEN 
GO  TO  9000 
ENDIF 

PRINT' (A/) ' ,  '  HAVE  A  NICE  DAY' ' 

STOP 


E.  Convert  Input  Vector  "U" 

From  Radians  to  Degrees 

1.  Description.  After  completion  of  option  26 
the  user  is  given  the  option  of  converting  several  of  the 
data  arrays  from  radians  to  degrees  prior  to  plotting  the 
data.  Previously  this  option  did  not  include  the  control 
input  vector  "U" .  The  routine  now  includes  this  conversion 
as  well,  to  account  for  plants  in  which  the  input  matrix 
is  given  in  terms  of  radians.  The  original  code  is  the 
work  of  Major  Terry  L.  Courtheyn  (5:C-1).  Courtheyn's  work 
is  merely  copied  to  accomplish  the  additions,.,  conversion. 


2.  User' s  Guide .  The  prompts  for  this  option  are 
identical  to  the  prompts  originally  programmed  by  Courtheyn 
with  the  addition  of  a  similar  prompt  for  the  conversion 
of  the  "U"  vector.  Both  the  use  and  programming  of  this 
change  are  self-explanatory  and  the  programmer's  guide  is 
omitted. 

F.  Plot  Combination  of  States 
and  Inputs 

1.  Description.  Often  it  is  desired  to  plot  not 
only  a  state  but  its  derivative  as  well.  In  the  case  of  an 
aircraft,  it  is  often  convenient  to  plot  the  normal  accelera 
tion  as  a  function  of  time,  requiring  a  combination  of 
states  and  state  derivatives.  Since  in  a  linear  system  of 
equations  any  state  derivative  can  be  described  in  terms  of 
the  states  and  inputs,  all  that  is  required  is  to  be  able 

to  combine  state  and  input  data  to  obtain  any  function  of 
states  and  state  derivatives  as  a  function  of  time.  This 
change  expands  the  existing  capability  of  plotting  combina¬ 
tions  of  states  to  the  option  of  plotting  a  user  definable 
combination  of  states  and  inputs. 

2.  User ' s  Guide .  Following  is  the  interactive 
dialog  that  the  user  will  encounter  after  selecting  one  of 
the  six  plotting  options  (31-36) .  This  particular  example 
is  a  terminal  plot  option  (31).  To  obtain  a  plot  of  some 
combination  of  states  and  inputs,  the  user  selects  plot 
choice  "4"  at  the  prompt. 


THIS  OPTION  PRODUCES  A  PLOT  AT  YOUR  TERMINAL 
PLEASE  CHOOSE  ONE  OF  THE  FOLLOWING: 


FOR  A  SINGLE  SAMPLING  TIME 

1.. .A  PLOT  OF  UP  TO  2  INPUT  AND  OUTPUT  PAIRS 

2.. . A  PLOT  OF  UP  TO  4  INPUTS  OR  OUTPUTS  OR  STATES 

3.. . A  PLOT  OF  UP  TO  4  DIFFERENT  SIMULATIONS 

(FOR  ANY  SINGLE  INPUT  OR  OUTPUT) 

OR  4... A  PLOT  OF  UP  TO  4  COMBINATIONS  OF  STATES  ' 


ENTER  CHOICE  DESIRED  > 
’  4 


CHOICE  44... YOU'VE  CHOSEN  TO  PLOT  COMBINATION  OF  STATES 
ENTER  THE  NUMBER  OF  COMBINATIONS 
OF  STATES  AND  INPUTS . > 


ENTER  "I"  MATRi:(...i  ROWS  WITH  4  ELEMENTS  EACH 
ROW  1  > 

•  «*«■•«** 

COMBO  MATRIX  C... 

.lOOOE-01  .COOOE-01  .3000E-01  .iC00£*01 


Ij  "HIS  C3R"EC'‘...''E3.N0.J...> 


Up  to  this  point,  the  user  sees  no  change  in  the  inter¬ 
active  prompts.  Now  the  program  requires  entry  of  the 
inputs  to  be  included  in  the  combination.  As  with  states, 
the  user  enters  the  matrix  which  adds  the  weighted  inputs 


into  the  desired  combination. 


CURVE  X  ABOVE  IS  COMBO  1 


DO  YOU  WANT  A  LIST  OF  POINTS  USED  IN  PLOTTING"^ 

ENTER. . .YES  OR  NO.  .  .  > 

^  N 

:?"IQN,  PLEASE  >  » 

3.  Programmer ' s  Guide .  Ali  of  the  code  to  iMOyduce 
pLjta,  either  at  the  terminal  or  files  for  CALCOMP  opjcting, 
IS  located  in  three  overlays:  OPTPLOT,  OPT31,  and  34. 
Altliough  there  are  six  options  (31-36)  that  require  the 
conujining  of  states  and  inputs  for  plotting,  there  is  only 
one  routine  to  accomplish  the  calculations  and  it  i.i  located 
in  ’r'TPTjOT.  The  code  to  combine  the  inputs  is  nearl” 
ideniL',;ai  to  the  code  to  combine  the  states  that  was 
OIL  i.niLLy  in  .MULTI.  This  code,  Listed  below,  requi.;'.-  r.he 
add.-.  :.on  of  two  variables. 

’/  .  c;l^  De  script  ion 

Array  containing  the  coefficients  used  ■  - 
combine  (I)  inputs  into  (J)  combination::. 

IIMl'  Logical  character  indicating  the  presei-.'-' 

inputs  in  the  combination  tc.'  be  plotted. 


1406  PRINT  '(/A)','  ENTER  THE  NUMBER  OF  COMBINATIONS 

PRINT»,  'OF  STATES  U  INPUTS . )  ' 

READ*,  K 
DO  1  499  1  =  1, K 
DO  1  497  J=1  ,M 

1497  Z i I , J) =0 
DO  1499  L=1 , INPS 

1498  ZU ( I ,L) =0 
1495  CONTINUE 

PRINT*,' 

PRINT*,  'ENTER  '' Z  "  M  A  TR  I  X  .  .  .  '  ,  K  ,  '  ROWS  WITH  ■,M,'  ELEMENTS 
DO  1  490  1  =  1, K 
PRINT*,  'ROW  ' , I ,  '  >  ' 

1  490  READ*  ,  {  Z  (  I  ,  J  )  ,  J  =  1 ,M) 

P5INT*,' 

SO  TO  190 

14R5  CALL  FIX  (Z  .K,M) 
l^O  5=;,'JT*,  'COMBO  MATRIX  Z.  .  .  ' 

:al_  matpr  (Z ,k,m) 

CALL  ANSWER  (*1RR2, *9010) 

5  =  I  N  T  *  ,  ■ 

5RINT*,  •  does  THE  OOMSO  INCLUDE  INPUTS' . Y 

"AD  ■  !  A)  ■  ,  I  INP 
;  5  :  INP.  ES.  ■  r*  '  )  THEN 

so  TO  i5o: 


y  OR  N 


30  ■^0  1303 


>«««,'•«  j'j  xATsr;  "HC  COMBO  OF  INPUTS  ♦♦♦♦♦♦»♦♦♦■*•* 

««««*>4«*«t«««*44***t*******«*4**4*«********T*********4.r 

3:;  -iNT*, 

=  -I'lT*,  ENTER  'Z'J"  MATr;;;...  •,k.  R0'.*S  WITH  ,INP=,  £ 

0 :  4  0 1  i  = '. .  A 

■ "  I  ■)  r  * ,  R  0  'w  ,  I , 

.  -5 :  .2* ,  ;  ZU  ‘  I ,  j '  ,  J  =  1 .  INPS) 

■-  I'l:*, 


1  A  5  _  I  . _  -  I  <  (  Z  U  .  A  ,  I  N  P  S  ) 

15'  -=;nt*,  ■  COMBO  MA'RIX  ;u. 

: , _  <1 A :  ?  R  : :  u ,  ' . ;  N  p  s  I 

: ANSWER  (  *  1  AQC  ,  *‘'0  ' 


1  joi  s  =  s  +  X  V I  ,L+n  *1 ;  J ,L) 
1505  PLMAT(I,J)=S 
1510  CONTINUE 
ICLN=K 

CHQICE='  COMBO' 


;i5  IF  (IFLTCN.EQ.O)  IFLTCN=1 
IFLTCN=IFLTCN-1 
IF  (  IFLTCN. NE. 0)  GO  TO  1320 
IF  (ICQDE.EQ.3)  ICLM=LINES 
GO  TO  aOlO 


'  j  .  1 


:iiii  1 .1 1  ion  of  lion  1  inea cities 
'  " ‘  ^  to  Aircra  f  t 

1.  Description  .  Linear  mode  is  of  aircraft  a 
[uj.te  accurat«9,  provided,  of  co'urse,  tlie  assLii; 
'  .  ide  in  ‘ibtaining  the  linear  m(jd(el  are  not  vkdI  ■ 
■  cLmuiation.  One  of  the  key  assumptions  is  that 
;  ::u!  mfiments  on  the  aircraft  are  linear  with  cfui 

1' ■  f  1  i.'o  t :  on  .  If  control  surface  >ieflections  ir 

is  i.n  maMimum  performance  maneuvers,  ijr  in  tine 
c’ntiy  nonlinear  control  surfaces  like  vectoriui 
ji  -■  tiirust,  a  Linear  model  is  inadequate.  For  t‘c. 
ioscribt^d  in  Chapter  III  the  principal  nonline 
:  ■  1  t'n'.;  Ltudinai  control  surfacu;  deflecti'ins  is  ■ 

;.;i  of  the  siqn  of  the  partial  derivative  of  vel. 

■  to  the  deflection  when  tiie  surtace  passe.; 

.•.ero  anqie  of  attack.  This  phenomenon  is  easi 
:  in  the  simulation  and  is  implemented  in  a  S[)e’  ■ 
j !:  ML'LTI  customized  for  aircraft  inodels.  Ln  • 
"Ctitorc’d,  vari.ible  thrust  a:ul  nou  i  .:.:u}a  i:  ,i  l  le  s  i 


more  complex.  A  rigorous  development  of  the  nonlinear 
effects  of  this  type  of  input  is  contained  in  Reference  15. 
These  effects  are  also  simulated,  at  the  user's  option, 
in  the  customized  version  of  MULTI. 

2.  User ' s  Guide .  Unfortunately,  to  implement 
these  two  features  it  is  necessary  to  place  additional 
requirements  on  the  allowable  form  of  the  model  used  in 
the  customized  version  of  MULTI.  These  constraints  ajc  as 
follows : 

a.  The  plant  must  be  longitudinal,  body  axis, 
linearized  model  of  an  aircraft. 

b.  The  states  are  defined  as  THETA  (pitch 
angle) ,  U  (X-axis  velocity) ,  Q  (pitch  rate) ,  and  ALPHA 
(angle  of  attack) .  These  are  all  perturbation  values  and 
must  be  arranged  in  that  order.  Additional  states  arc 
allowed  but  must  be  after  ALPHA  in  the  state  vector. 

c.  The  first  two  inputs  must  be  aerodynamic 
surfaces,  like  stabilators  or  canards.  The  third  input 
must  be  a  variable  thrust  input  like  a  throttle  or  revorser 
vanes.  If  a  two-dimensional  nozzle  is  desired,  its  deflec¬ 
tion  angle  must  be  the  fourth  input. 

d.  The  equilibrium  angle  of  attack  for  each 

of  the  aerodynamic  surfaces,  and  the  equilibrium  deflection 
of  the  two-dimensional  nozzle  must  be  known  and  entered  in 
option  #  3 . 
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It  is  important  to  note  that  when  a  two-diinunsionaL 


nozzle  is  used,  not  only  is  the  simulation  nonlinear,  but 
the  solution  is  not  unique.  This  of  course  means  that 
there  are  an  infinite  number  of  steady  state  solutions  and 
it  may  be  difficult  to  obtain  one  that  is  satisfactory.  It 
is  left  to  the  user  to  determine  how  one  finds  a  satisfac¬ 
tory  solution.  A  sample  of  the  interactive  prompts  for 
option  #3  follows: 


OPTION,  PLEASE  >  * 

N 

E'lTEP  EQUILIBRIUM  VALUE  FOR  EACH  INPUT 
INPUT  I 
^  )1 
INPUT  : 

INPUT  : 

-HERE  A  -'HC-jIMENS:  jNAL  NOZZLE  INPUT  ON  'HIE 
POP  AFT-  ;y  or  N) 

■I  I 

E'trZR  THE  NOZZLE  MOMENT  ARM  FROM  03  .'FTi' 

ENTER  =IT0:h  moment  OF  INERTIA  lYY  (SLUG-*P'««:! 

-  lEsiHSP 

ENTER  THE  AIRCRAFT  MASS  (SLUGS)  ) 

.  1  —  1 

•  to 

ENTER  THE  EQUILIBRIUM  VELOCITY  (FT, SEE:  ■ 

EN-ER  THE  DERIVATIVE  Z-ALPHA-DOT  (F'"  SEE:  • 

-  ;  . EOS 


OPTION,  PLEASE  '■  * 


3.  Programmer  '  s  Gu ide .  There  arc  two  bloi;,ts  ol 
added  to  MULTI  for  this  feature.  The  first  bl is 
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located  in  PROGRAM  OPTO  under  option  #3  and  is  the  inter¬ 
active  portion  where  the  user  enters  the  necessary  data. 
The  second  block  of  code  is  located  in  PROGRAM  OPT26  and 
is  where  the  nonlinearities  are  actually  computed  during 
the  simulation. 


Variables  Description 

EV(I)  Vector  containing  the  equilibrium  angles  of 

attack  for  aerodynamic  surfaces  and  initial 
nozzle  and  thrust  input  values. 

EVA (I)  Vector  containing  time  varying  angles  of 

attack  for  the  aerodynamic  surfaces. 

N0Z2D  Logical  character  indicating  presence  of  a 

two-dimensional  nozzle. 

LX  Local  real  variable,  nozzle  moment  arm. 

lYY  Local  real  variable,  aircraft  pitch  moment  of 

inertia . 

MASS  Local  real  variable,  aircraft  mass. 

UEQ  Local  real  variable,  equilibrium  velocity. 

ZAD  Local  real  variable,  body  axis  coefficient  of 

force  in  the  z  direction  with  respect  to  the 
time  derivative  of  the  angle  of  attack. 

BNOZl  Real  variable,  nonlinear  input  matrix  coeffi¬ 

cient. 

BN0Z2  Real  variable,  nonlinear  input  matrix  coeffi¬ 


cient. 
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-=:nT*,  'ENTER  EQUILIBRIUM  VALUE  FOR  EACH  INPUT' 

:  a  :  1 1  i  =  i .  .i 

—  'NT*,  'INPUT 
=  EAD».  EVd) 


31 i  ::ntinue 


^^INT*,  ' 

13  THERE  A 

t'nq-dihensional  ncz:le  inpu'"  jn  this 

. 

-:nt*,  ' 

AIRCRAFT"'  1 

Y  OR  N)  • 

-EAO  '(Ai 

.nq::d 

'NozcD 

.EQ.  Y'l  THEN 

-'PINT*. 

'ENTER  THE 

NOZZLE  HCMENT  arm  -pom  03  =-*1  ' 

5EAD» . 

L  A 

4 

-=:NTr, 

-:nte5  =:t 

CH  MOMENT  OF  INER'IA  I'/"' 

t 

'  -  L  . 

•  7  7 

3no::=l 

;  / 1  y  y 

4 

INTf. 

'ENTER  THE 

AIPCRAF'  MASS  ,3L'-33, 

=  EAC»  , 

1A33 

-  P  I  N  T  <  , 

ENTER  THE 

EQUILIBRIUM  VELOCI"/  F'  SEC: 

4 

='EAD», 

UEQ 

‘  P  1  N  T  ♦  . 

ENTER  THE 

OERI'VATIVE  Z-ALRHA-OOT  .F'.'SEZ:  ■' 

=  EAD»  , 

:ao 

1 

3No::=i 

,  /  'HASS*  fUEQ-ZAD) ) 

£-  i£ 

Na::D= 

N  ' 

$ 

£NDIF 

:-LAG(3) =1 
;=lhnt=i 
33  ’■Q  3001 


The  remaining  code  for  this  feature  is  located  in  PROGRAJi 
OPT26 ; 


CODE  TO  HARDWIRE  NQNL I NE AR I T I E3  FOR  STOL  F-15 

EVA(n=EV(l)+X  (4) 

E7A(2)=E7(:)+X(4) 

EVA(3) =EV(3) 

IFiNOZDD.EQ. 'Y')  EVA(4)=EV(4) 

DO  1  204  1  =  1, N 
SfiL  ( I  .  3)  =S  ( I  ,3) 

DO  1203  J  =  1  ,2 
3NL( I .J) =B( I  ,  J) 

!-  EVA( J) .GE.O)  THEM 

IF-'  (U(J1  +  EVA(J)  )  .LT.O)  THEN 
BNt.;:,J)=-3(2,  J) 

EMDIF 
E  -  3E 

IF!  (U(J)  +  EVA(J) ) .3T.0)  THEN 
BNL!:,J)=-3(;,J) 

E'IDIF 

E'IDIF 

?:  lOMTIiNUc 

;-i  :::it:noe 

I-  nq:;d.eo.  y' ;  them 
ENL  : , =8 ( : . 3) *9No:  1 ■  3) 

3)  =3(4.:)-3N0Z:»X';:) 

END  IF 

»  +  4**»*44-*  N0M-LIMEARITIE3  44444444444444- 


>444444444444 


;i.  Ca  Lata  Initial  In  tegrator 
3  ^  ?  Vec  tor  ZO 

L,  Description .  MULTI  requires  two  vectors  o  i: 
Initial  --onditions  to  specify  an  initial  system  state. 

The  fii  .C:  is  the  initial  conditions  desired  on  each  of  :ha 
:;Lant  acates.  In  the  case  of  an  aircraft  this  specifies 
tlie  initial  orientation  and  motion  of  the  aircraft.  The 
second  ’•■".•tor  is  the  initial  conditions  imposed  by  th>.' 


inti_"jril  of  the  error  vector  Z(0).  If  initial  control 
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surfaces  deflections  are  desired  they  must  be  specified  b' 
the  Z|0)  v'ector.  The  relationship  between  Z(0)  and  U(0) 
is  given  by: 

Z  (0)  =  g  Kl~^  U  (0) 

g  =  forward  loop  gain  (1/SAMPT) 

K1  =  controller  integral  gain  matrix 
[J(0)  =  initial  control  surface  deflections 

Option  46  now  includes  a  routine  that  calculates  the  Z(0) 
vector  using  the  current  values  of  Kl  and  g  as  well  as  a 
user  specified  U(0). 

1.  User's  Guide .  This  feature  is  invoked  by 
selecting  option  6  and  making  the  appropriate  choice  from 
the  menu  (shown  in  the  example  below) .  The  program 
requests  the  desired  initial  control  inputs,  calculates  t! 
Z(0)  vector  and  stores  it  in  the  appropriate  memory  locat 

'h:2  :z^!P'jte3  the  transfer  functions  of  the  syste;^ 

FOR  :FE')-.J0P  'ransfer  .-unction  enter  1, 

-OR  :t:£E:-'.3CP  transfer  '^unction  enter  2  > 

FOR  j  1  -ND  5-0)  INVERSE  NATRICES  ENTER  0  '■ 

*0  :.-L:,jt^'E  I'O)  VEC'OR  ENTER  A  -■ 

^  1 

En:e=  ■'■'E  :  I.EiENTS  OF  THE  DESIRED  U(0)  VECTOR 

THE  INI'I-l  OONDI’ION  10(1)  HAS  SEEN  SET  TO: 

szi.  itessssses  -;oo. /loesssasR  49.21010101011 


OP'ION, 


^tSASE 


4 


mmmmmmmmm 


3.  Programmer  s  Guide .  The  code  for  this  featut 
resides  in  PROGRAM  XFERFN,  the  overlay  for  computing  the 
system  transfer  functions. 


362  PRINT' (/A)',  '  FOR  OPEN-LOOP  TRANSFER  FUNCTION  ENTER  1,' 
PRINT*,  'FOR  CLOSED-LOOP  TRANSFER  FUNCTION  ENTER  2  >  ' 
PRINT*,  'FOR  G(0)  AND  G(0)  INVERSE  MATRICES  ENTER  3  >  ' 
PRINT*,  'TO  CALCULATE  ZtOt  VECTOR  ENTER  4  >' 

READ*,  TFTVPE 


-t-t-****-******* 


*•***«*«* 


C*  ’'HIS  SECTION  CALCULATES  THE  NECCESARY  INITIAL  CONDITION  VEC'C-  Z' 
C»  ■rC  SPECIFY  A  VECTOR  OF  INITIAL  CONTROL  INPUTS  U;0).  THE  ECJATIjN 
C*  use:  “OR  THIS  CALCULATION  IS: 

:♦  ;jiO)  =  G*Ki*:(0)  +  g*ko*zdot(0) 

C*  ASSUMING  THAT  THESE  INITIAL  CONDITQNS  ARE  IMPOSED  TO  ESTABLISH  AN 
C«  ESUILIBRIUM  CONDITION  AITH  NQN-IERO  CONTROL  SURFACE  DEFLZC'IDNS. 
;d;*'0)  =  D,  AND  THE  EQUATION  REDUCES  TO: 

U (0)  =  G*KI*Z (0) 

Z'O)  =  (K1  INVERSE)*U(0)/G 

********* 

'TFTYPS. £Q, 4)  THEN 

=  R:nT*,  'ENTER  I’HE  ',.1.'  ELEMENTG  OF  'HE  DESIRED  U(V' 

'E.AO* ,  u  M) 


:AL-  INVERT'Vl  I  ,P.PD,  *4S69) 

:c 

: ' ' : : =0, 

;C  43 iT.  J=l,? 

:  )  <  I  !  =10  I  I ) *K1 I f I . J) *U( j) ♦SAMPT 

iontinue 

MNTINUE 

=  =’:NT*,  THE  INITIAL  CONDITION  ZO>::  HAS  SEEN  SET  rC; 
■PINT  *  ,  '  :o (I )  . I  =  !  .?' 

SO  TO  iS'O 

=RINT*.  K1  MATRIX  IS  SINGULAR  AND  CANNOT  3E  INVERTED 


m 


e.  SUBROUTINE  INVERT.  This  subroutine  formats 
a  matrix  and  its  associated  parameters  for  inversion  by 
the  IMSL  routine,  LINV2F. 

f.  SUBROUTINE  FIX.  If,  in  SUBROUTINE  ANSWER, 
the  user  desires  to  change  a  matrix,  this  subroutine 
accepts  the  changes  and  updates  the  matrix. 

g.  PROGRAM  OPTO.  This  overlay  contains  the 
routines  for  the  plant  input  options  (options  0  through  ‘J)  . 
However,  option  6  is  a  separate  overlay  called  XFERFN . 

h.  PROGRAM  OPTlO.  This  overlay  contains  all 
the  design  parameter  routines  (options  10  through  19)  with 
the  exception  of  options  14  and  18. 

i.  PROGRAM  OPT14U.  OPT14U  calculates  the 
controller  matrices  KO  and  K1  for  designs  in  which  the 
plant  parameters  are  unknown. 

j.  PROGRAM  OPT14R.  KO  and  Kl  are  calculated 
in  OPTL4R  for  regular  plants  (first  Markov  parameter  non¬ 
zero)  . 

k.  PROGRAM  OPT14I.  Irregular  plant  controller 
matrices  are  calculated  in  this  overlay. 

l.  PROGRAM  OPT18.  In  the  case  of  an  irrctjular 
plant,  a  measurement  matrix  is  required.  Option  18  (con¬ 
tained  in  OPT18)  provides  several  utility  routines  that 
can  bo  useful  in  choosing  an  appropriate  measurement 


matrix . 


m.  PROGRAM  OPT20.  With  the  exception  of 
options  2b  and  28,  all  simulation  options  (20  through  29) 
are  contained  in  OPT20 . 

n.  PROGRAM  OPTPLT.  OPTPLT  is  the  first  of 
four  routines  (three  overlays  and  a  subroutine)  written  to 
generate  plots.  OPTPLT  is  the  interactive  portion  in  whicli 
the  user  selects  the  type  of  plot  and  the  necessary  param¬ 
eters  (options  30  through  39). 

o.  PROGRAM  OPT31,  Upon  selection  of  one  of  thi 
terminal  plot  options  (31-33) ,  OPT31  interactively  asks 

for  data  specifically  required  for  terminal  plots.  The 
data  is  tlien  formatted  for  use  by  the  terminal  plot  sub¬ 
routine  PI.OTIT. 

p.  SUBROUTINE  PLOTIT.  This  subroutine  is  an 
adaptatioti  of  the  generalized  routine  used  to  produce  plots 
on  tlio  1  i  no  printer.  It  produces  a  plot  at  the  user's 
terminal  using  non-graphics  characters. 

q.  PROGRAM  OPT34 .  OPT34  transforms  the  data 
for  plotting  into  the  form  required  by  the  CALCOMP  plotting 
routine . 

r.  PROGRAM  ERROR.  This  overlay  contains  all 
messages  t)iat  result  from  errors  that  are  neither  fatal 
nor  terminal.  These  errors  are  usually  a  result  of  attem^jt 
ing  to  perform  calculations  requiring  data  that  has  net  yet 


been  enter  e-d. 


s.  PROGRAM  MEMORY.  Upon  selection  of  option 
99,  EXEC  routes  execution  to  MEMORY  for  generation  of 
memory  files  MEMO,  MEMIO ,  MEM20 .  Section  3  contains  the 
format  of  the  files  generated. 

t.  PROGRAM  PRINT.  PRINT  contains  all  of  the 
100  series  options  that  print  the  current  values  of  the 
data  generated  in  any  of  the  input  options. 

u.  PROGRAM  OPT14B.  This  overlay  computes  the 
controller  matrices  when  the  BSTAR  method  is  chosen  in 
option  14. 

V.  PROGRAM  XFERFN.  This  is  the  overlay  that 
executes  option  6.  This  option  includes  computation  of  any 
open  or  closed  loop  transfer  function,  steady  state  transfer 
functions,  and  initial  controller  integrator  states. 

w.  SUBROUTINE  PHOFS .  This  siJoroutine,  called 
by  XFERl'll,  calculates  the  transfer  function  denominator 
polynomial s . 

X.  SUBROUTINE  CADJB .  CADJB  is  also  called  by 
XFERFN  and  computes  the  transfer  function  numerator  poly¬ 
nomials. 

y.  SUBROUTINE  POLYRT.  POLYRT  calculates  the 
roots  of  the  polynomials  generated  by  PHOFS  and  CADJB. 

z.  SUBROUTINE  CLMAT.  This  subroutine  calcu¬ 
lates  tlu’  closed  loop  matrix  used  by  XFERFN  to  compute 
the  closi'il  loop  transfer  functions. 
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a  a . 


PROGR/U^  OPT26.  OPT26  is  the  overlay  that 


performs  the  simulation.  It  is  important  to  note  that  the 
simulation  integrates  one  calculation  step  at  a  time,  allow¬ 
ing  the  introduction  of  noise,  nonlinear  effects  like  con¬ 
trol  surface  saturation,  and  data  packing  for  plotting. 

bb.  SUBROUTINE  CLPASS .  CLPASS  is  the  first  of 
four  subroutines  called  by  OPT26  to  form  the  differential 
equations  rmior  to  invoking  the  library  routine  ODE  to  solve 
them.  CLPASS  is  used  to  form  the  equations  when  both 
actuator  and  sensor  dynamics  are  present. 

cc.  SUBROUTINE  CLPSSl.  CLPSSl  is  used  to  form 
the  dif fcnintial  equations  when  only  actuator  dynamics  are 
present . 

dd.  SUBROUTINE  CLPSS2.  CLPSS2  is  used  to  form 
tlie  differential  equations  when  only  sensor  dynamics  are 
pre  sen  t . 

ee.  SUBROUTINE  CLPSS3.  CLPSS3  is  used  to  form 
the  differential  equations  when  neither  actuator  nor  sensor 
dynamics  are  present. 

ff.  SUBROUTINE  GPNML.  This  subroutine  uses 
the  IMSL  library  routine  GGNML  to  produce  a  zero  mean, 
gaussia-.  ..andom  vector  with  a  standard  deviation  of  1. 

GPNML  uses  this  normalized  random  vector  to  obtain  a  random 
vector  with  the  mean  and  standard  deviation  required  by  tlie 
various  noise  inputs. 
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gg.  SUBROUTINE  YOUT.  This  subroutine,  called 
by  OPTIG ,  computes  the  output  vector  from  tlie  state  vi-ctor 
and  "C"  matrix. 

hh .  PROGRAM  OPT28 .  OPT28  executes  the  figures 
of  merit  calculations  of  option  28.  The  figures  of  merit 
are  based  solely  on  the  empirical  data  calculated  during 
the  sunulation.  No  theoretical  techniques  such  as  the' 
LaPl.ice  final  value  theorem  have  been  implemented.  Option 
28  be  executed  only  once  for  each  simulation,  aft'.-r 

whit.'h  li  I  figure  of  merit  data  is  inaccessible. 

3.  Memory  Files.  MULTI  generates  four  local  imunory 
files  t..)  prevent  the  user  from  having  to  enter  all  tlv,: 
requi.'oi  data  for  each  e.xecution.  These  files  have  i  spe- 
cifi.-  •'crTnat  that  must  be  maintained  if  the  user  chouses 
t;  Miu'.aally  rreate  or  edit  the  data.  The  user  sliculd  bi} 
a'.va;-'  tiiar.  any  file  he  intends  to  use  must  be  a  local  fi,  be 
pr:'.'!  '.o  entry  into  the  MULTI  pro<jram.  Following  are 

exam:  i  ■  of  each  of  the  data  files. 

a.  MF.MO  .  This  file  '.'ont  ains  the  plant, 
ac  *. 'a  i;  ,,r  and  sensor  data. 


.  3, 

■“1  - .  Ic; -y' i; <5 1  -40. '.so'’  .315'.^ 

O  O  c  u  3  f  0  -  .  ?  ®  I  3  5  1  .  3  a  0  8  9 

-.o39;i3 

I  , 
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Y  .  n  tor  :  i  i  n  t  :  i 

3,  Ij’ 

4  H  at  r  ;  ( 


9  1  a  t  Y- ;  ( 


V*.  ill 


C  rt  a  t  r  1  :< 


-1 . 


Is  thers  a  D  Matri:(' 

D  M  a  t  r  1  ■< 

Equilibrium  surtacs  Dssitiqn  ♦♦ 
Is  thera  a  2-d  nozrls' 

N  02:13  at-fect  on  3  'i:  ♦♦ 

A  r  j  c  n  a  r  j  a  c j  a  t  o  r  : .  ■  j  ,t.  ;  i  ’ 

^  3  :  a  1 3  3  in  3  a  c  n  i  c  ;  r 
Actuator  1:  A  latr; 

A  ■;  t  J  a  t  0  r  1  :  3  1  a  t ■ 

Actuator  1  ;  C  Mat"'. 

Actuator  2;  A,  E, 
inO  Z  lat'ioaa 


Actuator  I:  A,  3. 

1  n  0  2  1  a : '  1 : a 
A  r  a  n  3  r  a  3  3  n  a  o  r  :  -  ;  :  i  ' 

3  i i  1 0  3  .1  }  a  c  n  5  a '  ■  :  ■ 

;  3  r  :  c  r  .  ;  A  ,1  a  t  . 

i  3  ^  5  •:  r  ;  ;  3  t  i  t  ' 

:  U  0  r  .  ;  ;  1  i  r  ■ 

i^nsf  2:  A.  3, 
an:  2  :i  a  t "  1 ;  3  3 


oanso'"  1;  A,  3. 
an:  2  1  s  t ' :  :  o  o 


11':  nos  ire 


found  only 


in  ML2M0 


filos  intondoi 


;  ■ .  L‘ 


I  vo-rsion  dosiqnod  for  .in'oraft. 


MrL'r 


b.  MEMlO .  This  file  contains  the  design  data. 


I 

.01 

15.  .4  .4 
.  05 

-.000444036307  .0411717192  .252760573 
.000301990294  -.0279978440  .0992128160 
-.0344416265  .00205523127  -.0261815939 
-.0000044408630  .0004117171  .0025276057 
.0000030199029  -.0002799794  .0008921281 
-.00034441626  .00002055231  -.0002619159 
0. 

C 

0. 

0.  1 . 

0  .  0 . 

1 .  0 . 


).  0. 

.  25  1 


Type  oT  design  (R,I,'J,3) 

Alpha 

Sigma  matrix  diagonal  element 
Epsilon 

KO  Matrix 


K 1  .Matrix 


Measurement  Matri 


F  Matrix 


c.  MEM20  .  This  file  contains  the  simuiati  .ii 


da  ta . 


) .  1 .  ' .  ) . 

) .  : .  1 . 

0 .  3 .  1 .  3 . 
0.  3. 

0.  0.  3.  0. 
3.  3. 


20 . 

.  025 


0 .  3 .  3 . 

0  .  0  .  ; . 

0.  3.  3.  0. 

).  3.  0. 

. 


initiai  state  'ector 
I  n  1 1 ;  .a  1  i  .n  t  e  q  r  a  t  o  r  ^  -m  : :  r  3 
Casto.m  or  Stanoaro 


Data  ‘or  sta  .rl 


DontroJ  aur-ace 

Sample  Time 

Total  simulation  tine 

Cai cui at  1  on  atep  3i : a 

Computational  delav 

Output  noise  meana 

Output  noise  stanoars  seviat: 

Disturbance  noise  near,  s 

Distjroance  loiae 

Measurement  noise  nea.ie 


\  JV 


yii4)iiiill|jMUiiw*MM 


mmfci'iiailii^iiii'ita  miiatJuiaaaiL 


0. 

I 

0 . 

Measurament  noisa  da/ii 

0. 

0. 

0 . 

0. 

0. 

0. 

0. 

0. 

Distarbanca  Noise 

0. 

0. 

0. 

0. 

G  Matrix 

0. 

0. 

0. 

0. 

1 

1  1 

Noise  flags 

d.  MEM30.  This  file  contains  the  data  us'vl  i.n 
the  Monte  Carlo  noise  simulation.  It  is  recommended  tha  r. 
the  user  not  tamper  with  this  file  since  it  contains  a 
great  deil  of  raw  plot  data  with  little  apparent  meanir.g. 
Thus,  .\v.  example  is  not  shown  here. 


I  I 

^  ■ 


K 
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Appendix  B:  Multivariable  Control  Theory 
(Edited  and  reproduced  from  Reference  2) 

This  thesis  uses  the  multivariable  design  method 
of  Professor  Brian  Porter  of  the  University  of  Salford, 
England  (13).  The  design  method  employs  output  feedback 
with  high-gain  error-actuated  controllers.  Output  feedback 
is  advantageous  since  state  variables  may  be  difficult  to 
measure  while  system  response  data  are  more  readily  avail¬ 
able  . 

System  State  Equations 

Porter's  method  works  equally  well  for  either  con¬ 
tinuous  or  discrete  systems,  but  it  is  often  easier  to 
first  excimine  a  system  in  the  continuous  time  domain. 

This  is  because  of  the  numerical  accuracy  problem  with 
designing  in  the  z-plane.  A  continuous  time  system  is 
represented  by  the  state  space  model; 

X  =  Ax  +  Bu 

y  =  Cx  (B-1) 


where 

A  =  continuous  plant  matrix  (n  x  n) 

B  =  continuous  input  control  matrix  (nx  n) 
C  =  continuous  output  matrix  (Hxn) 


X  =  state  variable  vector  with  n  states 


u  =  input  vector  with  m  inputs 
^  =  output  vector  with  £  outputs 

The  system  inputs  for  an  aircraft  are  the  control  deflec¬ 
tions  or  actuator  input  commands,  and  the  system  outputs 
are  aircraft  responses  affected  by  the  inputs. 

The  method  does  not  allow  for  a  feedforward,  D, 
matrix.  If  such  a  matrix  is  present  in  the  original  state 
space  model,  the  control  inputs  must  be  redefined  as  states 
so  that  the  D  matrix  is  absorbed  into  the  C  matrix.  This 
can  be  accomplished  by  incorporating  the  actuator  dynamics 
into  the  plant  model.  Actuator  inputs  then  become  control 
inputs. 

To  employ  Porter's  method,  it  is  desirable  (but  not 
necessary)  to  partition  the  system  state  equations  as 
follows: 

+ 

(B-2b) 

The  equations  are  partitioned  so  that  B2  and  £2  are  square 
(m  x  m)  and  (2,  x  £)  matrices,  respectively.  The  method 
requires  that  the  number  of  inputs  to  the  system  equals 


^1 


u 


(B-2a) 


the  number  of  outputs  which  means  m=£,  and  therefore  the 


dimension  of  equals  the  dimension  of  C2'  is  always 

possible  to  form  the  state  equations  so  that  Some¬ 

times,  however,  a  transformation  matrix  T  is  necessary  to 
achieve  [0,  B2]  form.  In  this  case,  the  transformed  states 
no  longer  have  the  same  physical  significance  that  the 
original  states  once  had. 

For  the  discrete  case  the  system  equations  are 
written  as  follows: 

x[(k+l)T]  =  ^(kT)  +  jjm(kT) 

y(kT)  =  rx(kT)  (B-3) 


where 


0  =  exp (AT)  =  discrete  plant  matrix 


=  f  exp(AT)Bdt  =  discrete  input  control  matrix 


£  =  C  =  discrete  output  matrix 


In  the  above  equations  T  is  the  sampling  period,  and  k 
takes  on  integer  values  from  zero  to  plus  infinity. 

System  With  Output  Feedback 

Figure  B.l  shows  the  block  diagram  for  a  continuous 
output  feedback  system,  where  v  is  the  command  input  vector 
and  y  is  the  desired  output  vector.  The  blocks  for  the 
plant  are  derived  directly  from  the  system  state  equations. 


Equation  (B-1) .  The  proportional  plus  integral  controller 
has  three  parameters,  and  g,  which  must  be  deter¬ 

mined  by  the  designer.  The  output  signal  of  the  con¬ 
troller,  u,  is  given  in  the  following  control  law  equa¬ 
tion: 

M  (B-4) 

where 

u  is  the  output  signal  of  the  controller 
e  is  the  error  signal  at  the  input  of  the  controller 
Kq  is  the  proportional  gain  matrix 

is  the  gain  matrix  for  the  integral  term 
g  is  the  scalar  forward  path  gain 

Figure  B.l  is  the  depiction  of  a  system  with  only  first- 
order  integration  in  the  controller  design.  The  theory 
allows  for  a  q-dimensional  bank  of  integrators  in  which  case 
the  controller  is  made  up  of  (q  +  1)  K  matrices,  thru 
K  .  A  measurement  matrix  M  is  included  in  the  system  if 
the  plant  is  irregular.  Regular  and  irregular  plants  are 
discussed  later. 

The  discrete  system  block  diagram,  shown  in 
Figure  B.2,  is  similar  to  the  continuous  system,  but  Equa¬ 
tion  (B-4)  becomes 

u(kT)  =  (1/T)  [K^eCkT)  +Kj^z(kT)]  (B-5) 
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where  the  forward  path  gain  g  equals  the  sampling  frequency 
(1/T) .  The  £(kT)  matrix  is  derived  from  the  backward  dif¬ 
ference  equation, 

2[(k+l)T]  =  z(kT)  +  Te(kT)  (B-6) 

The  steps  to  be  taken  next  in  the  design  method 
depend  on  whether  or  not  the  first  Markov  parameter  [CB] , 
has  full  rank,  i.e.,  does  it  have  an  inverse.  If  the 
matrix  [CB]  has  full  rank,  the  plant  is  called  "regular" 
and  no  measurement  matrix  M  is  needed.  However,  if  [CB] 
does  not  have  full  rank,  the  plant  is  called  "irregular" 
and  M  is  needed  to  form  a  new  matrix  [FB]  (See  Equations 
(B-12)  through  (B-14))  which  does  have  an  inverse.  This 
is  explained  in  more  detail  in  the  next  sections.  When 
the  partitioned  B  matrix  in  Equation  (B-2a)  has  the  form 


(B-7) 


then 


[CB] 

=  ICjBjl 

(B-8) 

[FB] 

(B-9) 

As  in  the  continuous  case,  a  q-dimensional  bank  of 


integrators  applies  equally  well  to  the  discrete  design 
(Figure  B. 2) . 


As  the  gain  factor  of  the  system,  g  (or  1/T  for  the 
discrete  case) ,  approaches  infinity,  the  system  transfer 
function  matrix  G{s)  assiames  the  asymptotic  form 

r(X)  =  f(X)  +  r(X)  (B-10) 

where 

£(X)  is  the  slow  transfer  function  matrix 

£(X)  is  the  fast  transfer  function  matrix 

The  roots  of  the  asymptotic  closed-loop  transfer 
function  may  be  grouped  into  three  sets:  ^3* 

Table  B,1  gives  the  equations  for  finding  these  asymptotic 
roots.  Sets  and  correspond  to  the  slow  modes  of  the 
system,  where  the  modes  associated  with  the  roots  in 
become  uncontrollable,  and,  for  regular  designs,  the  modes 
associated  with  the  roots  in  Z^  become  unobservable  as  the 
gain  increases.  Set  Z ^ ,  the  infinite  roots,  are  associ¬ 
ated  with  the  fast  modes  of  the  system  which  become 
dominant  as  the  gain  increases. 

The  roots  in  set  ^2  correspond  to  the  transmission 
zeros  of  the  system  which  are  not  altered  by  output  feed¬ 
back.  As  the  gain  is  increased,  the  closed- loop  roots  of 
the  system  tend  to  migrate  toward  the  transmission  zeros. 
This  may  adversely  affect  the  system  stability  if  the 
location  of  these  zeros  is  in  the  unstable  region. 


TABLE  B.l 


ASYMPTOTIC  EQUATIONS  FOR  ZERO-B2  FORM 


■itii  represented  by: 


-11 

^12 

1 

il 

22 

s 

-21 

i  -22 

1  ^ 

-2 

Continuous  Case 
(s-plane) 

Gain  Factor  »  g 


Discrete  Case 
1  r-piane ) 

Gain  Factor  •  1/ 


C„(XI„  -  A.)  -3- 
— Q  — n  —0  —u 


r  ( X) 


(XI 


i„  *  gc,3-:<^)  -gc.a.x-. 


■'X)  •  (x: 


—r\  — Ti 


:a  Soots 


C;^iSo  *  iiii 
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*.rta  Scots 


=  01 


t '  : 
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Secular  Desicn 


Irregular  Desicn 
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Reference  15  gives  a  procedure  for  locating  the  transmis¬ 
sion  zeros  of  a  system. 


For  a  regular  design,  as  the  gain  increases,  the 
system  output  responses  become  increasingly  decoupled  and 
dominated  by  the  infinite  root  characteristics.  The 
asymptotic  closed-loop  transfer  function  for  the  continuous 
case  has  the  form 


. 1555^1 


(B-11) 


For  the  discrete  case  the  form  is 


r(X)  =  drag 


(B-12) 


where  the  (i  =  l,  ...,2.)  are  determined  by  the  weighting 
matrix, 

For  certain  irregular  designs  where  the  structure 
of  the  output  vector  creates  a  diagonal  r  matrix,  the  system 
will  exhibit  increasingly  decoupled  behavior  (Chapter  IV) . 

In  other  cases,  the  r  matrix  contains  off-diagonal  terms 
which  prevent  full  output  decoupling  as  the  gain  approaches 
infinity.  In  all  irregular  designs,  the  transmission  zero 
always  appears  as  a  finite  asymptotic  root  in  at  '  jt  one 
position  on  the  diagonal  of  F  and  may  appear  on  the  off- 
diagonal.  This  characteristic  places  an  upper  bound  on  the 
time  responses  of  these  particular  outputs  (13) . 


Regular  Plant 


For  the  system  to  be  classified  as  "regular"  the 
first  Markov  parameter  [CB]  must  have  full  rank.  If  this 
is  true,  the  gain  matrices  can  be  found  from 

Kq  =  [^]~4  (B-13) 

and 

=  a[^]“^E  (B-14) 

where 

a  is  a  constant  which  assigns  the  ratio  of  propor¬ 
tional  to  integral  control 

^  is  the  diagonal  weighting  matrix 

The  diagonal  weighting  matrix,  Z  =  diag  } , 

is  specified  by  the  designer.  Each  (i  =  l,  ...,)l) 
determines  the  weighting  of  the  effect  of  a  particular  error 
signal  on  each  control  input.  This  is  the  methodology 
used  in  the  MULTI  design  program  and  is  a  simplified  ver¬ 
sion  of  the  complete  Porter  method.  In  theory,  the  total 
number  of  finite  (slow)  roots  of  the  system  is  equal  to: 

=  n  +  q2.  -  2.  (B-15) 

which  also  equals  +  Z^  (Table  B.l). 

The  roots,  equal  to  (q2.)  in  number,  are  assigned 
by  the  relationship  between  the  proportional  and  integral 
matrices.  If  the  matrices  differ  by  a  simple 


Proportionality  constant,  a,  then  all  of  the  roots  are 
assigned,  under  conditions  of  infinite  gain  (asymptotically) 
to  the  value  of  -a  in  the  s-plane.  By  replacing  a  with  a 
diagonal  matrix,  these  roots  can  be  individually  assigned 
as  the  negative  value  of  its  diagonal  elements. 

Irregular  Plant 

If  the  first  Markov  parameter  [CB]  is  rank  defi¬ 
cient,  then  the  plant  is  called  "irregular."  In  this  case, 
the  C  matrix  must  be  replaced  by 

Z  =  1  Z2^  (B-16) 

where 

ll  =  (B-17) 

Z2  =  ^^2  —12^  (B-18) 

The  matrix  M  in  the  above  equations  is  a  measurement  matrix 
which  is  chosen  such  that  the  matrix  [FB]  has  full  rank. 

The  designer  chooses  the  measurement  matrix  so  that  it  is 
as  sparse  as  possible,  thus  the  smallest  number  of  addi¬ 
tional  measurements  are  required.  Reference  18  gives  an 
approach  for  selecting  the  measurement  matrix  to  achieve 
optimal  decoupling.  Once  M  is  formed,  Kq  and  are  com¬ 
puted  by 


=  a[FB]  I  (B-20 

which  are  similar  to  Equations  (B-13)  and  (B-14) .  As  in 
the  regular  design  case,  the  same  conditions  of  Z-^  root 
assignment  apply  here. 

For  irregular  plants  the  error  vector  e  is  defined 
as 

e  =  V  -  w  {B-21 

where 

w  =  y  +  (B-22 

For  step  inputs  the  values  of  the  rates,  x^,  become  zero 
in  the  steady  state  because  they  represent  kinematic 
variables  (no  B  matrix  entries) . 

The  computer  program  MULTI  greatly  reduces  the 
time  required  to  achieve  a  satisfactory  design.  The  MULTI 
User's  Manual  (  9  )  describes  the  program  and  its  operation 


Appendix  C:  Actuator  and  Sensor  Data 


Introduction 

This  appendix  contains  the  data  provided  by  MCAIR 
for  modeling  the  actuator  and  sensor  dynsunics  for  the 
STOL/F-15  (10) ,  and  the  simplifying  assumptions  made  in  the 
course  of  the  thesis.  Also  the  data  used  for  the  sensor 
noise  model  is  presented. 

Actuator  Models 

The  following  linear  actuator  transfer  functions 
are  provided  by  MCAIR  as  models  of  the  actuator  dynamics: 

Stabilator /Canard 

_ (30.62)  (272.9)^ _ 

(s  +  30.62)  [s^  +  2  (.502)  (272.9) s  +  (272.9)^] 

Aileron/ Flaperon 

_ (40.37)  (322.2)^ _ 

(s  +  40.37)[s^  +  2  (.371)  (322.3) s  +  322.3)^] 

Rotating  Vanes 


39 

s  +  89 


With  the  exception  of  the  rotating  vanes  (C-3)  the  actu¬ 
ator  dynamics  are  third  order  in  "s."  The  current  version 
of  MULTI  has  provisions  for  no  more  than  two  actuator 
states  and  therefore  Equations  (C-1)  and  (C-2)  must  be 
reduced  in  to  second  order  approximation  of  these  equations 
The  approximations  used  retain  the  slowest  root  and  model 
the  complex  pair  as  a  single  real  root  located  at  the 
natural  frequency  of  the  original  poles.  This  results  in 
the  following  second  order  transfer  functions: 

Stabi la tor /Canard  (Approx.) 


(  < 


(30.62)  (272.9) 

(s  +  30.62) (s  +  272.9) 


(C-4) 


Aileron/Flaperon  (Approx.) 


(40.37)  (322.3) 

(s  +  40.37)  (s  +  322.3) 


(C-5) 


The  validity  of  these  approximations  is  demonstrated  by 
comparing  the  low  frequency  characteristics  and  step 
responses  of  corresponding  third  order  transfer  functions. 
It  is  apparent  in  Figures  C.l  through  C.4  that  the  dynamics 
are  reasonably  well  modeled  by  second  order  transfer  func¬ 


tions  . 
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Combination  of  Surfaces 

As  discussed  in  Chapter  III/  it  is  necessary  to 
combine  mathematically  the  ailerons,  flaps,  and  canards 
into  one  equivalent  surface.  The  method  used  is  to  weigh 
the  contributions  of  each  of  the  surfaces  (B  matrix  entries) 
according  to  their  respective  deflection  limits  and  summing 
their  effects  to  obtain  a  single  column  in  the  B  matrix. 
Since  the  deflection  limits  are  not  symmetric  about  the 
equilibrium  position,  the  negative  limit  is  used  for  the 
weighting  calculations.  Following  is  an  example  of  the 
calculations  for  flight  condition  1: 


B  Matrix 

Element 

Deflection 

Limit  (-) 

Canard 

u/S^  .  . 

-2.45679 

-30 

deg 

Aileron 

u/6u  = 

2.192 

-15 

deg 

Flaperon 

u/6^  = 

1.486 

-5 

deg 

u/6^  =  -2.45679  +  (15/30)  (2 .192)  +  (5/30 )  (1 . 486 )  =  -1.1142 

(C-6) 

This  weighting  is  performed  for  each  of  the  affected  ele¬ 
ments  in  the  B  matrix  for  flight  condition  1.  For  simpli¬ 
city,  the  same  relationships  are  used  for  all  of  the  flight 
conditions.  The  data  presented  in  Appendix  D  contains  the 
resulting  B  matrices  for  each  flight  condition.  There  is 
no  accurate  method  of  combining  actuator  dynamics  since 
the  output  signals  are  not  electrical  and  are  physically 
separated.  As  an  approximation,  the  slowest  actuator  (that 


of  the  canard)  is  used  for  the  combined  canard,  flap,  and 
aileron  surface. 

Sensor  Dynamics 

The  sensor  dynamics  provided  by  MCAIR  include  the 
dynamics  of  measurement  as  well  as  an  aliasing  filter. 

Airspeed  Sensor 


_ 1200 _ 

(s  +  30) (s  +  40) 


(C-7) 


nr  y ’T 


i,' 


( 

-  a 

-  .  .  •  ■ 


Angle  of  Attack  Sensor 

_ (14) (209)^ _ 

(s  +  14) [s^  +  2(.74)(209)s  +  (209)^] 


(C-8) 


In  MULTI,  sensor  dynamics  are  also  restricted  to  second 
order  transfer  functions  so  the  angle  of  attack  sensor  model 
is  reduced  by  the  same  procedure  used  for  the  actuators. 

The  resulting  transfer  function  is  given  by  Equation  (C-9) 
and  Figures  C.5  and  C.6  demonstrate  the  validity  of  the 
approximation . 

Angle  of  Attack  Sensor  (Approx. ) 


.  j  .j 

--.A 


(14)  (209) 

(s  +  14) (s  +  209) 


(C-9) 


4 


Fig.  C.5.  Second  and  Third  Order  Angle  of  Attack 
Sensor  Step  Response 


Also,  since  no  data  is  available  for  the  flight  path  angle 
sensor  dynamics,  the  transfer  function  is  arbitrarily 
chosen  to  include  the  slowest  root  of  the  angle  of  attack 
sensor  and  the  pitch  rate  sensor  dynamics.  The  resulting 
transfer  function  is  shown  in  Equation  (C-10)  • 


Flight  Path  Sensor 


(14) (80) 

(s  +  14) (s  +  80) 


(C-10) 


Sensor  Noise 

The  noise  associated  with  measuring  the  outputs  is 
assiamed  to  be  zero  mean,  independent,  white,  gaussian  noise 
injected  into  each  of  the  measurement  channels,  including 
the  minor  loop  pitch  rate  feedback  measurement.  Realistic 
values  for  the  standard  deviations  for  the  various  measure¬ 
ments  are  obtained  from  Mosley  (12)  and  are  presented  in 
Table  C . 1 . 


TABLE  C.l 


REALISTIC  NOISE  VALUES 


Measured  Quantity 


Standard  Deviation 


0 


tv: 


Appendix  D:  Aircraft  Data 


Introduction 

This  appendix  contains  the  original  data  obtained 
from  MCAIR  (10) ,  a  brief  discussion  of  the  conversion  to 
linearized  body  axis  state  space  form,  lists  the  longi¬ 
tudinal  state  space  models  obtained  for  each  flight  condi¬ 
tion  and  finally  the  resulting  open  loop  transfer  functions. 

MCAIR  Data 

The  following  computer  products  (Tables  D.l  through 
D.6)  are  listings  of  nondimensional  body  axis  stability 
derivatives  provided  by  MCAIR  and  the  dimensional  equivalent 
of  that  data.  The  variables  are  defined  as  follows: 


CZA  =  C 

z 

a 

CMA  =  C„ 
m^ 
a 

CXA  =  C 

X 

a 

CZQ  =  C 

CMQ  =  C 

"‘q 

CXQ  =  C 

X 

q 

CZH  =  C 

CMH  =  C 

CXH  =  c 

^h  "'h  ^h 


CZDn=  C 


CMDn=  C 


CXDn=  C 


TABLE  D.l 

_  FLIGHT  CONDITION  1  AERO  DATA 

♦  ♦•»****-»***-»*******-»***********'»-»K-*****-»***-*iHHm-)*-**'«' 

AIRCRAFT  PARAMETERS 

Q  (DYNAMIC  PRESSURE  -  LBS/FT**2)  =  43.0600 

S  (WINS  REFERENCE  AREA  -  FT**2)  =  608.000 

C  (WING  MEAN  AERODYNAMIC  CORD  -  FT)  =  15.9399 

B  (WING  SPAN  -  FT)  =  42.7000 

VT  (TRIM  VELOCITY  -  FT /SEC)  =  200.000 

THETA  =  11. 5840 

W  (WEIGHT  -  LBS)  =  335/6.0 

I XX  (GLUG-FT**2)  =  23634.0 

IVY  (SLUG-FT**2)  181837. 

IZZ  (SLUG-FT**2>  =  199674. 

IXZ  (SLUQ-FT»*2)  =  -3086.00 

*♦■***■****•)(■•»•■)»■■»****■)»■■*■)»■■»******•**  ■*-)»■*♦  •»!••» *■»  *-*■«  ***•)<■**■**** ** 
*•»•**■»•»•»  »■»****  »•»■■»*  »-******1*-****  *  ****#♦»»■)•■»**♦*  »******* 
ALPHA  =  1 1 . 5840 


LONGITUDINAL  NON-DIM  BODY  AXIS  COEFFICIENTS ( 1 /DEG) 


CZA  = 

-.796233E-01 

CMA  = 

.931356E-02 

CXA  = 

. 208763E-03 

CZQ  --= 

0. 

CMQ  = 

-. 169490 

CXQ  = 

0. 

CZH  = 

181866E-04 

CMH  = 

38859 lE-04 

CXH  = 

.681123E-03 

CZU  = 

.6584188-02 

CMU  = 

-. 140683E-01 

CXU  = 

-.246580 

CZDl  = 

257164E-02 

CMDl 

= 

. 572880E-02 

CXDl 

= 

-. 153014E-02 

CZD2  = 

-.955232E-02 

CMD2 

= 

-. 107546E-01 

CXD2 

= 

-.201656E-02 

CZD3  = 

4884278-02 

CMD3 

= 

. 1 12899E-02 

CXD3 

= 

. 136541E-02 

CZD4  = 

451559E-02 

CMD4 

=S 

-.21421 lE-02 

CXD4 

= 

. 925632E-03 

CZD5  = 

.  13502QE-02 

CMD5 

= 

. 129075E-02 

CXD5 

=: 

-.340353E-02 

CZD6  = 

-. 135028E-02 

CMD6 

-. 137616E-02 

CXD6 

-. 340353E-02 

CZD7  = 

. 135028E-U2 

CMD7 

i;; 

. 129075E-02 

CXD7 

= 

-. 340353E-02 

CZDS  = 

1 350288-02 

CMDS 

~. 137616E-02 

CXD8 

-. 340353E-02 

*♦♦♦■»*■»*■*******  »****■*■»■*****  »****•#■«■*■)<•******■*■*****■*•**■»•■»»■ 
*  it  •»***■»•**♦*»»***■»**  HP*"****#****  ***■»»**♦*■«•***■»*•»****•»** 

LONGITUDINAL  AXIS  DIMENSIONAL  DERIVATIVES 


BODY  AXIS  (1/RAD) 


ZA  - 

-127.843 

MA  = 

1 . 36688 

XA  = 

. 335190 

IQ  = 

0. 

MQ  - 

-.991250 

XQ  = 

0. 

ZH  = 

-. 2548218-05 

MH  = 

-. 497684E-06 

XH  = 

. 954355E-04 

ZU  = 

. 184508E-02 

MU  = 

- . 360356E-03 

XU  = 

-. 69099 lE-01 

ZDl 

= 

-4. 12902 

MDl  == 

. 840770 

XDl  = 

-2. 45679 

ZD2 

-15.3372 

MD2  = 

- 1 . 57837 

XD2  = 

-3.23779 

ZD3 

= 

-7.84218 

MD3  = 

. 165693 

XD3  =■■ 

2. 19230 

ZD4 

= 

-7 . 25022 

MD4  = 

-.31 4380 

XD4  = 

1 . 48619 

ZD5 

= 

2. 16801 

MD5 

. 189433 

XD5  = 

-5. 46470 

ZD6 

-2. 16801 

MD6  == 

- . 20 1 968 

XD6  = 

-5. 46470 

ZD7 

= 

2 . 1 680 1 

MD7  = 

. 189433 

XD7  = 

-  5. 46470 

ZD8 

= 

-2. 16801 

MD8  = 

- . 20 1 968 

XD8  = 

-5. 46470 

♦  ***♦*■»*  *■»**»  •»*■»*****»*  «■*********-)(■  ■»******■)<•*■«■**■<(■■«**** 


TABLE  D.2 

FLIGHT  CONDITION  2  AERO  DATA 
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TABLE  D.3 

_  FLIGHT  CONDITION  3  AERO  DATA 

********■»**•*********■»**■»*****■»*•»*************•»****** 

AIRCRAFT  PARAMETERS 

0  (DYNAMIC  PRESSURE  -  LBS/FT**2)  =  48.8700 

S  (WINS  REFERENCE  AREA  -  FT**2)  =  608.000 

C  (WINS  MEAN  AERODYNAMIC  CORD  -  FT)  =  15.9390 

B  (WING  SPAN  -  FT)  =  42.7000 

VT  (TRIM  VELOCITY  -  FT/SEC)  201.100 

THETA  =  15.4400 

W  (WEIGHT  -  LBS)  .=  4351 1."-!;* 


IXX 

(SLUG-FT**2)  = 

35215 

.  0 

lYY 

(SLUG-FT**2)  = 

1 90800 . 

IZZ 

(SLUG-FT**2)  -- 

219105. 

IXZ 

(SLUG-FT-»*2)  = 

-2881 . 

00 

LONGITUDINAL  NON 

-DIM  BODY  AXIS  COEFFICIENTS ( 1 /DEG) 

CZA  = 

-. 784900E-01 

CMA  = 

. 957400E-02 

CXA  = 

. 1 50900E-02 

CZQ  = 

0. 

CMO  = 

-. 169000 

CXQ  = 

0. 

CZH  = 

— . 167600E— 03 

CMH  = 

. 176600E-03 

CXH  = 

. 666200E-03 

CZU  = 

- 1 . 06500 

CMU  = 

. 639400E-01 

CXU  = 

-.619300E-02 

CZDl  = 

-. 263600E-02 

CMDl 

=  . 557600E-02 

CXDl 

= 

-. 155200E-02 

CZD2  = 

-.8315008-02 

CMD2 

=  -. 102000E-01 

CXD2 

-. 274900E-03 

CZD3  = 

-.  559100E-02 

CMD3 

=  .8521008-03 

CXD3 

= 

.  1 15700E-02 

CZD4 

- . 450a00E-02 

CMD4 

=  -.2111 OOE-02 

CXD4 

■c= 

.942100E-03 

CZD5  = 

. 189600E-02 

CMD5 

=  . 255400E-02 

CXD5 

=s 

-. 312000E-02 

CZD6  = 

-.  742200E--03 

CMD6 

=  -.130100E-02 

CXD6 

=: 

-. 359500E-02 

CZD7  = 

.  109600E-02 

CMD7 

=  . 255400E-02 

CXD7 

=: 

-.31 2000E-02 

C2D8  = 

-. 742200E-03 

CMD8 

=  -.  130100E--02 

CXD8 

-. 359500E-02 

LONGITUDINAL  AXIS 

DIMEN 

SIONAL  DERIVATIVES 

BODY  AXIS  (1/RAD) 

ZA  = 

-98.8871 

MA  ^ 

1.36158 

XA  =■ 

1.90114 

ZQ  = 

0. 

MQ  -■= 

-.952482 

XQ  = 

0. 

ZH  = 

-. 1B3259E-04 

MH  = 

.217975E-05 

X.H  = 

. 728443E-04 

ZU  = 

- . 232900 

MU  = 

. 157841E-02 

XU  = 

-.  135432E-02 

ZDl  ^ 

-3. 32101 

MDl  = 

. 793002 

XDl 

-1 . 95532 

ZD2  = 

-10.4758 

MD2  = 

-  1 . 4506 1 

XD2 

-.346338 

ZD3  = 

-7.04393 

MD3  = 

. 121183 

XD3 

zz 

1 . 45767 

ZD4  = 

-5. 67949 

MD4  = 

- . 300220 

XD4 

1 . 18692 

ZD5  = 

2. 38871 

MD5  = 

. 363222 

XD5 

= 

-3. 93079 

ZD6  = 

-.935075 

MD6  = 

-. 185024 

XD6 

-4. 52923 

ZD7  = 

2. 38871 

MD7  ^ 

-7-  y  rf  ^ 

V  O 

XD7 

-3.93079 

ZD8  = 

-.935075 

MD8  = 

-. 185024 

XD8 

= 

-4.52923 

*****■»■)»■•»■-*■(»■  it-** **■»*■***♦***■***•»*  It- ***•»******•»•*•****■»*♦** * 


TABLE  D.4 


FLIGHT  CONDITION  4  AERO  DATA 

**********-*********'****************-)HHHH»-************ 

AIRCRAFT  PARAMETERS 

Q  (DYNAMIC  PRESSURE  -  LBS/FT**2)  =  109.970 

S  (WING  REFERENCE  AREA  -  FT**2)  =  608.000 

C  (WING  MEAN  AERODYNAMIC  CORD  -  FT)  =  15.9390 

B  (WING  SPAN  -  FT)  =  42.7000 

VT  (TRIM  VELOCITY  -  FT /SEC)  =  304.000 

THETA  =  5.04000 

W  (WEIGHT  -  LBS>  =  43511.0 

I XX  (SLUG-FT-»»2)  =  35215.0 

lYY  (SLUG-FT**2)  =  190800. 

Ill  (SLUG-FT**2)  =  219105. 

IXZ  (SLUG-FT**2)  =  -2881.00 

♦  •d-****-********************************************** 


LONGITUDINAL  NON-DIM  BODY  AXIS  COEFFICIENTS ( 1 /DEG) 


CZA  = 

-. 690500E-01 

CMA  = 

. 733300E-02 

CXA  = 

. 854200E-02 

CZQ  = 

0. 

CMQ  = 

-. 157000 

CXQ  = 

0. 

CZH  = 

-.  106400E-03 

CMH  = 

. 1 49500E-03 

CXH  = 

. 290600E-03 

CZU  = 

-. 477000 

CMU  = 

. 358200E-0 1 

CXU  = 

-. 354400E-01 

CZDl 

234500E-02 

CMDl 

.562100E-02 

CXDl 

= 

- . 692400E-03 

CZD2 

= 

-.792100E-02 

CMD2 

= 

-.887400E-02 

CXD2 

= 

434000E-04 

CZD3 

= 

-. 546300E-02 

CMD3 

= 

. 698900E— 03 

CXD3 

= 

. 305B00E-03 

CZD4 

= 

-. 4a2800E-02 

CMD4 

= 

- . 2S8900E-02 

CXD4 

. 328300E-03 

CZD5 

. 824600E-03 

CMD5 

r: 

. 1 1 2400E-02 

CXD5 

- . 1 48600E-02 

CZD6 

= 

-. 1 74800E-03 

CMD6 

= 

409300E-03 

CXD6 

-. 170000E-02 

CZD7 

. 824600E-03 

CMD7 

. 1 1 2400E-02 

CXD7 

-. 148600E-02 

CZD8 

= 

174800E-03 

CMD8 

- . 409300E-03 

CXD8 

=s 

- . 1 70000E-02 

***  ***********  ***************************************** 


LONGITUDINAL  AXIS  DIMENSIONAL  DERIVATIVES 
BODY  AXIS  (1/RAD) 


ZA  = 

-195.759 

MA  = 

2.34674 

XA  = 

24.2168 

ZQ  = 

0. 

MQ  = 

-1.31717 

XQ  = 

0. 

ZH  = 

-.  173182E-04 

MH  = 

. 27468 lE-05 

XH  = 

. 472995E 

ZU  = 

-. 155278 

MU  = 

. 131626E-02 

XU  = 

-. 1 1536BE 

ZDl 

= 

-6. 64814 

MDl  = 

1 . 79886 

XDl 

= 

-1.96297 

ZD2 

= 

-22.4563 

MD2 

-2.83990 

XD2 

= 

-. 123040 

ZD3 

-15. 4878 

MD3  = 

. 223665 

XD3 

. 866951 

ZD4 

= 

-13. 6875 

MD4  = 

-.828543 

XD4 

= 

. 930739 

ZDS 

2.33776 

MD5  = 

. 359707 

XD5 

= 

-4.21285 

ZC  ') 

= 

-.495563 

MD6  = 

- . 1 30986 

XD6 

-4.81955 

ZD7 

2.33776 

MD7 

. 359707 

XD7 

= 

-4.21285 

ZD8 

-.495563 

MD8  = 

-. 130986 

XDB 

= 

-4.81955 

**************************************************** 
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TABLE  D.5 

FLIGHT  CONDITION  5  AERO  DATA  — 

AIRCRAFT  PARAMETERS 

Q  (DYNAMIC  PRESSURE  -  LBS/FT**2)  =  B1.1400 

S  (WING  REFERENCE  AREA  -  FT**2)  =  608.000 

C  (WING  MEAN  AERODYNAMIC  CORD  -  FT)  =  15.9390 


B 

(WING  SPAN  -  FT) 

=  42.7000 

VT 

(TRIM  VELUCITY  - 

FI/SEC)  = 

304 .  <I>00 

THETA 

=  8.07600 

W 

(WEIGHT  -  LBS)  = 

435 11.0 

IXX 

(SLUG-FT**2)  = 

352 15.0 

lYY 

(SLUG-FT**2)  == 

190800. 

IZZ 

(SLUG-FT**2)  = 

219105. 

IXZ 

(SLUG-FT**2)  = 

-2881 . 00 

*****■»****■»•»*•»■*****  **  ****■*♦•»■**■»■*****■»****♦*****  ****  * 

LONGITUDINAL  NON-DIM  BODY 

AXIS  C0EFFICIENTS(1/DEG) 

;zA 

= 

-. 740500E-01 

CMA  = 

859500E-02 

CXA  = 

. 828500E-02 

:zQ 

= 

0. 

CMQ  =  -. 

1 62000 

CXQ  = 

0. 

:ZH 

= 

946a00E-04 

CMH  = 

1 6280t:>E-03 

CXH  = 

. 386600E-03 

:zu 

- . 657000 

CMU  = 

363200E-01 

CXU  = 

-. 122200E-01 

:zDi 

199400E-02 

CMDl  = 

. 586000E-02 

CXDl  = 

108900E-02 

:ZD2 

= 

-. 870000E-02 

CMD2  =  - 

. 955500E-02 

CXD2  = 

-. 452B00E-03 

;ZD3 

-. 525200E-02 

CMD3  = 

. 773900E-03 

CXD3  = 

. 554600E-03 

:ZD4 

= 

-.  464600E-02 

CMD4  =  - 

. 262200E-02 

CXD4  = 

. 507400E-03 

:ZD5 

. 900700E-03 

CMD5  = 

. 122300E-02 

CXD5  = 

157800E-02 

;ZD6 

= 

-.  257300E-03 

CMD6  =  - 

.518700E-03 

CXD6  = 

-. 180700E-02 

;ZD7 

.  9<l)07<l>0£-03 

CMD7  = 

. 122300E-02 

CXD7  = 

-. 157800E-02 

;ZD8 

= 

-.  257300E-03 

CMD8  =  - 

. 518700E-03 

CXD8  = 

-. 180700E-02 

*■)(•*■«■♦****•»*■»*****■)»■*•»**•»*•»**•»*********■*****■»*****•»*■** 

LONGITUDINAL  AXIS  DIMENSIONAL  DERIVATIVES 

BODY 

AXIS  (1/RAD) 

ZA 

= 

-154.897 

MA  = 

2 . 02950 

XA 

=  17.3305 

ZQ 

0. 

MG!  =  - 

1 .  (jC>280 

XQ 

=  <-') . 

ZH 

-. 1 13705E-04 

MH  = 

. 220700E-05 

XH 

=  . 464284E-04 

zu 

= 

157804 

MU  = 

. 984744E-03 

XU 

=  -.293510E-02 

;di 

zr 

-4. 17103 

MDl  = 

1 . 38370 

XDl 

=  -2.27796 

;d2 

= 

-18. 1986 

MD2  =  - 

2. 25618 

XD2 

=  -.947162 

;d3 

=: 

- 1 0 . 986 1 

MD3  = 

. 182738 

XD3 

=  1.16011 

:d4 

= 

-9.71845 

MD4  =  - 

.619123 

XD4 

=  1.06137 

;d5 

=: 

1 .  B8407 

MD5  = 

. 288782 

XD5 

=  -3.30084 

:d6 

= 

-.538217 

MD6  ^  - 

. 122479 

XD6 

=  -3.77986 

:d7 

= 

1 . 88407 

MD7  ^ 

. 288782 

XU7 

-3.3<j<I>a4 

:db 

:= 

-.538217 

MD8  =  - 

. 122479 

XD8 

^  -3.77986 

****■*«■■*•)•■♦*-»•*■**»-»■  x-******-***-******-*****-***************  ■  . 
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TABLE  D.6 

FLIGHT  CONDITION  6  AERO  DATA 


*****■»*****•*■»■»****■»*•»*****************•»•»*■**•»•***■)»■**** 

AIRCRAFT  PARAMETERS 

Q  (DYNAMIC  PRESSURE  -  LBS/FT**2)  =  81.1400 

S  (WING  REFERENCE  AREA  -  FT**2)  =  60B.000 

C  (WING  MEAN  AERODYNAMIC  CORD  -  FT)  =  15.9390 

B  (WING  SPAN  -  FT)  =  42.7000 

VT  (TRIM  VELOCITY  -  FT /SEC)  =  304.000 

THETA  =  6.23200 

W  (WEIGHT  -  LES)  =  33576.0 

I XX  (SLUG-FT**2)  =  23644.0 

lYY  (SLUG-FT**2)  =  181847. 

IZZ  (SLUG-FT**2)  =  199674. 

IXZ  (SLUG-FT**2)  =  -3086.00 

**■****•»*■*■****♦**■**********•»-»■-»*******■)••*  ****  ■****■»**■*■*  * 


LONGITUDINAL  NON-DIM  BODY  AXIS  COEFFICIENTS ( 1 /DEG) 


CZA 

= 

-. 740500E-01 

CMA  = 

. 859500E-02 

CXA  = 

. 828500E- 

CZQ 

= 

0. 

CMQ  = 

-. 162000 

CXQ  = 

0. 

CZH 

= 

-. 946800E-04 

CMH  = 

. 162800E-03 

CXH  = 

. 386600E— 

CZU 

= 

-. 657000 

CMU  = 

. 363200E-01 

CXU  = 

- . 1 22200E- 

CZDl 

-.  199400E-02 

CMDl  = 

. 5a6000E-02 

CXDl 

= 

-. 108900E 

C2D2 

= 

-. S70000E-02 

CMD2  = 

-.955500E-02 

CXD2 

= 

452800E 

CZD3 

= 

-. 525200E-02 

CMD3  = 

. 7  73900E-03 

CXD3 

= 

. 554600E 

CZD4 

= 

-. 464600E-02 

CMD4  = 

- . 262200E-02 

CXD4 

=E 

. 507400E 

CZD5 

= 

. 900700E-03 

CMD5  = 

. 122300E-02 

CXD5 

-. 157800E 

CZD6 

= 

-. 257300E-03 

CMD6  = 

“. 518700E-03 

CXD6 

= 

-. 180700E 

CZD7 

.  9(;)0700E  -Ci3 

CMD7  = 

. 122300E-02 

CXD7 

— . 157800E 

CZD8 

-.25/3O0E-03 

CMOS  ■- 

~. 518700E-03 

CXD8 

- . 1 80700E 

**-«•■**■»•*■»*■)*•■)*•**♦*■»*■)(••»*■*♦**■»**»*■»*-»■»***•»*•»**■»*♦**  »***♦* 

LONGITUDINAL  AXIS  DIMENSIONAL  DERIVATIVES 

BODY  AXIS  (1/RAD) 

ZA 

-200. 730 

MA  = 

2. 12942 

XA 

22.4585 

ZQ 

0. 

MQ  =  - 

1 . 052 1 8 

XQ 

0. 

ZH 

=i 

-.  147350E-04 

MH  = 

. 231566E-05 

XH 

. 601664E 

ZU 

= 

-.204497 

MU  = 

. 103323E-02 

XU 

-.380358E 

ZDl 

= 

-5. 40522 

MDl  = 

1.45182 

XDl 

-2 . 95200 

ZD2 

= 

-23.5835 

MD2  =  - 

2.36726 

XD2 

=Z 

-1 . 22742 

ZD3 

-14.2368 

MD3  = 

. 191735 

XD3 

= 

1 . 50338 

ZD4 

= 

-12.5941 

MD4  =  - 

. 649604 

XD4 

1 . 37543 

ZDS 

= 

2.44157 

MD5  = 

. 303000 

XD5 

= 

-4.27755 

ZD6 

= 

-.697474 

MD6  =  - 

. 1 28509 

XD6 

-4.89831 

ZD7 

= 

2.44157 

MD7  = 

. 303000 

XD7 

-4.27755 

ZDB 

-.697474 

MD8  =  - 

. 1 28509 

XDB 

-4.89831 

**■*•((■****•»■**  *■♦***  *♦■»*•»*♦*■»**■»  ii-**********  **■»**■**  *■»■■»■**♦ 
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where  n  is  an  integer  corresponding  to  a  particular  control 
surface  as  defined  below: 

1  =  Canard 

2  =  Stabilator 

3  =  Trailing  edge  flaps 

4  =  Ailerons 

5  =  Right  top  vane 

6  =  Right  bottom  vane 

7  =  Left  top  vane 

8  =  Left  bottom  vane 

State  Space  Models 

The  nondimensional  body  axis  stability  derivatives 
are  converted  to  dimensional  state  space  form  first  by 
calculating  the  dimensional  derivatives  and  then  the  state 
space  elements  by  the  equations  contained  in  References  2 
and  10 .  These  calculations  are  performed  by  the  Program 
STOLCAT,  an  adaptation  of  the  Conversion  and  Transformation 
(CAT)  program  written  by  Finley  Barfield  (2) .  The  resulting 
state  space  models  are  contained  in  Tables  D.7  through  D.12. 

Open  Loop  Transfer  Functions 

The  open  loop  LaPlace  transfer  functions  for  the 
state  models  at  each  flight  condition  are  listed  in  Tables 
D.13  through  D.18. 


TABLE  D.7 


FLIGHT  CONDITION  1  STATE  MODEL 


0 

-31.55 

0 

-3.233E-2 


A  (Plant  Matrix) 

0  1 
-.06909  -40.161 

-3.603E-4  -.9913 

9.225E-6  .97963 


0 

-1.1142 

.81 

-.0521 


B  (Input  Matrix) 
0 

-3.238 

-1.57837 

-.0766 


0 

-21.79 

-.02507 

0 


0 

.33519 

1.3668 

-.6392 


C  (Output  Matrix) 


TABLE  D.8 


FLIGHT  CONDITION  2  STATE  MODEL 


A  (Plant  Matrix) 

0 

0 

1 

0 

-30.72 

.0120 

-50.22 

22.20 

0 

2.160E-3 

-1.011 

1.020 

730E-2 

-2.410E-3 

.9540 

-.5249 

B  (Input  Matrix) 

0 

0 

0 

-.2829 

.21149 

-24.049 

.7399 

-1.376 

.  3430 

-.0596 

-.0898 

.0210 

0 

0 

1 


C  (Output  Matrix) 
1  0 

0  0 

0  0 


0 

1 

-1 


9 


X  = 


u 

q 


a 
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TABLE  D.9 


0 

-31.038 

0 

4.263E-2 


FLIGHT  CONDITION  3  STATE  MODEL 


A  (Plant  Matrix) 

0  0 
-1.354E-3  -53.53 

1.578E-3  -9.524 

-1.158E-3  .9639 


0 

1.9011 

1.3615 

-1.492 


B  (Input  Matrix) 

0  0  0 
-1.208  -.3463  -16.92 

.803  -1.4506  .356 

.03868  -.05208  .0144 


C  (Output  Matrix) 


0 

A  (Plant  Matrix) 

0 

1 

0 

-32.07 

-.0115 

-26.70 

24.21 

0 

-1.31E-2 

-1.317 

2.347 

9.305E-2 

-5.107E-4 

.9961 

-.6439 

B  (Input  Matrix) 

0 

0 

0 

-1.517 

-.123 

-18.06 

1.773 

-2.839 

.4562 

-.0612 

-.07386 

.01212 

C  (Output  Matrix) 

q 


TABLE  D.ll 

FLIGHT  CONDITION  5  STATE  MODEL 


A  (Plant  Matrix) 


0 

0 

1 

0 

-31.88 

-2.930E-3 

-42.71 

17.33 

0 

9.847E-4 

-1.003 

2.029 

.01488 

-5.191E-4 

.9901 

-.50953 

B  (Input  Matrix) 

0 

0 

0 

-1.52 

-.9471 

-14.16 

1.51 

-2.256 

-.3300 

-.0323 

-.0597 

8.84E-3 

C  (Output  Matrix) 


[•Viil 


A  (Plant  Matrix) 

0 

0 

1 

0 

-30.88 

1.007E-2 

-56.56 

21.01 

0 

1.167E-3 

-.7558 

.7577 

4.533E-2 

-1.465E-3 

.9591 

-.3034 

B  (Input  Matrix) 

0 

0 

0 

-1.482 

-.4642 

-15.35 

.6499 

-1.227 

0 

-.0245 

-.0679 

9.863E-3 

C  (Output  Matrix) 

q 


s  -  .2865)  (s  +  1.728)  (s  +  .04099  ±  .2763j) 
.2829(s  +  136.9) (s  +  .3181  ±  .4282j)  /  C.E. 
.0596(s  -  10.91) (s  +  .0290  ±  .2795j)  /  C.E. 
0596(s  -  .0625) (s  +  .8106  ±  2.056j)  /  C.E. 
2115(s  +  318.2) (s  +  .3051  ±  .4996j)  /  C.E 
.0898(s  +  15.57) (s  +  .02635  ±  .2777j)  /  C.E 
0898  (s  +  3.005)  (s  -  2.648)  (s  -  0565)  /  C.E 
.74.05(3  +  2.256) (s  -  .0118  ±  .3845j)  /  C.E 
0210(3  +  19.27)  (3  +  .0339  ±  2701j)  /  C.E. 
.0210(3  +  4.26) (3  -  1.402) (s  +  .1499)  /  C.E 


TABLE  D.15 

FLIGHT  CONDITION  3  TRANSFER  FUNCTIONS 
C.E.  =  (s  -  .04311  ±  .2615j)(s  +  2.387)  (s  +  .1451) 
u/6^  =  -1.028  (s  +  42.4)  (s  +  1.30)  (s  +  .629)  /  C.E. 

u/6g  -  -.0386  (s  -  19.11)  (s  -  .00964  ±  .1909j)  /  C.E. 

u/6^  =  .0386  (s  -  .0217)  (s  +  .8471  ±  5.346j)  /  C.E. 

a/6^  =  -.3463  {s  -  223.6)  (s  +  1.494)  (s  +  .6005)  /  C.E. 

a/6g  =  -.0521  (s  +  27.77)  (s  +  .01128  ±  .1939j)  /  C.E. 

a/6^  =  .0521(s  -  6.522)  (s  -  .02176)  (s  +  6.485)  /  C.E. 
Y/6^  =  -16.92(s  +  2.832)  (s  +  .3682  ±  4921j)  /  C.E. 
Y/6g  =  .0144(s  +  26.14)  (s  +  .00127  ±  .1971j)  /  C.E. 
Y/6^  =  -.0144(s  -  5.329)  (s  -  .0502)  (s  +  6.801)  /  C.E. 


TABLE  D.16 


FLIGHT  CONDITION  4  TRANSFER  FUNCTIONS 


TABLE  D.18 

FLIGHT  CONDITION  6  TRANSFER  FUNCTION 

C.E.  =  (s  +  1.413)  (s  -  .2945)  (s  +  .01514  ±  .2130j) 
u/6^  =  -1.482(s  +  25.78) (s  +  .2643  ±  .3803j)  /  C.E. 

u/6g  =  -.0245(s  -  24.81)  (s  +  .0146  ±  .2176j)  /  C.E. 

u/6^  =  .0245(s  -  .0498)  (s  +  .8948  ±  2.866j)  /  C.E. 
a/6^  =  -.4642(s  -  145.9) (s  +  .3198  +  .4050j)  /  C.E. 

a/6g  =  -.0680  (s  +  18.01)  (s  +  .01417  ±  .2181j)  /  C.E. 

a/6^  =  .0680  (s  -  2.678)  (s  -  .0510) (s  +  2.727)  /  C.E. 
Y/6^  =  -15.35(s  +  1.449) (s  -  .1515  ±  .1061j)  /  C.E. 
Y/6g  =  .00986(5  +  3.027)(s-  .000118  ±  .1982j)  /  C.E. 
Y/6^  =  -.00986  (5  -  .3757)  (5  +  2.42)  (5  +  .982)  /  C.E. 
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